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This thesis attempts to study the exchange-correlation 
(EC) effects in atoms using the powerful coupled-cluster 
many-body theory. Many-elec tron EC effects have attracted 
much attention in theoretical chemistry. One of the 
approaches to the study of EC effects is to start from an 
independent particle model and introduce refinements that 
account for the residual inter-electron interactions. The 
Hartree- Fock approximation has traditionally served as the 
independent particle model for most of these studies. In 
this thesis, we use the Slater Xa approximation, which is a 
local potential model, to provide the reference starting 
point for the study of EC effects in atoms. Our aim is to 
assess the effectiveness of this local potential model in 
representing EC effects in atoms, as compared to the Hartre 
Fock model. 

The coupled-cluster (CC) many-body theory of Coester, 
Kummel and Cizek has been used in this thesis. The CC 



methodology relies on the cluster expansion of the wavefunction 
and uses the language of second quantisation, the associated 
Feynman graphology and the powerful linked-cluster theorem. 
Further, the methodology encompasses several sophisticated 
many-elec tron methods like, the coupled electron-pair approxi- 
mation, Nesbet-Bethe-Goldstone approach and a hierarchy of 
approximations to many-body perturbation theory.^ Earlier 
studies have used the Hartree-Fock method to provide the 
reference independent particle approximation for the CC expan- 
sion. Surprisingly, this expansion method has been used 
merely for the purpose of obtaining accurate correlation 
energies although its versatility extends well beyond. In 
fact, the method can be used to advantage for studying the 
independent particle approximations that serve as reference 
points, 

''In this thesis, the CC methodology and its allied appro- 
ximations are used to evaluate EC corrections to the Xa 
potential and to provide an insight into the nature of the 
potential vis-a-vis its effectiveness in representing the 
true EC potential. Empirical and semi-empirical modifications 
to the Xa method are also examined in this light, with a view 
to gauge their efficacy in improving upon the parent independent 
particle approximation. ’ It is hoped that a study such as that 
made in this thesis will provide a rational basis for the 
choice of local potentials that describe the EC effects in 
atoms more accurately. 



The thesis consists of five chapters, a bibliography 
of the references cited in the text of the thesis and two 
appendices . 

Chapter I of the thesis reviews the state of the art 
with reference to the treatment of many-electron interactions 
in atoms. The independent particle approximations normslly 
in vogue and the methods to go beyond these are discussed. 

The linked-cluster expansion which forms the basis of the 
present work is discussed briefly, in the context of pertur- 
bation theory. The chapter concludes with the delineation 
of the scope of the thesis. 

Chapter II introduces the CC methodology and its basis. 
With the aid of Appendix A, this chapter attempts to provide a 
detailed explanation of the Feynman-Goldstone graphology, 
which is used to derive the CC equations. Expressions for the 
correction to the EC energies and equations that define the 
cluster coefficients are explicitly- given. A comparison of 
the CC methodology with various other many-elcctron theories 
is made with a view to highlight its powerfulness. 

Chapter III presents a model calculation on beryllium 
atom to illustrate the application of the CC methodology and 
the approximations derived from it, in the context of evalua- 
tion of corrections to EC energy. Computational details 
are also given. 



Chapter IV deals with the detailed analysis of the 
Xa potential through the EC corrections obtained as described 
in the previous chapter. For the purpose of this analysis, 
the ideal many~body atomic systems, beryllium and Be-like four 
electron ions are studied. The role of the exchange parameter, 
a, in the local potential model and refinements to the Xa 
potential are examined with reference to the EC corrections. 
Attempts are made to obtain physically meaningful insights 
into the nature of the Xa approximation. 

Chapter V reports the expectation values of some 
operators over the wavef unctions obtained through the CC 
expansion for the systems discussed in the previous chapter. 

The results are compared with the reference state expe’etation 
values. 

The thesis concludes with some salient remarks on the 

application of the CC methodology for the study of EC effects 

in atoms, A few interesting features emerge from the present 

analysis of EC effects. Firstly, ‘the CC methodology is seen 

to be extremely versatile and superior to its other many- 

electron counterparts. And, it can be used to assess the merit 

of the independent particle approximations that provide the 

reference states. With specific reference to the Xa method 

the analysis reveals that the Xa local potential describes 

very accurately the inner-shell EC effects, as is evidenced 
2 

by the 1s pair corrections in Be isoelectronic series. On 
the other hand, the valence- shell EC effects are poorly 



described by this local potential. It can be surmised from 
the analysis of EC effects that for obtaining accurate des- 
cription of many-elcc tron interactions in atoms with the Xa 
method, it will suffice to improve upon the valence-shell 
description also, to achieve this correction through Cl-like 
procedures, it will suffice to consider a very limited basis 
set of Xct orbitals. 

Appendix A presents a brief discussion on the language of 
second quantisation, Feynman- Goldston e graphology as applied to 
many-body perturbation theory, fourth-order Raleigh-Schrddinger 
perturbation theory and cancellation of the unlinked diagrams, 
the hole-particle formalism, VJick ‘ s theorem and contractions 
and illustration of the elementary Feynman-Goldstone diagrams 
used in this thesis. 

Appendix B lists some of the computer programs used in 


the present study. 
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1.1 THE SCHRODINGER EQUATION 

The non-relativi Stic Schrodinger equation, 

hU^ = eN> (1.1) 

provides a theoretical framework for the quantum mechanical 

I*. 

description of a many-electron atomic system with H as the 

Jl 

Hamiltonian operator^'^ and E, the energy. The many-clectron 
state function | t|i > in configuration space, obtained as a 
solution of the above equation contains, in principle, all 
possible information about the system. It is a function of 
the space and spin coordinates of all the electrons of the 
system. However, the equation is not exactly solvable except 
for the simplest case of one-electron atomic systems [1 ] 
owing to the presence of the interelectron interaction terms 
in the Hamiltonian. It becomes, therefore, necessary to 
resort to approximations that are mathematically tractable, 
physically meaningful and reasonably accurate. The problem 
simplifies considerably in the case of atoms due to the 
separability of nuclear and electronic motions. The first 

M." ■ ■ I ■ ■ 

’’'^All quantum mechanical operators will be represented with a 
"■*'* on top of the corresponding symbol. 

* 

Numbers in square brackets pertain to literature references 
which appear at the end of this thesis. 
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step in the solution of the many-electron Schrodinger 
equation would be to reduce the N-electron equation to IM 
one-electron equations. This is the leit motiv behind all 
the independent particle approximations. These approximations, 
however, fail to portray correctly the many-body effects that 
originate from simultaneous interaction between many electrons. 
The exact treatment of this many-body problem is complicated. 
Although several sophisticated methods hove been developed 
for solving the problem, the state of the art can be said to 
be still in a nascent stage. An effort is made in this thesis 
to describe the many-body interactions between electrons in 
atoms using the coupled-cluster (CC) method? this appears to 
be the most promising in the many-body theory horizon. 

1,2 THE INDEPENDENT PARTICLE APPROXIMATION (IPA). 

1.2.1 The Hartree and Hartree-Fock Methods 

The underlying premise in all those approximations is 
that every electron moves independent of the other electrons, 
i.e. in an averaged field produced by all other electrons and 
nuclei. The choice of description of the averaged field and 
its proximity to the exact description varies from model to 
model. In the case of an atom, the averaged field can be 
conveniently located at the centre and the ensuing spherical 
symmetry of the potential reduces the complexity of the 

A 

problem considerably. The averaged potential, V(r), then 
depends on the radial part of the wavef unction , R(r), which in 
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A 

turn depends on V(r) and so both \/(r) and R(r) have to be 
determined together in an iterative way. 


The problem of setting up a self-consistent field ( SCF) 
was first outlined by Hartree [2], He assumed that the field 
of all electrons on a particular electron can be approximated 
by the electrostatic action of their averaged charge density. 
The one-electron equations of Hartree can therefore be 
written as [3] 

7^(1) + 0^(1) + v^id)] d(i) = d-2) 


where. 


v^d) 


■2Z 


^AZ) 'J>..(2) 


^ T-X*-/ T- 


d T- 


and 


V,id) = 


-n. /-^ 


^.(2) 


^ 1-^2 


dx , 


(1.3) 


(1.4) 


Atomic units have been used alongwith the usual notation [3]. 

A. 

\J^ is the field arising from the nuclear charge and the total 
charge density of all the electrons in the system and 0^^ 
is the term that corrects the coulombic term for the self- 
interaction of the electron. The decomposition of the N- 
electron function into a simple product of N one-electron 
functions is tacitly implied by the above method. This 
coupled with the application of the variation theorem for 
obtaining the energy extremum yieldsthe one-electron equa- 
tions [4]» that were obtained through intuitive arguments by 
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Hartree. The Hartree methDd fails to depict the Fermion 
characteristic of the many-electron system; also it does 
not account for the instantaneous interaction of the 
electrons, a shortcoming inherent in all the independent 
particle models. The Fermi-Dirac Statistics [5], when 
incorporated into the SCF method through the use of deter- 
minantal wavef unction [6] yields the Hartree~Fock one- 
electron equations [7] 


[- « .( 1 ) 


1 "j/ «*(2) *^(2) 

<3 


1 




^.( 1 ) 

J 


e . 

1 


4)^(1) 


(1.5) 


The last term on the left-hand side (LHS) of the equation 
(1.5) now includes in addition to the self -correcti on term 


* 

- / (2) ^^(2) 


dfo <!>■ ( 1 ) the terms, 


I n. / <J*(2) <!, (2) 
j^i ^ 


r" -r ‘=‘■'2 

^2 ^ 


( 1 . 6 ) 


characterised as exchange terms. Tho only non-vanishing 
contributions in eqn. (1.6) arise from the summation over the 


orbitals <|>j which have tho same spin as owing to the 

orthogonality constraints. The term can be elegantly recast 
in the form of a potential multiplying the orbital 

- i:nj<^*(1) ^*(2) «,(1) ^A2) (1/lr,-rJ) dt 

\! { 4 \ ^ .1,- -...I.-... 1 X C. 

“ «t(i) «i(i) 


(1.7) 
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This form yields greater insight into the above, since 
this is the potential of a charge whose magnitude is one* 
electronic charge removed from a hole surrounding the loca- 
tion of electron 1 [S], This hole, termed as the Fermi hole, 
accounts for the correlation of like-spin electrons but 
cannot account for the correlation of unlike-spin electrons. 
Therefore, the last term on the LH5 of eqn. (1.5) should 
be corrected for the coulombic correlation of unlike-spin 
electrons. 

1.2,2 Local Potential Methods 

The solution of the Hartree-Fock (HF) equations is 
beset by two main problems associated with the exchange- 
correlation (EC) potential. Firstly, the potential is 
different for different orbitals. Secondly, the EC 
potential is non-local in the sense that it cannot be 
written as a simple product of the type V(r) <^{ t ) . In 
fact, it assumes the form \/(r) <l>(r), where 

V(r) =/ dt^ v(r,r') Mr') (I.S) 

These two factors complicate the HF equation enormously and 
any further simplification should focus attention on these 
two problems. The first problem can be overcome by using 
the average value of the potenbial over all the orbitals, so 
that a single EC potential is common to all one-particle 
equations [ B j , The second problem is more complicated; 
concomitant to the non-local nature of the potential is the 
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non-uniform distribution of the electron density. This 
non-uniformity prevents the treatment of the many-olectron 
system as a homogeneous electron gas, a problem that has 
been salved elegantly by Pines [ 9 ] . However, inhomogeneous 
electron gas problems have been solved, albeit approximately, 
using local EC potentials [8,10,11], 

All these methods assign a key role to one-clectron 
density P(r) and utilise a local density potential, thus 
eliminating the need to solve many-dimensional integrals, 

A rigorous theoretical justification for such local— density 
potential methods is provided by the theorem of Hohenberg 
and Kohn [l2]. The theorem states that the ground-state 
energy of an inhomogeneous electron system is a unique 
functional of density p(r). For a given external potential 
0(x) 

E [P(r)] = / v(r) P(r)d^r + F[ p(r)J (1.9) 

This theorem merely proves the exi stence of a unique 
functional and offers no insight into its exact nature. 
However, a number of approximate functionals have been pro- 
posed to suit various types of many-elec tron systems [13,14]. 
The energy functional can be decomposed into various terms 
and eqn, (1,9) rewritten as 

E[p(r)] =/ v(r) p(r)d^r + d^r d^r’ 

+ Tg[p{r) ] + [p(r) ] 


( 1 . 10 ) 
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where in the right-hand side (RHS) the second term corres- 
ponds to the classical coulombic interaction, the third to 
the kinetic energy of the non-interacting many-electron 
system with density p(r} and the last term, to the 

EC energy term. The function E [p(r)] attains its minimum 
value for the exact ground-state density. Exploiting the 
formal similarity between the Hartree equation and eqn. 
(1.10), a set of one-particle equations can be obtained [15] 

\ 7 ^\ = "i 

( 1 . 11 ) 

The Lagrange parameters arise from the charge normali- 

sation constraint. In the above equation the potential 
Vj^^, is the functional derivative of the EC energy 

'^XC 


(r) = 


--^TjTu'ir 


( 1 . 12 ) 


All the many-body effects are included in this EC potential, 
Vxc(^)» which in contrast to the HF potential, is local. The 
construction of potentials Vyj,(r) is central to all local- 
density potential schemes. 


As a first approximation, the many-electron system may 
be regarded as a homogeneous electron gas. The EC potential, 
v^^ in eqn, (1.12) may then be approximated as the exchange- 
only potential of the homogeneous electron gas [10,11,16]. 

The exchange energy of the electron gas [17] 
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E^c [ d^r (1.13) 

gives the exchange potential 

[Ptr)] = -2 [|;- p(r)]'''3 (1.14) 

The total energy aquation (1.10) may then be solved varia- 
tionally to determine the lowest energy [l5,10], A similar 
approximation to had earlier been obtained by Slater [ B ]. 

He replaced the exact HF potential in the one-electron 
equation by the free-electron exchange potential and obtained 
a value. 

'^XC = 47 (1.15) 

This differs from the Gaspar, Kohn and Sham (GKS) [15,18] 
value by a factor of 3/2. This is so because, the order of 
the two steps taken, one variational in nature and the other 
an approximation, differ in the two methods. In the GKS 
method the electron gas exchange approximation is made first 
and then the variation is taken, while in Slater's method [B] 
the variation is taken before the approximation is made. Since, 
the two steps do not commute with each other, the potentials 
differ by a multiplicative constant. It is worth noting that 
in atomic calculations, local EC potentials mentioned above 
cannot be expected to reproduce the HF exchange potentials, 
owing to the non-uniformity of the electron density . 
distribution in atoms. 
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To account for slight non-uniformities in charge 
density, vis-a-vis the Fermi-hole shape, in a semi- 
phenomenological manner, Slater [19], advocated the 
introduction of a parameter a , such that 

"xc = -3“ [j; (1.16) 

There is no minimum principle associated with the dependence 
of the total energy on a and consequently a cannot be 
determined variationally . However, various methods for 
determining a have been postulated and a review of these 
methods is given by Slater [20]. The method, known as the 
Xa method, is non- vari ational to the extent of ot being 
an arbitrary parameter. Also, it only retains the spirit of 
the theorem of Hohenberg and Kohn [12] and cannot be 
deemed to follow directly from the theorem. Nonetheless, 
the Xa method has numerous merits and successes to its 
credit [21 ]. At the theoretical level, the method satisfies 
both the Hellmanj-Feynman [ 22] and the virial theorems [ 23] 
and the energy expression rigorously leads to Fermi-Dirac 
Statistics [2D]. At the computational level, it has decided 
superiority over the HF method, where the ease of computation 
is concerned. Also it reproduces with atleast as much 
accuracy [24], the double-zeta basis atomic HF results of 
Clementi [25], 
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1.3 EXCHANGE AND CORRELATION 


Since quantum chemical literature is replete with 
confusion about the definition of correlation energy and EC 
energy, it is necessary to clearly delineate these terms. 
Correlation energy is defined. [26] as the energy difference 
between the exact non-relativistic energy , , and the HF 

sxngle determinant energy Epp» 


E 

corr 



(1.17) 


The HF energy includes contributions due to the correlation 
between like-spin electrons originating from the Pauli 
exclusion principle. On the other hand, approximations to 
the HF potential, like that defined by eqn, (1.14) cannot 
be strictly termed as ’’exchange" potentials. By virtue of 
these not being variationally determined potentials, they 
include some correlation although a clear-cut demarcation 
between exchange and correlation cannot be made. Any corre- 
ctions to methods using such potentials should be rightly 
termed as EC corrections and not as correlation energy 
corrections, 

1.4 BEYOND THE INDEPENDENT PARTICLE APPROXIMATION 


1.4,1 Drawbacks in IPA 

The independent particle models satisfy admirably, 
most quantum chemical computational needs and are rarely 
in error of greater than 1^, where total energies of atoms 
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and molecules arc concerned* However, by their very 
nature, they suffer from the drawback of approximating the 
instantaneous, inter-electron potential by an averaged one- 
particle potential, thereby inexactly accounting for the 
entire correlations between electrons [26]. This short- 
coming becomes particularly serious where small energy 
differences between different electronic states or different 
geometries (for molecules) are concerned [27], In such 
cases the error in energy duo to electron correlation may 
parallel the energy differences themselves. Also, where 
two-electron properties are concerned, the role of electron 
correlation may be very important [20]. There is a need, 
therofore, to treat the motion and interaction of electrons 
exactly; and the independent particle models which so 
nearly mimic the actual interactions can serve as excellent 
starting points for more exact treatments. 

As mentioned earlier, the HF method takes into account 
a large portion of the dynamic correlation between the 
electrons with like-spins, whereas it neglects entirely 
the correlations between unlike-spin electrons. The lodal 
density potential methods also neglect a part of electron 
correlation although a clear-cut distinction, like the 
neglect of unlike-spins electron correlation, as in th® 
case of HF method cannot be made. In these models the 
correlation errors can be attributed to the approximation 
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of non-uniform electron density distributions by uniform 
density distributions; and hence this warrants a different 
kind of analysis of correlation effects. Besides this, 
there is also a non-dynamical correlation error peculiar 
to the HF-like models. In these models, the set of single 
particle basis states, may give rise to several degenerate 
configurations, and choice of any one configuration alone 
may cause the model to break down [29]. 

The attempts to go beyond the IPA can be broadly 
categorised into two divisions; viz., traditional and 
modern methods. These methods include (a) the method of 
configuration interaction (Cl) involving the expansion of 
the wavefunction in terms of several configurations [30]; 

( b) the electron-pair function treatments [31]; (c) the 
method of incorporating the inter-particle coordinates 
directly into the wavefunction [32]; (d) the classical 
perturbation expansion methods which are derived from the 
Cl method [33]; (e) the many-body perturbation methods [34] 
and (f) other field-theoretic methods like the CC [35] and 
Green's function methods [36]. In the ensuing discussion, 
no efforts are made to compartmentalize these methods and 
due emphasis is laid only on their conceptual development 
and inter-relationship. 

1.4.2 The Cl Method 

The method of configuration interaction [37] owes its 
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origin to two well established principles; one, the 
expansion of the configuration function in terms of a 
complete or near complete basis of one-electron functions 
[32] and the other, the Ritz method of linear variations 
[38]. Conceptually, this implies the expansion of the 
cusp-like behaviour of the exact wavefunction in a Fourier 
series. The method employs the linear expansion ansatz for 
the trial wavefunction 

N 

= I Cy (1.18) 

v=1 

where are predetermined expansion functions and 

are the variational parameters cons brained to minimise the 

energy functional 

E (1.19) 

A formally exact solution to the N-particle problem 
is obtained, if and only, if all the N-tuple excitations 
are considered [39] , i.e. 'i' is chosen such that 


^ = ^ * I * I I c-."" + ... 

ar ci<6r<s ®pcip 

All the O' s are orthonormal Slater determinants. 


( 1 . 20 ) 


denotes the Slater determinant obtained by replacing the 
orbital a by r in 0^; 0^^ p denotes the Slater deter- 

minant obtained by replacing a and 8 by r and s in 0^, 
and so on. In general, we shall use the label a , 6 ... 
for occupied orbibals and r, s, ... for unoccupied orbitals. 
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Indices i,j ... will be used to represent either type of 
orbitals. 


Inserting eqn. (1.18) into H'? = E H' and left-multiplying 


by yields 


I = EC^ 


( 1 . 21 ) 


where 


Hvu = » 


<S) > 


When Op = it immediately fallows from eqn, (1.21) 


that, 


<0jHh>C <ojti|o^>rJ + I I < \lHl 4a 

° aa„_u a a<6r<s'^ 


or 


= 'EC 


This equation can be rewritten as, 


E * Eq + I e + I r 

° a ^ a< 


( 1 . 22 ) 


(1.23) 


A 

with E * H 0 > 

m 


0 r 

o r < s 


The coefficients are determined by projecting the Hamiltonian 
with various excited states 4>y . 

With 0^= oj 

° ° sP® ^ 'p<fi s<u 


V V ^ ut aSUW . pSUW ppT 


P<5<e s<u<w 


P6e 


o 


(1.24) 
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with , 
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u<w 
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uw 
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+ 1 
6<e<n 



.uwx ^ _uwx 
®6£n^ 


+ I I 

6<e<n<£ u<w<x<y 


rs 


< 'I* 
a(3 



uwxy 


pUwxy 

^ ^ 5 e n S 


EC 


rs 
a 3 


(1.25) 


This projection can be continued till <5^ is an n-fold 
substituted determinant. All the coefficients are then 
obtained and the hierarchy of equations yields the upper 
bound of the exact energy. The coefficients are coupled in 
a complicated manner such that, even though , 0^^ and ^ 
determine the energy, a knowledge of requires and 

*-035 that of Cpjg requires and The solution 

of the entire hierarchy is computationally unmanageable 
and this necessitates the truncation of the wavefunction 
expansion and concomitantly the hierarchy of equations. 

A reasonable truncation point is the set of double excited 
determinants and this leads to various electron-pair 
theories [31,40,41], This truncation, however, voids the 
separability criterion for energy. The energy of an 
N-particle system can be decomposed into contributions 
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arising from one-particle terms, simultaneous interaction 
of two-particles, simultaneous interaction of three parti- 
cles, etc, [42]. This decomposition in terms of energies of 
clusters of one, two, three, etc. particles is unique. All 
these clusters must be linked in the sense, that an N-body 
cluster is not decomposable into (N-r) and i body clusters, 
e.g., the energy arising from simultaneous interaction of 

4-particles cannot be written down as a product of two 2- 

s 

particle interactions. This is nBces|itated by the fact 
that energy, an extensive property, depends only on N, while 
the contribution from unlinked clusters would be proportional 
to higher powers of N [43]. While in complete Cl this 
linkedness is automatically taken into account, truncated . 
Cl incorporates several of the spurious unlinked cluster 
terms, thereby causing a breakdown of the variational energy 
extremum principle [44], This serious problem, termed size- 
extensivity problem common to perturbation theories as 
well, is highlighted with a view to emphasise on the unphy- 
sical nature of basis truncations. The magnitude of this 
problem becomes transparent in fourth-order perturbation 
theory, which is discussed later in Appendix A. Besides 



10 


the size-Bxtensivity problems, the Cl method suffers 
from other computational handicaps. Firstly, the number 
of reference functions needed to obtain meaningFul results 
is enormous. The complexity is apparent from the fact 
that the number of symmetry-adapted Cl configurations for 
a given level of excitation is approximately 
j^{ n N/2 ) ^/( 1 ♦ ) ^sj where 1 is the level of excitation, n 
is the number of basis functions, N is the number of 
electrons and s is a symmetry factor that should be between 
2 and 10. For a simple case of a lO-electron system with 
excitations upto 1=4, this amounts to about 10^ confi- 
gurationsl Besides, the convergence of the method with 
increase in the level of excitations is extremely slow. 

1.4.3 Electron-Pair Methods 

The Cl method can also be modified to deal with, 
instead of an N-electron problem, an N(N-1)/2 electron-pair 
problem [31,45]. A justification for this reduction origi- 
nates from two reasons : ( i) the Pauli Principle prevents 



more than two electrons from occupying the same point in 
space and (ii) the many-electron Hamiltonian consists 
of no more than two-particlc operators. In view of this, 
the correlation problem could be solved independently for 
each pair of particles in a many-particle system [ 46] , The 
correlation energy obtained in this case would be formally 
equivalent to that obtained with the approximation, 


'P 


I 

a3 


I 

rs 



rs 
a 3 


(I.26) 


!C S 

in the Cl expression. Here represents the excitations 

of a pair of electrons in orbitals ct , 3 to orbitals r,s 

in the reference state $ . The Cl coefficients are obtained 

o 

by a modified Schro'dinger equation for each pair (a 3), 

This equation for a given pair of electrons [ 40] takes the 
form ; 


(1.27) 

where the projection operator ) annihilates all com- 

ponents containing occupied orbitals , j a » 3 from the 

orbital expansion 7 (1,2), an arbitrary two-particle 

0.0 

function. The solutions of these sets of equations, referred 
to as the Bethe-Goldstone equations, corresponds to the 
exact solution of a two-particle problem, subject to the 
constraint that the wavefunction should be orthogonal to the 
remaining (N-2) orbitals of an N-particle Slater determinant. 
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The _ Method 

Instead of trying to expand the cusp-like behaviour 
of the wavefunction in a Fourier-like series, as in the Cl 
method, it is possible to incorporate it directly into the 
wavefunction. Slater [47] suggested the use of a factor 
B in the wavefunction to obtain the desired effect, 

Hylleraas [ 32] used an approximation to the above. For the 
ground-state of helium, he used the expression 


T( , r 2 ) = e 


-(r^+r2)/2 “ 


n, l,m=o 


\ ^ / \ 21 m 

^n,2l,m^^1’^^2^ ^12 


(1.28) 


where n,l,m are the set of quantum numbers, r^ and r^ are 
the electronic coordinates and r^ 2 in ter-electionic 

distances. With a mere six-term expansion, Hylleraas 
obtained a near-exact energy for helium atom. Many attempts 
[48,49] have been made to extend the ' r^ * method 

to many- electron systems, by multiplying the dcterminantal 
wavefunction by a correlation factor containing all the r. .*s. 
However, these methods based on explicit use of r. coordinates 

X J 

are very difficult both to implement [ 50] and to interpret. 

An interesting combination of Cl and Hylleraas r^ 2 niethod has 
been applied to atoms with a fair amount of success [50,51]. 

In this method the configurations in the wavefunction are 
chosen to be antisymmetrised, projected products of one- 
electron functions with powers of inter-electronic 
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coordinates. Recently [52] the above method has been 
elegantly recast in a form suitable for the deployment of 
the powerful linkcd-cluster expansion of many-body theory. 

1.4,5 Perturbation Methods 

Both the Cl and the Hylleraas method are beset with 
problems associated with large computations, for in both 
cases the number of parameters to be variationally deter- 
mined are large. A means of reducing the computational 
effort, is to choose carefully limited and meaningful 
terms in the variational procedure. An insight into this 
is achieved by partitioning of the model space into appro- 
priate subspaces [ 53] , Adaptations of the variational form 
then, leads in a straightforward manner to various forms 
of perturbation theory. An arbitrary partitioning of the 
model space into a single function and a remainder X , 

orthogonal to yields 

E = < $^|h1 h1 X><x1 E-H !x>‘"'*<x1H I > (1.29) 

■A A 

If H is partitioned into a model Hamiltonian and a per- 

A 

turbation V 

H = H + V (1.30) 

o 

H 5 = E $ 

O O 0 0 


(1.31 ) 
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then 

E = E^+E^+ <O^101x><x1E-Hq-V|x> <X Ivl > (1.32) 

Expressions for various forms of perturoation theory can 
be obtained from eqn, (1.32) by using the operator expansion 

(A - a”’' + A’’''b + ... (1.33) 


for the inverse term in the above equation. 
When 


A = <xiE - jx> 

B = <X |vl X > 


(1.34) 

(1.35) 


the Brillouin-Wigner (BW) form of perturbation theory is 
obtained [54], The Raleigh-Schrodinger form of perturbation 
theory (RSPT) is obtained when 


A = <X 1 > (1.36) 

B =<xl^ ~ AeIx> (1.37) 


where 


4 E 
A E 


is defined by the equation 

- E - E„ 

0 


CO 



E 

n 


(1,38) 


On expanding and regrouping the terms, 


AE 


I <«o *^1^0 ^ <I*39) 

k=o 
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where P is the projector for the orthogonal complement to 
and defines the partitioned subspace 

P = lx><xlx>~^ <xl ( 1 . 40 ) 


Truncation of the summation at various k's leads to different 
orders of perturbation theory and the order-by-order corre- 
ctions can be obtained. For any order m. 


E . = < 0 

m+1 o 


$ > 
m 


( 1.41 ) 


? = lx><xi E-S ix>“‘' <xlC(I-E, ) h , > - E . I ® >] 

^ o I m**! rn*“k n 

( 1 . 42 ) 

Defining R = P/(E - H ), 

O Q 

the energy expression upto first four orders is 

o o 

= < V > 


ttVt V' 


o O ■ 


A E 


( 3 ) 


= <V R V> 

A A 

= ‘*0 rrVr “ TrV)V i v- 

a o 0 o 

A 

<<5 jv i« ><<f IV rr -- g -T2 V M > 

o' ' □ o ( E -R ) ^ o 

a o 

= <VRVRV> - <V><V R^V > 


( 1 , 43 ) 


( 1 . 44 ) 


= <VRVRV> - <VR <y > RV> 


( 1 . 45 ) 
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IV R(\7-E^ )R(\/-E^ )RV I > -E 2 < h > 
= <\/ R V R VR \J > •- <y R<V> R V RV>- 
<0 R V R<V'R V> + <V R <V>R<V>R V>- 


A. 4k A 


▲ ▲ 4i A * A 4- 

<V R<V R V>R V> 


( 1.46) 


Irrespective of whether the model space is complete or 

finite, the energy series converges as long as the perturbs- 

1 

tion theory itself converges. When the unperturbed 

A 

Hamiltonian is chosen to be the ground-state HF operator, 
the Mj6ller-Plesset [33] version of the RSPT is obtained. 
This form has been applied extensively by Pople [ 55] to 
carry out correlation energy calculations upto second and 
third order. 


The Ralei gh-5chrodinger perturbation theoretic method 
has two inherent drawbacks; firstly, it shares in common 

with Cl, the limitation on the size of the space defined by 

* 

P. Usually, an intuitively appealing approximation like 
limiting 1 x > "to single and double excitations of the 
reference wavefunction , j > > is made. For perturbation 

expansion upto three orders this approximation is exact, 
since only single and double excitations give non-vanishing 
matrix elements. The first-order correction to the wave- 
function is then : 



-I 

a<B r< s 


Ivl<5„p > 


^r"'’ ^s~ ^3 


(1,47) 
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I I 

a<3 3c < s 



~rs(1) ^ ^0 » 

( % - c„ - e^ ) 

2 ? 5 

where > are obtained by substituting the occupied 


where C 


(1.48) 

(1.49) 


orbitals a and (? in the HF ground-state determinant 
by the virtual orbitals r and s; e‘ s are the eigenvalues 
of the one-electron Fock Hamxltonian operators. Even 
with the restriction of the space to single and double 
excitations, it is not practically feasible to go beyond a 
few orders. The astronomical growth in the number of 
terms as the order of perturbation increases is the second 
limitation of RSPT. 


1.4.6 Linked-Clu ster Theorem and Manv-bodv Perturbation 
Theory 

The application of field-theoretic techniques to the 
above problem heralded a new epoch in the solution of many- 
body problems. The first of the genre is the famous 
linked-cluster theorem of Goldstone [ 56] . This theorem 
uses the time -dependent perturbation theory in the interaction 
representation, the language of second quantisation, and 
Feynman's graphology (vide Appendix A), to prove rigorously 
the cancellation of unlinked terms to all orders in RSPT. 

There are two very important features associated with the 
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linked-cluster theorem. Firstly, it can be formulated in 
an entirely time -independent framework [ 57,58,59 ] . 
Secondly, its applicability extends well beyond perturbation 
theoretic methods [ 6o] . Numerous methods, built on this 
theorem, eventhough non-v/ariational in nature, are devoid 
of the si ze-extcn sivi ty problems dealt with earlier. With 
the application of the linked-cluster theorem eqn. (1,39) 
becomes 

00 

AE = I V]*"! (I.5D) 

n=::Q ^ 

where the subscript indicates that the summation is 

limited only to linked terms. In other words, the proje- 
ction operator P in eqn. (1.39) which excludes the 5^ 
component of any function it operates on, has the same 
effect as the restriction to summation over linked diagrams. 
Also, the myriad linked terms thab survive are elegantly 
handled, with the graphical representation of Feynman [61], 
The Feynman diagrams, in addition to handling the book- 
keeping also have provided insight for selective summation 
in eqn. (1.50) [62,63], 

The power of the linked-cluster theorem is highlighted 
by the manner in which it gets rid of the si ze-ex tensi vity 
problem mentioned earlier with reference to Cl. RSPT shares 
the si zB-extensivity problem with Cl. On the one hand, it 
is necessary to limit the basis to double-excited functions 
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only, while on the other this approximation introduces 
unphysical, spurious energy corrections. This aspect 
becomes transparent in tho fourth-order of R5PT and has 
been illustrated in Appendix A. The application of the 
linked-clus ter theorem, however, immediately solves the 
sizc-Gxtensivity problem. There ore no unlinked terms 
that figure at all in the summation and all the linked 
terms have correct behaviour for N dependence. The eli- 
mination of the SI ze-ext en SI VI ty problem portrays only one 
facet of the linkod-cluster thoorern. The theorem is much 
more versatile as will become apparent from the formula- 
tion of non-perturbative theories and this has been 
dealt with in Chaprer II of this thesis. The above des- 
cribed many-body approach to RSPT, that goes under the name 
of many-body perturbation theory (MBPT) has been applied 
widely to atomic and molecular correlation problems [64], 

1.4.7 Other Field Theoretic Methods 

Notwithstanding the theoretical elegance, MBPT suffers 
from some major drawbacks. Firstly, the perturbation series 
for extended systems involving coulombic forces is itself 
divergent. And, despite the elimination of all the unlinked 
terms there still remain the numerous linked diagrams to be 
taken into account. However, in practice, this is not 
feasible in the infinite summation. Although, one of the 
powers of the linked-cluster theorem is to provide topological 
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arguments for selective summation of diagrams to infinite 
order [62], physical insight into the meaning of the per- 
turbation terms is lost. Also, if only partial summation 
IS done, infinite number of diagrams vjhich may contribute 
significantly to the perturbation correction are ignored. 

This casts doubt on the choice of diagrams invoked in the 
partial summation. Lastly, the MBPT, can boast neither of 
the conceptual simplicity nor the physical insight into 
correlations of electrons provided by Cl. 

A method that retains a3.1 the advantages of the linked- 
cluster theorem, but at the same time avoids the short- 
comings of the MBPT was developed by Coester and Kummel 
[65,66] and extended to many-electron systems by Ci'zek [35]. 
This method known as the CC or the Bxp(5) method [67], is 
based on the cluster expansion of the exact wavefunction 
using an exponential ansatz. The linked diagram property 
of the exponential operator and the effective interaction 
IS exploited and sets of explicit CC equations for the 
components of the cluster operator obtained [35], The CC 
method is the most powerful off-shoot of the linked-clust er 
theorem. Though, non-perturbative in origin, it displays 
the theoretical rigor of true many-body treatment and bears 
bhe conceptual simplicity of a physically meaningful 
theory. The details of the method and the diagrammatics 
are presented in Chapter II of this thesis. The other field- 
theoretic method, the Green's Functions method [36], has 
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only been sparsely used in quantum chemistry, due to its 
complicated nature. 

1.5 BEYOND IPA WITH SPECIFIC REFERENCE TO LOCAL POTENTIAL 
METHODS 

The methods described for going beyond the IPA, are not 
restrictive as far as the choice of the initial approximation 
IS concerned. However, the HF approximation is the preferred 
starting point for most calculations involving these methods. 
This IS so because, the HF operator is the exact Hamiltonian 
for the single determinant trial function. Besides, the 
terms exchange and correlation are traditionally defined with 
reference to the HF Hamiltonian [ 26] , which fact is used to 

advantage in decomposition of the Hamiltonian into a model 
Hamiltonian and a well-defined perturbation. Above all, 
the HF solutions rigorously satisfy Brillouin’s theorem, 

<0^ I H > = 0 (1.51 ) 

This considerably simplifies Cl and perturbation calculations 
with the HF basis. 

The local potential methods, on the other hand, do 
not share the above merits of the HF method; however, they 
score over the HF method in the computational domain and are 
better suited for large system calculations [60]. Also, these 
methods are conceptually very attractive since they are based 
on a physically well-defined, measurable quantity, the one- 
electron density. However, unlike in the HF method a 
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clear-cut demarcation of the exchange and correlation 
effects cannot be made in these methods; nor is it clear 
how exactly the local potentials simulate the non-local 
electron-electron interaction potential. Attempts hav^e 
been made to highlight the nature of these potentials by 
introducing local and non-local corrections to these poten- 
tials [ 69] . Main amongst these attempts are (a) the 
introduction of ompirical correlation potentials based on 
phenomenological or conceptual arguments, (b) the incor- 
poiation of gradient corrections to kinetic and exchange 
energies and (c) the modification of the local potential 
to reproduce the correct shape of the Coulomb and Fermi 
holes. In the class of methods (a), the EC potential 

in eqn. (1.10) assumes the form 

^XC " ^ o(r)dr + / e|^P(r)dr (1.52) 

where e is a functional of density and defines the 
correlation energy per electron. e has been variously 
approximated by correlation energy of a f ree-electron 
gas [ 70 ], correlation energy of a low density [71] and 
high density [72] electron gas obtained using perturbation 
theory, correlation energy of a homogeneous electron gas 
obtained through a dielectric formulation [73,74], corre- 
lation energy of a spin-polarised electron liquid [75], 
etc. All these approximations are valid only for slowly 
varying uniform density distributions [l3] and are incapable 
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of taking into account density fluctuations due to atomic 
shell structure [iS]. These potentials also have the in- 
herent error of inclusion of self-correlation. This arises 
from the treatment of the electron as a continuous electron 
density and the subsequent calculation of the correlation 
between different parts of this density [76]. This self- 
correlation leads to a non-zero correlation energy for the 
hydrogen atomi [75], Method (b) involves the incorporation 
of gradient corrections to kinetic energy [77] and to exchange 
energy [72] and leads to fairly accurate valence charge den- 
sities, but divergences are encountered in the inner regions. 
Also, this has the theoretical flaw of not obeying the sum 
rule, which states that the EC hole contains one electron [69], 
The last mentioned of the methods, (c), attempts to correeb for 
the shape of the Fermi and the Coulomb holes [75,78,79,80] . 
Although, these methods are conceptually very attractive, they 
are difficult to implement. 

All the above mentioned attempts that correct the local 
potential approximation are still within the domain of the IPA. 
Zare [81] has performed a Cl calculation on magnesium starting 
with a basis obtained from the X® method and has obtained 
excellent results. However, few attempts have been made to 
go beyond the local potential based IPA using the other 
methods described in Section 1.4. In particular, it may be 
pointed out, that none of the more powerful approaches such 
as the linked-cluster methods, have been employed in this 


context . 
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1,6 SCOPE OF THE IHESTS 

Although the above mentioned attempts tend to refine 
the local potontial method to higher degrees of sophisti- 
cation, they fail to provide insight into the nature and 
efficacy of the local EC potential. In view of the 
growing importance of these methods for large-scale compu- 
tations, it is not only desirable, but necessary to examine 
these models with a view to understanding the nature of and 
the basis for the approximations used. This thesis, an 
outcome of this motivation, uses the powerful CC many-body 
theory to study the many-body EC corrections to a member 
of the local potential family, namely, the Xet method of 
Slater [ 8 ] . 

Cizek and Paldus [35,82] have employed the CC approach 
to study the many-body correlation effects in atoms; the HF 
approximation served as the starting point for their cal- 
culations. However, to the best of our knowledge, there 
have been no attempts made in literature to use the CC 
method for analysing atomic EC effects, starting from a 
local potential-based reference function. Such a study 
would, in addition to incorporating many-body correlation 
effects, provide an insight into the nature of the IPA 
chosen for the reference function. 

In the present thesis, we examine in detail the role 
of the Slater exchange potential and the choice of the 
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exchange parameter a y in the treatment of many— body 
exchange correlation effects. The use of gradient corre- 
ctions to the X c* potential is also examined in this light. 

The present study entails, in addition to the use of CC 
approach, the use of several hierarchy of approximations 
derived from it, namely, variants of coupled— electron pair 
theories [31,46], Cl [37], Sinanoglu's decoupled-pair many- 
electron theory [4l], IMesbet-Bethe-Goldstone method [ 4D] 
and PS second— order expansion [55] > Our results are compared 
with those of Cizek and Paldus ; such a comparison provides 
a framework for assessing the relative merits of the HF and 
Xct potentials in the description of EC effects. The use 
of different local potentials to generate the reference fun- 
ctions for the CC calculations, has permitted a detailed 
study of the role of these potentials vis-a-vis, the EC 
effects. The diagrammatic language facilitates the analysis 
of the importance of the terms that contribute to EC 
corrections and thereby provides an insight into the nature 
of the local potential used. An assessment of the efficacy 
of the potentials, in representing EC effects in various 
regions in the atom, has been made. It is hoped that a 
study such as the present one will provide a rational basis 
for the choice of local potentials that describe the EC 
effects in atoms more accurately. 

Chapter II of the thesis discusses the CC method, its 
basis, diagrammatics and the cluster equations, in the 
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context of evaluation of corrections to EC energy in atoms. 
Chapter III presents a model calculation on beryllium atom 
to illustrate the application of the CC methodology and the 
approxima bions derived from it. Chapter IV presents the 
results of calculations of corrections to EC energies in 
four-electron atomic systems obtained by the CC approach. 

The X a method end its refinements have been examined 
through these calculations. Chopter V reports the expecta- 
tion values of some operators over the wavef unctions obtained 
through the CC expansion for the systems discussed in the 
previous chapter. The results are compared with the reference 
state expectation values. The thesis concludes with a few 
salient remarks on the application of CC methodology for the 
study of EC ofFects in atoms. 

Appendix A presents a brief discussion on the language 
of second quantisation, Feynman-Goldstone graphology as 
applied to MBPT, fourth-order RSPT and cancellation of the 
unlinked diagrams, the hole-particle formalism, Wick’s 
theorem and contractions and illustration of the elementary 
Feynman-Goldstone diagrams used in this thesis, A listing 
of some of the computer programs used in the present study 
is given in Appendix B. 
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THE COUPLED-CLUSTER FORMALISM 


11. 1 Introduction 

11. 2 The Nature of T 

11.3 The Schrodinger Equation 

11.4 Diagrammatics and CC Equations 

Relationship of the CC Method to othe 
Many-Electron Theories 


II. 5 



IK1 INTRODUCTION 


The CC methodology^ which is the most significant off- 
shoot of the linke d*-clu ster theorem is used in this thesis 
for the analysis of EC effects in dtoms. This method is 
based on the cluster expansion of the exact many-Fermion 
wavefunction in an exponential Form [ 35 , 65 , 66] . This 
exponential ansatz, which is a facsimile of the linkcd- 
cluster expansion, was first inferred from the structure of 
many-body perturbation theory [ 63] , A diagrammatic Raleigh- 
Schrodinger expansion of the exact wavefunction elegantly 
proves this point. The wavefunction expansion can be equi- 
valently expressed by diagrams that contain one or more open- 
linked pieces (cf. Appendix A for diagram terminology) . Some 
of these diagrams have topologically identical linked pieces, 
differing only in their relative arrangement. It turns out 
that diagrams with n such topologically equivalent linked 
pieces are counted nl times thanks to the various orderings 
of the linked pieces and independent orbital summations. It 
is therefore necessary to multiply the contribution from the 
diagram by (ni) . This is illustrated for the case of a 
diagram containing two open-linked parts in Fig. 11,1, 

For the general case of a diagram with n linked parts 
of topological structure r, , n linked parts of topological 
structure X 2 t etc., the multiplication factor for the contri- 
bution is [ n (n^l)]~^. The net outcome of these factors is 
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that the total wavefunction can be expressed in the "form [57] 


'i' = e 


(II. 1 ) 


where, T represents the sum of all open diagrams consisting 

A 

of just a single linked piece. As implied above, T can be 
written in the form. 


T = 


(II. 2) 


where, is the fully linked operator for exciting n- 
particles out of the Fermi sea I represents the single 
determinantal ground state. This result can also bo verified 
independent of perturbation theory [84]. The transformation 

A 

exp( T) is non-unitary, since the expansion satisfies the 
intermediate normalisation condition 


<'y 0 > = 1 

o 


( II.3) 


II. 2 THE NATURE OF T 


The nature of the operator t is best understood 

A 

through the second quantisation formalism. The * s 
defined such that, 


r. *+ 

""1 1 


„ ^<^2 *t *t *t 

I ^ “•'“2 ^2 

a. 112 2 

12 12 


Ji r^r^... r, ^-j- 

t 4 I ^ ^ ^ ^ 

k. “102*** “k „ „ a ^ 

r^r^... k 
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where, the b ’s are the particle-hole creation operators 
(vide Appendix A) defined on the complete set of states, 

^ 2 *** represent hole states while r ^ jX ^*** T'spresent 
particle states. The function of the operators in Tj^ is 
to excite k particles out of the Fermi sea. In the hole- 
particle formalism, this is equivalent to saying that k 
particle-hole pairs are created. The are the corres- 

ponding amplitudes of excitation. The following two pro- 
* 

pertios of T oporaiors become apparent from their defini- 
tion. Firstly, T's are non-local operators and have non- 
vanishing amplitudes between hole and particle states. And, 

A 

the operators commute with each other, since they can 
only create and never annihilate porticle-hole pairs. The 

A 

operators can be represented graphically using Feynman- 
Goldstone diagrams [56], This is shown in diagrams A-F of 
Fig. II. 2, A solid lino joining the n particle-holc linos 
IS drawn to show that the n particle-hole pairs are created 
simultaneously. The arrows indicate that each upgoing 
line represents a Fermion creation operator and a down- 
going line represents a Fermion annihilation operator. 

The diagrams D,E and F show the disconnected diagrams 

A 2 ^ ^ A 2 

Ti , "^2 "^2 I'espectively . 

The incorporation of the many-body effects through 
the exponential operator becomes transparent on comparison 
of the cluster expansion with Cl expansion [05]. Writing 
the Cl expansion of the exact wavefunction '5' , in the form 
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< 


FigJl.2 A diagrammatic representation of t operators. 
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n 

I '?> = Uo> + ^ l^n> (II. 5) 

S=1 ® ° 

where, is an s-fold excitation operator such that 

A 

I ^ yields the sum of all the k s-fold excited con- 
figurations, 


Es l»n > = , J: “k |5''> 

kes 

It is seen an comparing eqns. (II. 6) and (II. 1) that, the 

A A 

operator can be decomposed into a connected part and 

^ A 

a disconnected part The can be written as product 

.f f^-S. 

^ (II. 7) 


where. 


U. = 


>^2 = 21 


U. = 


4. 


i_ f 

31 '1 


U, = 


k "^2 + k ^1 k ^2 


( II.B) 


etc. The exponential ansatz then is seen to be entirely 
equivalent to a Cl expansion [86], 

II. 3 THE SCHRODINGcR EQUATION 

The Schrodinger equation, eqn. (1.1), on substituting 
for T from eqn. (II.1), assumes the form, 

A A 

H b"^ I 0 >= E e"*" I 5)„> 

» o ’ o 


(II. 9) 
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Left-multiplying eqn. (II. 9) successively by each of the 
s Q yields the set of eqns. (11.10), which are 
the constraints on T. 


< |H|e^ t > = 0 


|H|e^ > = 0 


( II .10) 




rs 


T 


a6 




Cf > 

o 


= 0 


These constraxnts (determine T completely. The set of 
constraints yield sets of coupled non-linear equations 
linking the various n-body amplitudes t^ to each other. By 

jk. 

suitably truncating the series , computationally feasible 
approximations are obtained. The energy E is obtained by 
left-multiplying eqn. (II. 9) by and substituting the 

appropriate ' s determined from the constraint equations, 

A pedestrian derivation of the constraint equations is 
tedious and fails to reveal the inherent properties of the 
CC formalism [65,36], Use of field-theoretic methods 
greatly simplifies the problem and provides insight into 
the nature of many-electron interactions [35], For the 
derivation of the CC equations, we then employ the tools 
of field theory, viz., the hole-particle formalism, Wick’s 
algebra and Feynman-Goldstone diagramma ties (vide 
Appendix A) . 
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The Hamiltonian can be expressed in the second 
quantised form as (vide Appendix A) 


= I 

ik 


< i z 


k 


> a , 


y2 l < ij 

ijkl 


* t 1 •, *+ At 

V kl > a . a . 

1 J 


=1 =k 
( 11 . 11 ) 


or alternately, in the normal product form as 


H = < I H I > + I < i I f I k > N [ ] 

ik 

+ y2 I < ij jv|kl> N [at a, a. 1 (11.12) 

ijkl r 3 1 k 

where , 

<ilfjk>= <ilzlk> + < i j gj k > (II. 13a) 

<i Ig |k > = j;<i C |vj k? > ( II .1 3b) 

K A 

z refers to the sum of kinetic energy and nucleus-electron 
attraction operators and v the electron-electron repulsion 
operator. The subscript A on the second term in Gqn.(II.13b) 
stands for an trsymmetri sation of the matrix element. When 
i,k and C refer to orbitals, 


I <iC lvjkC>. 2<i5 |0lkC>- <i?|v|ck> (II. 13c) 

S a ^ 

The advantage in expressing the opBrat(3rs in the normal 
product form lies in the obviation of the need to perform 
contractions of the operators within the normal product 
(vide Appendix A). 
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Eqn* (11,9) can now be recast in the second quanti- 
sation formalisin and expressed diagranmatically , To obtain 
eqn, (II, 9) in a compact form, the operator exp(T) is re- 
placed by the formal laylor series expansion, 


* 

T 

e 


1 + T + 


II 

2i 


(11.14) 


then , 

* * *2 

H{1+f+ — + ...) ^ = E(1+f+ + ••♦) $ (11,15) 


The Hamiltonian has to be contracted with the T ^ s of the 

n 

Taylor series expansion; in other words, the creation- 

A. 

annihilation operators in H are contracted with creation- 


annihilation operators in , After these contractions are 
performed, the Schrodinger equation, eqn. (II. 9) assumes 


the form, 


{fl e"^} 


y o o 


( 11 . 16 ) 


The subscript implies that only topologically linked 
diagrams need be considered. This result, the non- 
perturbative analogue of the linked-clus ter theorem, is 
easily proved [ 57] . 


On contraction with exp(T), the diagonal part of the 

A 

one-body operator H can link itself diagrammatically to 

no mors than two of the In the most general case, this 

can be illustrated by contracting the creation operator of 

"fc h 

the one-body operator with one of the k particle-hole pair 
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and contracting the annihilation operator with one 

th . “*■ 

of the 1 particle-hole pair of T^. 


■f 


1 ^ r / .*t -r t 


(SI S,) ((e; b; )(£+ s; ) . . ) 

112 “2 k k 


wki* jk ^ 

• « ( t} ^ rv 

r* (X 

1 1 


). 


* t *t ^ 
) 

n n 


e {H} 


e{T> 


The diagram obtained from this contracted term has none, one 
or two open~paths; each open-path denotes a particle-hole 
creation pair. The operator that corresponds to this diagram 

A 

commutes with all the T^'s. Owing to the commutation property 

A 

of the T 's within themselves and the T 's with operators 
n n 

obtained on contraction, there factors out on the left, a 

A 

factor exp(T). Similarly the two-body operator of the 
Hamiltonian can be contracted with one, two, three or four 

A 

of the The diagram obtained on contraction can again 

have none, one or two open-paths? and the operator corres- 
ponding to this also has particle-hole creation operator 

A 

pairs and hence commutes with all the T^^'s. The result is 
again the factoring out of exp(T) on the left. Eqn. (II. 9) 
can, therefore, be written as, 


e”^ t H e"^ ^ = e"*” E <5 


(II. 17) 

A 

-T 


Left— multiplying both sides of eqn. (11.17) with e yields 
eqn, (11.16), 

A diagrammatic version of the proof [B7], for the case 
of f = is presented in Fig. II. 3, for the sake of clarity. 
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Tha quantity 
diagram can 


below the dotted line yields the factor and the combined 
be written as the product of eT and a linked port. 


Fig.n,3 A diagrammatic 


pjcturlsation of the linkcd-cluster result. 
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This result can also be obtained, without the aid of 
diagrammatics, by rewriting eqn. (II. 9) in the form, 

^ A 

e“^ H e''' f = E 0 (II.1B) 

o o 

and using the Baker-HaUsdorff- Campbell formula [88] . 

H = H + [H,f] + ~ [ [ H , T] , T] + ... (11.19) 

The commutator series terminates after the fifth term, since 
H is atmost a two-body operator and expansion of the result- 

A 

ing commutator series yields eqn* (11*17) [ 67 ] • When H 

IS expressed in the normal product form, eqn* (II. I 6 ) becomes 
an eigenvalue equation for A E, where 


AE=E-<0 H $> 

o O 


( 11 . 20 ) 


The eigenvalue equation is given by, 

A 

|h1 0 >?> + {H e"'"} ) $ = Eo 


o 

i * e • , 


( H e^)',4>= A E * 


The prime on LHS expression indicates that. 


( 11 . 21 ) 


{ H e"’" f = { ( H- <C ! H I ?■ > ) B ' } 


( 11 . 22 ) 


Eqn- (11.21) can be rewritten as, 

Jk 

< U H e'''}^! = AE (11.23) 

SincB, the RHS in eqn* (II.23) is a scalar, only closed— 

^ T 

linked diagrams of the expression {H e need be taken into 
account for AE evaluation; the open-linked diagrams correspond 
to particle-hole pair creation operators. The diagrams 
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corresponding to eqn. (11,23) are presented in Frg . 11,4. 

To obtain explicit expressions for AE, the normal 
product form of the Hamiltonian may be rewritten as, 

( I H 1 0^>) = < ijh^ jk> N aj^] + J <1] g-v^l k>N[a'^ 

ik ' ik 

+ 1/2 I <ij [v| kl>iM[ a| a't a, a.] (11.24) 

ijkl 1 j 1 k-i 

A 

whore is any one-partrcle Dperator diagonal in the given 
□ibital basis 

[i > - ( z + v^) li> 

= £^li> (11.25) 

The potential v is any single particle operator and in the 
present study it is defined through the X^i local potential 
(vide Chapter I). The one-body part of the operator in the 
LH S □ f eqn. (II. 21), f, is given by the first two terms 
on the RH5 of eqn. (11.24), while the two-body part of 
the operator is given by the last term. Using this defini- 

A 

tion for H and the standard rules for constructing Feynman- 
Goldstone graphs, an explicit equation for A E is obtained. 

A 

If h^ happens to be the Fock operator, AE represents the 
correlation energy. In our case, AE actually represents the 
corrections to EC energy (vide Section 1.3). The expression 
for A E from Fig. II. 4 is then, 

AE = I <alg-0 jr> tj + I <aBlvlrs> . tj t| 
ar aers 

+ I xaep! > A (11.26) 

aBrs 
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Diagram A of Fig. II. 4 corresponds to the first term in the 

RHS of eqn. (11.26), B and C correspond to the second term 

and D and E to the last term. A knowledge of the amplitudes, 

which are the cluster coefficients, t^ and t^p,is necessary 

and sufficient to obtain AE. The amplitudes t^ and t^^ 

a a 0 

are obtained from the cluster equations, that result from 
projection or eqn. (11.21) with and as functions 


l{H e^}^ I 0 


Cr e }’ I * >.= 0 


(11.27) 


II. 4 DI AGRAMMATICS AND CC EQUATIONS 


The closed-linked diagrams of the LH5 of eqn. (11.21) 
provide expressions for the correlation energy; on the 
other hand, their open counterparts with one , two , ... etc. 

open-paths represent one or more particle-hole pair creation 
opearators and therefore have no corresponding terms on the 
RHS of the equation. Hence, they must independently be 
equal to zero. Eqns. (11.27), therefore, may be expressed 
die grammatically and the resulting diagram topologies pro- 
vide meaningful insight into the nature of the diagrams. In 
bhe ensuing description Feynman-Goldstone diagrammatic s , 
which represent non-antisymmetrized matrix elements, are 
employed. The rules for obtaining these linked diagrams are 
simple (vide Appendix A). The "H” diagrams are linked with 

A 

the "T” diagrams of each term in exp(T) in accordance with the 



contraction rules for the creation-annihilation and particle- 
hole pair creation operators. The resulting diagrams are 

V 

referred to as "R" diagrams by Cizek [35,89] . 

The R diagrams have connected lines, free lines, 
closed loops and H and T vertices. The rules for obtaining 
matrix element expressions with proper weight factors and 
signs are the same as that employed in Goldstone perturba- 
tion theory (vide Appendix A). However, since there is no 
projection operator implied in those diagrams, energy 
denominators don’t arise. By virtue of contraction rules 
for creation-annihilation operators, the following assignment 
can be made for the lines in the R diagrams; connecting lines 
that begin from a H vertex and enter a T vertex carry hole 
labels and those that begin from a T vertex and enter a H 
vertex carry particle labels. In the case of free lines, 
those entering the vertices of T or H diagrams carry hole 
labels, while those leaving the vertices carry particle 
labels. The amplitudes t^ , t^f, are also referred to, 
as the t matrix elements <ct|t|r> , <oPjt|rs> to preserve 
a ^old stone-li ke diagram correspondence. 

Diagrams with one open-path «. 

All the topologically distinct one open— path linked 
diagrams obtained by linking H and T diagrams are given in 
Figs. II. 5 and II. 6. The diagrams, each representing a 
particle-hole creation operator pair are obtained by,' tH'^'lljK 
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Fig. 11*5 (continyed) 
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R 

( 2 ) ( 2 ) 


K 


Fig.n.5 On« open-path diagrams yielding 
terms linear in "t". 
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“following linking bBiwcen H and T diagrams i 

a) f diagrams with none of the T diagrams (diagram A oi 
Fig. II. 5) 

A 

b) f With a diagram (diagrams B and C of Fig. II. 5) 

c) V with a diagram (diagrams D and E of Fig. II. 5) 

d) f with a diagram (diagrams F and G of Fig. II. 5) 

e) V with B diagram (diagrams H-K of Fig. 11.5) 

A 

f) f with two di agrams ( di agram A of Fig. II. 6) 

g) V with two diagrams (diagrams B-E of Fig. II. 6) 

h) V with a and a T 2 diagram (diagrams F-M of Fig. II. 6). 

The contributions to the one open-path linked diagrams from 
are represented by diagrams A-E of Fig. II. 5, while those 
for T 2 are represented by diagrams F-K of Fig. II. 5. All 
those terms are linear in the coefficients and give the 
linear contributions to the cluster equation. Diagrams 
A-E of Fig. II. 6 arise from linking (two diagrams) 

with f and v diagrams, while diagrams F-M of Fig. II. 6 arise 
from linking of T^ T 2 (a and a T 2 diagram) with v dia- 
grams. These give the non-linear contributions to the 
equations containing one open-path linked diagrams. No 
other 'one open-path' diagrams result owing to the linked- 
eJuster theorem represented by eqn. (11.21). The equations 
corresponding to these diagrams are given in Eqn. (11.28). 
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'^(2)(2,1)= ^ (4 <e 6 1 v|wu> - 2 <eSl v! uw > - 

c uw 

I (2 < 6 e j v j uw> - <5 ej v| wu > ) - 

e uw 

1 [2 <(Sejvlwu> - <(Se|vluw> + 

' ' ' ' 0 cxe 

e uw 

^ (2<5elvluw>-<6ej<?jwu>)tgt^g . 
e uw 

The following notation is used to label the terms from the 
Various R diagrams. Each class is represented by 
superscript indicating the number of open-paths, the subscripts 
indicating the elementary diagrams involved in the linking; 
i is equal to 1 or 2 depending on whether an f or a v diagram 
is involved in the linking; k indicates the operator invol- 

A. 

ved in the linking process. If more than one T operator is 
involved in the linking it is indicated in the second subscript 
as (k,l), etc. 

The one-particle operator f (eqn. II. 13a) is diagonal if 
HF basis is used and the matrix elements automatically vanish 
if f connects different single-particle states. For any other 
basis ^ is not diagonal and the matrix elements connecting 
different states are given by 

<i|flj> = <i|g-Vg|j> for i j 

<ijf |i> = e. + <i ! g~Vg|i> 


(11.29) 
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The single-particle states i,j in this case are the Bigen-* 
functions of eqn. (II,25). 

Di_a 9 rams with two ooen-paths : 

All the topologically distinct two open-^path linked 
diagrams obtained by linking H and T diagrams are given in 
Figs# 11,7, II#8 and II,9# Each diagram represents a pair 
of particle-hole creation operators* These diagrams are 
obtained by the following linking between H and T diagrams* 

a) V with none of the T diagrams (diagram A of Fig* II#7) 

b) V with a T^ diagram (diagrams B and C of Fig# II. 7) 

A 

c) f with a T 2 diagram (diagrams D and E of Fig. II. 7) 

d) V with a diagram (diagrams F-K of Fig. 11,7) 

b) V with a ^ (two ' s) (diagrams A-C of Fig. II. 8) 

f) V with a (two s) (diagrams D-*N of Fig. II. B) 

g) V with a and a T 2 (diagrams A-N of Fig, 11,9) 

Diagrams B and C of Fig. 11,7 arise from diagram linkings 

while diagrams D-K of Fig. 11,7 arise from T 2 diagram linkings 

All terms represented by these diagrams are linear in and 

2 

T 2 . Diagrams A-C of Fig. 11,8 originate from linking of 

2 

diagrams, D-N of Fig. II. 8 from linking of T 2 diagrams and 
A to N of Fig. 11,9 from linking of T 2 diagrams. These 

give the non-linear contributions to the cluster equation 
with two open-paths. The corresponding equations are given 
in eqn. (11,30). The topologically identical diagrams with 
the index permutations a >> U , r<— $s are not drawn 
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Fig. II. 8 Two open-path diagrams (involving Tj) 
yielding terms non-linear in 't’ 
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explicitly; however, the terms arising from these diagrams 
are written explicitly in eqn. (11.30), 

The two open-path cluster equations are given by : 


p^2) . r(2) , A2) .(2) 

'^( 2 )( 0 ) ( 2 ){ 1 ) ^ ^{ 1 )( 2 ) ^ ( 2 )( 2 ) 


+ R 


(2) 

( 2) (1 , 1 ) 


+ 


r(2) 

( 2 )( 2 , 2 ) 


+ R 


( 2 ) 

( 2 )( 1 , 2 ) 


0 


(11.30) 


where, 
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I [ <36|''|xv.> tjt“^ + 

I [<rei;|ux> + 

e ux 

I [ <5rilv| wp > + 

I [< en| V |uB > + 

en Li 

I [<se|vjxu> + 

e ux 


<r6|vlxw> - 

<se|Oj ux> J + 

n! vj wa > ] + 

f -^1 . u . sr 

<en|vl ua > ] - 

< re Ivj XU > 1 

'I P ea J 


The coupled non-linear equations, eqns. (II. 28) and 
(11.30) yield as solutions the cluster coefficients, which 
can then be substituted in eqn. (11.26) to obtain A E, 

1 1. 5 RELATIOtMShIP OF THE CC METHOD TO OTHER MAM Y-ELECTRO N 
THEORIES 

The relationship of the CC coefficients to the Cl 
coefficients was outlined in Sec. 11,2. The CC method is 
an extremely elegant way of incorporating the many-electron 
correlations. In contrast to Cl, the method is rigorously 
size-extensive and avoids the messy problem of dealing with 
Slater determinants. However, the method is non-variational 
and can lead to energies lower than the exact energies [ 90] . 
Variational counterparts of the CC method, although, are very 
complicated, still involve less computational effort than 
full Cl methods [31]. 
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The relationship of the CC method to the perturbative 
methods is even more obvious; after all, as shown in 5ec.II.1» 
the basic form of the wavefunction in the CC method is a 
direct result from R5PT, The perturbative flavor of the 
method becomes apparent from the comparison of the Goldstone 
diagrams in the two versions of the theory. The cluster 
equations, eqns. (11.20) and (11.30) can be written in such 
a manner that a diagrammatic expansion is obtained for each t. 
When this is substituted into the energy diagrams (vide 
Fig. II. 4) taking into account the appropriate label matching 
of the lines, a Goldstone-like expansion is obtained for the 
energy. Iteration of the equations leads to generation of 
Goldstone diagrams to higher orders. For instance, inclusion 
of the diagram K of Fig. II. 7 in cluster equations is equi- 
valent to performing an MBPT partial summation to all orders 
of the corresponding ladder diagram [62], Similarly, the 
diagram F of Fig. II. 7 generates the corresponding ring dia- 
grams of perturbation theory to all orders. This is illu- 
strated in Fig. 11.10. Alternatively, this can also be 
pictured from perturbation theory. Snipping off the top-most 
vertices of the Goldstone perturbation energy diagrams results 
in the open-path diagrams that have topological structures 
similar to the diagrams in Figs. II. 6 - II.9[60]. It is 
obvious from this that the iE obtained by solution of the 
coupled equations, eqns. (11.26), (11.28) and (11.30), 
invokes perturbation theory diagrams that far exceed the 
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number that can be included by MBPT calculations. Also, by 
choosing only certain diagrams in eqns, (11.28) and (11.30) 
it is possible to obtain selective perturbation summations 
to infinite orders. For instance, retaining only the diagrams 
A and K of Fig. 11,7 leads to a result equivalent to selective 
summation of particle— particle ladder diagrams to all orders 
in perturbation theory. Likewise, consideration of all 
ladder diagrams in this theory yields equations similar to 
the Bethe-Gcl dstone equations of nuclear many-body theory [ 91] • 
On the other hand, retention of only the diagrams D of 
Fig. II. 5, A and F of Fig. II. 7 and D of Fig. II. 8 in the 
cluster equations yields solutions equivalent to those obtained 
by inclusion of ring diagrams to all orders in perturbation 
theory . 

As shown by Paldus et al. [ 89] , the T 2 approximation 
to the CC method is. a generalisation of Sinanoglu’s decoupled- 
pair many-electron theory (DPMET). When the occupied orbitals 
in the CC equations for the'T 2 only' approximation are restri- 
cted to a pair of spin-orbitals , and the inter-pair neglected, 
each surviving term assumes the form of a product of a pair 
correlation energy and the corresponding t 2 amplitude. The 
CC equations, therefore, reduce to a number of independent 
equations one for each pair of occupied orbitals. These 
equations are exactly those of Sinanoglu’s DPMET theory. 

The CC method, it is apparent from the above discussion, 
is superior to both Cl and MBPT methods at least for many- 
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electron systems. The closed-shell version described here 
has been generalised with the aid of graphical methods of 
angular momentum to open-shell systems by a number of 
workers [ 59,60,92,93 ]. This open-shell version of the 
theory has also been shown to be very superior to all the 
existing open— shell many-electron theories. An excellent 
review of the CC method and its applications has been 
presented by Kummel et al [67] » The applications have 
included the analysis of correlation effects in atoms and 
molecules [35,89,94], analysis of exchange effects via 
inclusions of ring diagrams in electron gas [95], and the 
analysis of long-range forces in nuclear matter 167]. 


I n the 
application 
corrections 


next chapter, we present as an example 
of the CC methodology, the calculation 
in beryllium atom. 


of 

of 


the 

EC 
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III.1 GENERAL 

This chapter presents the procedural details of 
the computations performed in this thesis and a model 
calculation to illustrate the CC methodology. The X ot 
approximation of Slater[8] is our starting point. The 
self-consistent field orbitals obtained by this method are 
used to generate the matrix elements and to calculate the 
total electronic energy of the atom. The CC equations 
for various hierarchy of approximations are set-up and 
solved using iterative algorithms. For all the cases a 
partial wave analysis of the EC energy corrections is 
presented. 

The beryllium atom has been chosen as the model 
system for this study in view of the following : beryllium 
is the simplest example of a true many-ele ctron system. 
Almost every tool in the theoretical chemists’ repertoire 
has been employed for the study of this atom with the 
result that it is the best understood many-ele ctron 
system [96] . And the tour de force calculations that have 
been performed on beryllium [97,98] make it an ideal probe 
to test new theories and determine their range of validity. 
Above everything else, the spherical symmetry coupled with 
the presence of a mere four electrons renders it the most 
attractive choice for model studies. 
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III. 2 DETAILS DF COMPUTATIONAL METHOD 
111*2.1 Th e Atomic Q rbi tal s 


The occupied and virtual atomic orbitals of the 
beryllium atom are obtained as solutions to the self- 
consistent field equations of the Xa Hamiltonian. 


{- 


-n! 

2m 



-3ci [|^ P(r)]^^^} <J)^(r) 


8. 


i 


(^i(r) 


( III .1 ) 


In this study, the virtual orbital basis is choson to 
consist of the three 2p, the 3s and the throe 3p orbitals 
bounded by the self-consistent local potential of the 
occupied orbitals 1s and 2s. Slater’s exchange parameter a 
is chosen to be 0.768, keeping in view the observation, that 
orbitals obtained with this a in eqn. (III.1) closely 
mimic the corresponding HF orbitals [24]. The orbitals 
which are functions of coordinate space are written in the 
central field approximation as 


4). = R , (r.) Y (e,$) 

n.l* X l.m. 

XX J. J. * 


(III. 2) 


The R ,(r) form the radial part of the wavef unction and 
nl 

Yj ( 6 f (j:) STB the normalised spherical harmonics : n,l,rn^ 
are respectively the principal, orbital and magnetic 
quantum numbers* Wxthin the central "Field approximation, 
the radial part of eqn* (111*1) assumes the form, , 

[. ^ . V(x)] p„,(r) = P„,(r) ( 111 . 3 ) 

dr r 
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where , 




( I II .4) 


ore the normalised radial wavef unctions , The potential 
v(r) in eqn. (III. 3) is given by 

v(r) = - — + f / I [P^i(t)]2 dt + 2 j”l dt - 

0 nl r nl 

2 

3« -f — ^ I ^ *^nl^^^J > for r < r (III. 5) 

4it r nl ° 


and 

v(r) = = iLZz:~^.±l l 
r 


for r>r|^ (III. 6) 


Z being the atomic number and N the total numoer of 
electrons. Eqn. (III. 6) defines the Latter modification 
to the Xa potential [99], introduced to correct the un- 
physical behaviour of the local potential u(r) of eqn, 

(III. 3) at large r values. In the absence of Latter 
correction the virtual orbitals remain unbounded. A 
modified Herman-Skillman program [100] is used to obtain 
the solutions to eqn. (III. 3). In this program a linearly 
scaled 441-point mesh is used to obtain the numerical wave- 
function at these points. 

1 1 1 . 2 . 2 Matrix Elements and the Total Energy 

The total ground-state electronic energy of the atom 
is determined with the basis set of orbitals of 



78 


eqn. (III. 2) whose radial part is obtained as solutions to 
eqn. (111,3). The Fock Hamiltonian is constructed with 
these Xct orbitals for determining the total energy 


E = ^ 1(a) + I [J(a,b) - K(a,b)l (III. 7) 

3 apb 

pairs 


where , 

1(a) = / (f.g*(1) f(1) (!>^(1) dT^ (III. 8) 

where f( 1 ) defines the sum of one-electron kinetic energy 
and nucleus-electron attraction energy operator, and 

2 

J(a,b) = / J 'i»*(1) — ^h^2) dT dx (III. 9) 

" hi-2l 

2 

K(a,b) = / J (^/d) (!>b^(2) ^ ^ J2) d^(1) dx d x^ 

I ~^2 i 

' ^ (III. 10) 


are the Fock matrix elements that represent the coulomb and 
the exchange energies. The matrix elements J(a,b) and 
K(a,b) are given by the usual expressions [3,101] 


J(a,b) = I 3*^(3, b) F'^(a,b) (III. 11) 

k 

K(a,b) = 6(m ,m ) \ b^(a,b) G*^(a,b) (III. 12) 

®a ®b k 

k k 

where the factors a and b can be expressed in terms of the 
Wigner 3-j symbols 
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a, b) 


{m +m J 

(-1) ^ °(21^+1)(21^+1) 






( III. 14) 


k k 

and F and G are the Slater-Condon parameters. 


F^a,b) 


G ( a,b) 


e^R^(a,b| a,b) 


e^R*^(a,b; b,a) 


(III .15) 


(III. 16) 


The radial matrix elements R (a,b| c,d) are given by, 


R‘^(a,bj c,d) = / / 


k+1 

T 

D D > 


Pad) Pad) P^(2) dx, dr^ 


(III. 17) 


The general r.| 2 matrix eleiments are given by 

< a b t ■ 


2 (m +m.) ry 

— ^ |cd>= 6(m ,m ) 6 (m ,m )(- 1) ^ e 

r. -r-i I a c b d 


•1 "-2 


X I( " 

k \0 0 0 


X [(21 +1 )(2lj^+1 )( 21^+1 )( 21^+1 ) ] 

\ 

/ 


^r- Ad ^b Aa /^d ^b 


0 0 0 / \m^ ’^c r'^d "^b 


X R^(a,b; c,d) 


j where m =(tng )= (m^-m|^ ) . 


(III. 18 ) 
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k . • 

The F and G integrals are evaluated with a modified 
version of the computer program of Zars 102J and the ground” 
state energy of the atom is computed. Programs have also 
been developed for the evaluation of the general matrix 
element given by egn. (III. 18), An extended Simp son* s rule 
integration formula [103], is used For the evaluation of the 
r'^ integrals. A FORTRAN listing of the relevant program 
IS presented in Appendix B. All calculations presented in 
this thesis have been performed with a DEC 1090 system. 

III. 2. 3 The CC Equations 

The CC equations described in eqns. (11,28) and (11.30) 
are coupled non-linear equations, whose unknowns are the 
excitation coefficients, t's (vide Chapter II). These 
equations define the full CC many-electron theory, where the 
operator T contains T^ terms, T 2 terms and all the non-linear 
terms involving the products and squares of these. The non- 
linear CC equations are very complicated and tedious to solve 
As is true of any system of non-linear equations, existence 
of multiple minima can cause innumerable computational 
di Fficulties and lead to spurious solutions. Keeping this in 
view, we first attempted the solution of the linear part of 
the equations. 

The linear CC equations contain terms arising from T^ 
and T^ operators. The diagrammatic representation of these 
equations has been given in Figs, II, 5 and II. T ■ The one 
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open path and two open-path diagram equations independently 

but are coupled due to the presence of t”^ and t^® 
occurang in both the equations, eqns. (11.28) and (11.30).”" 

The linear part of these equations is given in eqns. (III.19) 
an(d { III.2D) . 


<r|f‘|s>+ I<r|f|a>t? . J <f|f|„ ,, 

o n 


I (2 <rp. |vjas> - < rp j v | s a > ) t® 


Gs 


+ I (2<B|f I s> t=,J - < 


p s 


T|f i s>)t;j' 


>u (2<rp v us> -<r,0 


(B|v! su>)t®^ 


p^g - <6P|0l sa > )t|J =0 (I II. 19) 


<rs|v| «p>+ I [<Pu|0!sr> t"' +<au 


V rs > to ] - 
I [ < Pa[0| > tj + < ac|^i > tf J + 

I [< r|f|us t- f<s|f|u> t”]- 

I [ <a|f 1 (S > tj® + <pjf I 6 > t®^] + 

1 [(2<xs [v[otu > - <r6 [v I u a > ) 
u (5 ' ' ' 6 . ' 


(2<s6 |v[ 3u> - < s6 |v| - 


<r(S|vj au > t®p - <5s|v|au> - 

<s« |v| Pu> -<6r|vlru> + 
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+ 72 ^ [<rs |vj uw > + <srlv{uw> t'^|^ ] + 

1/2 I [ <aclvl ^ e>t^® + < 3 ^j\?| 5 ^ ] 

^ e 

=0 ( III. 20) 

where a , f ... represent the occupied Is and 2s orbitals 

while r,s, ... represent the virtual 2p5 3s and 3p orbitals. 

III. 2. 4 Numerical Methods 

The linear equations, eqn. (III. 19) and (III. 20) 
contain, besides the unknown CC coefficients, the general 
2 matrix elements and the orbital energies, which have 
been evaluated as described in Section III. 2. 2. For our 

calculation on Be, we chose a basis consisting of the three 
2p, 3d and three 3p orbitals in addition to the ground -state 
Is and 2s orbitals. The number of cluster coefficients, i.e, 
unknowns of eqns. ( III .1 9) and (III. 20), total 210. These 
contain the 14 single excitation coef f icients, t^ and 196 
double excitation coefficients, t^p . The t^ p possess the 
permut ational symmetry t^^ = t^^^ and, therefore, the number 
of unknowns reduces further to 119. Our linear equations 
then are of the form, 

A X = B (III. 21) 

where A has the dimension 119x119, These equations are 
solved using the triangular decomposition procedure [l04] . 
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The matrix ^ is decomposed into an upper triangular and a 
lower triangular matrix such that 

A = L U ( III. 22) 

The system of eqns. (111,21), is then solved through the 
solution of the system, 


J= y 

= 1 

(III. 23) 

and then the system 




= 

(III. 24) 

since, _A x 


(III. 25) 


Since both L and U are triang ular ^ eqn s . (III. 23) and 
(111,24) are much more easily solved than eqn . (111.22), 

A routine to improve the solutions has also been introduced. 
The entire solution takes less than 10 secs of CPU time in 
the DEC 1090 computer. A listing of the computer program is 
presented in Appendix B. 

While the linear equations were well-behaved and yielded 

— »8 

solutions accurate to within 1 in 10” , the solution of the 
full non-linear equations turned out to be computationally 
very complicated. Several efficient algorithms [105] were 
used in an attempt to solve the non— linear equations. The 
Newton-Raphson procedure, which has quadratic convergence, 
yielded solutions which however corresponded to incorrect 
minima. We then used a modification of the Newton's method, 
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the powerful Brown's algorithm [106], with the hope of 
obtaining the correct minima. In Newton's iterative method, 
all members of the system of equations are expanded simul- 
taneously about a point, x^, assumed to be close to the 
solution. For the system 

•f(x) =0 ( III. 26) 

the expansion is, 

f(x) = f(x") -f J(x")(x-x") (111.27) 

where J is the Jacobian matrix; the higher order terms in 
the expansion are neglected. The expansion is iterated to 
convergence. In Brown's algorithm a successive substitu- 
tion scheme is used rather than the simultaneous treatment 
of all f^ as in Newton’s method; i.e., the modified root of 
f^ is used in f 2 and that of f 2 in f^^ etc. The derivative- 
free version, where the derivatives in the Jacobian were 
replaced by first-difference quotients, was used. The 
search based on Brown's method also failed to converge. 
Further work is in progress. Only results obtained from the 
solution of the linear equations are reported in this thesis. 

III. 3 HIERARCHY OF APPROXIMATIONS TO THE CC EQUATIONS 

A. In this approximation the entire linear part of the 
cluster equations, eqns. (III. 19) and (III. 20) is set-up and 
solved. These equations result from inclusion of all the 
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diagrams involving the linear and terms. This 
procedure is formally equivalent to performing a Cl cal- 
culation with all single- and double- excited determinants. 

B. In this set of approximations only the terms arising 
from the linear part of T 2 diagrams have been retained in 
the CC equation. In the first case one and two open-path 
diagrams of the above mentioned type (diagrams F-K of 
Fig. II. 5 and diagrams A and D-K of Fig. II. 7) have been 
retained, while in the second case only two open-path 
diagrams of the above mentioned type have been retained 
(diagrams A and D-K of Fig, II. 7). Both cases lead to 
almost identical results and, therefore, the results of 
only one of these has been presented. The procedure is 
equivalent to a variant of the coupled electron pair 
approximation [ 31] . 

C. In this approximation all the one and two open-path 
ladder diagrams (diagrams A-C of Fig. II. 5 and diagrams 
A,D,E and J of Fig. II. 7) arising from the T^ and T 2 
approximation have been included in the CC equation. The 
procedure is equivalent to solving the generalised Nesbet- 
Bethe-Goldstone equations [46]. 

D. In this approximation, only the terms arising from the 
two open-path ladder diagrams (diagrams A,D,E and J of 
Fig. II. 7) have been included in the CC equation. This 
procedure is equivalent to a partial summation of ladder 
diagrams to all orders in MBPT [^2]. 



E. In this approximation only the particle-particle 
ladder diagram given by diagram J of Fig. II. 7 is included. 
This result is equivalent to a selective summation of the 
corresponding particle-particle ladder diagram in pertur- 
bation theory to all orders. 

F. In this approximation , the CC equivalent of the 
second-order R5PT equations are set up and solved. This is 
achieved by taking into consideration only the diagram A of 
Fig. II. 7 into account for generating the terms of the 
equation . 

III. 4 RESULTS OF THE CALCULATION AND DISCUSSION 

Table III.1 presents the pair-wise contribution to 

the energy corrections and the total EC energy corrections 

for the Approximations A to F mentioned above. In Table 

I II. 2, we present a partial wave analysis of the EC energy 

corrections wherein the ss and pp excitation contribu— 

2 2 

tions to Is , Is 2s and 2s pairs are determined separately. 
In Table III, 3, we give the individual contributions to the 
EC corrections, 

A cursory analysis of the results presented in 
Tables III.1 and III. 2, reveals the following : In all the 
approximations used, the contribution to EC correction ener- 
gies comes mainly from the 2s^ pair, whereas in the calcula- 
tions of Paldus et al [90], where the HF method is used 
as the starting point, 1 s^ pair contributions are found to be 



Table III.1 EC Energy Corrections in Beryllium Atom 


H 

ro 

-p 

o 


+ 

>5 

cn 

M 

0 

C 

LU 


* 

C 

o 

-P 

u 

0 

p 

p 

o 

u 


0 

JC 

•p 

o 

p 

c * 
o c 

p Oj 
P -ril 

n p 
M o 
•H ml 
p p 
p P 
c o 

0 CJ 
tJ 

0 
CO 

ri 

5 

1 

P 
H 
CD 

a. 


(M 

to 

CM 


to 

CM 


CM 


CO 


C 

o 

•H 

P 

ro 

E 

X 

Q 

P 

a 

a 

<c 



0 

0 

P 

p 

p 

ro 

X 

O 

U) 

p 

C 

a 

c 

•H 

to 

0 

□ 

H 

CD 

> 


+ 


a f 


Table III»2 Partial Wave Analysis of the EC Energy Corrections 
in Beryllium Atom 
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Table III. 3 Individual Contributions to EC 
Corrections in Beryllium Atom 


■^virtual 

occupied 

„ 2 

2p 

2p 3p 

to 

T 2 

3p 

1 

-8.1063 

2.5304 

-0.5343 

1 .6906 

1 s 2s 

-31 6.0383 

134.4522 

2.5562 

64.8626 ' 

2 

2s 

-16263.1 236 

4920.681 2 

-244.8591 

1956.9266 


“5 

All values in units of 10 Hartree. 
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equally important. The EC effects seem to be well represen- 
2 

ted in the 1s orbital by the X® method. Tong and Sham [107] 

have observed that the uniform gas model leads to over- 

estimation of correlation energy and the under-estimation of 

exchange energy in the inner regions. Our results seem to 

indicate that these two opposing effects balance each other, 

with the net result that EC corrections in the inner orbital 

2 

are negligible. The 2s orbital, on the other hand, seems 

to be poorly represented by the Xa model. It is a well- 

known fact, that valence regions in atoms are not correctly 

represented by this model [75], The local-density potential, 

therefore, appears to be well-suited for mimicing the exact 

electron-electron interaction potential in the inner regions 

of atoms only. Tne partial wave analysis given in Table 

III. 2 reveals that the dominant contribution to the energy 

corrections come from the pp excitations, amongst these the 

2 2 

most important are the 2s -»■ 2p excitations (see Table III. 3), 
This is due to the near-degeneracy of the 2p orbital, and its 
spatial distribution nearly parallelling that of the 2s 
distribution . 

Table III.1, also reveals the fact that T^ terms play 
a crucial role in the evaluation of EC corrections. In the 
calculation of Paldus et al. [ 90 ] , with linear coupled-pair 
many-electron approximation, only the T 2 diagrams with two 
open-paths were considered. By virtue of , ^ using the HF 
basis which obeys the Brillouin's theorem, they could 
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neglect all the single excitations. However, our calcula- 
tion uses the Xot orbitals as the basis; these orbitals do 
not obey Brillouin*s theorem and the excitations survive. 
Also , the matrix elements of the type <i|f|j> are non- 
vanishing, since Xa orbitals are used in the calculation. 

The terms however contribute only indirectly through 
their coupling with the terms. The approximations which 
neglect the T^ terms tend to over-estimate the EC corrections. 
This is similar to the results obtained by Paldus et al [ 9 D] 
where the quasi-degeneracy of the 2 s, 2 p orbitals has been 
shown to cause over-estimation of correlation energies. 

An interesting feature apparent from our calculations, is 
that contributions from ring diagrams are found to be very 
important for a proper estimation of the EC corrections. 
Approximations C and D, where only contributions of ladder 
diagrams are considered, tend to over-estimate the energy 
corrections. The importance of ring diagrams in represen- 
tation of EC • in electron gas approximations is well known 
[71,95] and is reflected in our calculations as well. For 
the purpose of comparison, we performed the second-order 
RSPT calculation (vide Approximation F) . This accounts for 
only about 50 ^ of the corrections. In tho following chapter, 
we present a detailed discussion on these and other results 
obtained by us in our CC analysis of the EC effepts. 
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IV. 1 GENERAL 

This chapter presents the analysis of the many-body 
EC effects in beryllium atom and the ions B"*", and 

N , which are isoelectronic with Be, The CC formalism 
and the hierarchy of approximations derived from it have 
been employed for this purpose. The Xo local potential 
model has been used to generate the reference functions for 
the CC expansion. The comput atio.nal methods employed have 
been already outlined in the previous chapter. The analysis 
presented in this chapter entails the study of the EC 
effects, for various choices of a in the Xa potential, with 
gradient corrections to the X“ potential and for Z-variation 
along the isoelectronic series. The results are discussed 
in the light of the electron gas approximation which serves 
as the basis of the Xa model. As has been shown through 
variational arguments [15] , the use of the uniform electron 
gas potential in the model leads to a valud of 2/3 fora • 
Use of any different value of a in the Xa potential would 
imply deviation from uniform charge description. In view 
of this, the role of a in the portrayal of EC effects, is 
examined. The use of gradient corrections to account for 
the inhomogeneities in the charge distribution is also 
examined from this view point. 

After commenting briefly on the CC equations in the 
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next SBctian, we present the results of our calculations 
in Sec* IV.S. The results are discussed in Sec, I\/,4 after 
a brief resume of the theoretical framework concerning 
many-body effects in electron gas models# The conclusions 
are then summarized in the background of this theoretical 
presentation # 

IV. 2 SOME COMMENTS ON THE CC EQUATIONS 

The CC equations, eqhs# (III. 19) and (111.20) in 
full, and under a hierarchy of approximations, have been 
solved for the Be isoelectronic series using the numerical 
methods outlined in Chapter III. The structure of the 
equations, eqns. (III.19) and (III#20) becomes clear if they 
are recast in the following manner : 

<cclflr> + e^)tj + A tj + = 0 (IV. 1) 

and 

<aP|vlrs> + (£„ +ep~e^ - + K 2 = 0 (IV. 2) 

Eqn. (IV, 1) corresponds to eqn. (III. 19) where A is the 
coefficient oF the diagonal term, t , and stands for 
all the non-diagonal terms with their coefficients. Like- 
wise, eqn. (IV. 2) corresponds to eqn. (III. 20) where B is 
the coefficient of the diagonal term, t^ g, and K 2 stands 
for all the non-diagonal terms with their coefficients. 
A,B,K^ and K 2 contain the hole-hole, hole-particle and 
particle-particle ladder and ring diagram terms(uide 
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Figs. II. 5 and II. 7). These diagrams represent hole state 
interactions occuring via creation of particle-hole pairs, 
and interactions arising from scattering between hole-hole, 
hole-particle and particle-particle states. In Brueckner's 
electron gas theory [71,72] ring diagrams contribute domi- 
nantly to the correlation energy correction while in 
Sinanoglu's many-electron theory [41] ladder diagrams 
contribute dominantly. In the work of Singal and Das [ 74] 
the non-diagonal terms were replaced by a constant times 
the diagonal coefficient and the resulting equations were 
solved. The importance of ring diagrams was again gleaned 
from their calculation. In our calculations, we have solved 
the linear part of eqns. (IV.I) and {l\l.2) with no approxi- 
mation to the non-diagonal terms. The effect of various 
diagrams on EC corrections is studied through the various 
Approximations (vide Chapter III), For instance, Approxima- 
tion B neglects the T^ excitations altogether, while Approxi— 
m'ations C, D and E neglect ring diagrams entirely and include 
some of the ladder diagrams. A comparison is made with the 
second-order RSPT calculation made by us to highlight the 
importance of higher order corrections. 

IV. 3 RESULTS 

Tables IV.I to IV, 3 present the results of pair EC 
corrections in beryllium atom for the choice of various a 
values and under the different Approximations outlined. Total 
EC corrections and the calculated total energy values for 



Table IV.I 1 Pair EC Corrections in Be - Effect of a Variation 


CO 

MO 

a ^ 

ch 

ro 

CO 

(T\ cn 

r-* 

cn 

cn T- 

O 


'XD 

CD 

a MO 

r- 

^ VD 

cn 

t- 

r- 

m 


VO 

Q 


CJN 

OD VO 

tn 

tX) 

m 

VO 


o\ 

CM 

in cn 


CD 

CTN 

a\ 

a 

C\J 

CD 


ID VO 

0- 

C^ 

O <n 

VD 

a\ 

* 

* 

« « 

« 

» 

• 

9 

• 

• 

• 

• • 

• 

9 


X— 

2 

2 

OJ 

C\i 

oo 

OJ 

OJ 

CM 

OJ 

cn cn 

cn 

m 

r 

1 

1 1 

1 

1 

1 

1 

1 

1 

1 

i I 


9 6 


UJ 


h—r— VOCD ON 

cr^a^^^Occl r"C^c^Jt-•vovoc^Jc\Jh-•v£) 

t^OD CJ\r-'^T-aC\l 

cDChaO't-tMm^invoGDcrvfnLnco 

• • •••••«»«««••• 

•r-'-r- CMc\jf\jr\)C\Jc\ic\Jc\ic\jc\irorofn 

1 t I t I I I t I t I I I I I 


c 


OJ 

a 

LO 

cn 

VO CM r- 


cn 

CD in ^ 

VD 



o 


in 

CM 

ro 

VO cn 

^ (M ON 

ON CM 

VO 

CM m «r“ 

ir* 



•H 

n 

m VO 

03 

-t- in 

QD ^ m 

VO OJ 

o> 

h- m 

cn 



+> 


0\ OJ 

OJ 

in ^ 

CO VO VD 

03 CM 

VO 

CM VO CM 

CO 



m 


• 

» 

• 

♦ • 

« « » 

• « 

9 

9 » 9 

» 



B 


LO CM ov m cn 

'T- a o\ 

CD CO 

r- 

7 

6 

6 

6 





CM 

CM 



r** 







X 


1 

1 

1 

t 1 

1 1 t 

I 1 

1 

1 1 1 

1 



o 













u 


























Q 













< 
















VO CM ■n* VO 

CD CD 

0\ CD 

VO 

CO a c?v 







o\ 

T-' 

CO 0\ 

o- in a 

CD VD 

CD 

in OJ cn 

cn 


« 






a\ OJ 

cn cr\ o 

cn <d- 

On 

CM VD VO 

in 

Q) 

C 


C-J 

* 

a cn ^ 

h" GO CO 

ON CM 

VD 

CM iD 'T- 

cr\ 

03 

o 




• 

« 

» t 

* • « 

* • 

• 

* • f 

9 

U 

•H 




O CO VO ^ 

CM o CO 

CD CD 

r-* 

7 

6 

6 

in 

-P 

-P 




CM 








M 

3 




1 

1 

1 1 

1 1 i 

t 1 

1 

1 ( 1 

1 

OJ 

H 












X 

O 













CO 












in 













t 

03 












a 

JZ 






Q 

crv CD O 

VD 

o 

in h- VO 

VD 

IT- 

-p 






CO 

a in ^ 

in OJ 

OJ 

OJ CM in 










^ 

'M' m 


'P-* T- 

o 

cp 

c 



* 

* 

* 

ID 

r-- cn -r- 

0\ CD 

CD 

CO CJv 


O 







VD m 

^ ^ 

cn cn 

cn 

m m 

<d- 

CO 

n 






1 1 

1 1 1 

1 I 

1 

1 1 1 

1 1 

+> 

0) 











1 

•H 

M 












c 

0) 












3 

-P 













c 












C 

3 












'H 

O 






cn in 

VO a -f- 

CO 

r- 

CM VO 

in 


u 






a VO 

0“ CM o 

CM in 

cn 

^ -r- 

OJ 

CO 

c 






VO ^ 

VD m CO 

'r- •n' 

in 


c7 

D 

03 


< 



3k 

CD cn 

0- O VO 

in ^ 

'n- 

in h-* 

cn 

3 







• « 

9 9 9 

• • 

• 

» • 

* 

H 

(0 






O VO 

^ Kt cn 

cn cn 

m 

3 

3 


CD 

03 












> 

•H 






1 1 

f \ \ 

1 1 

1 

t I 

1 


-P 











>> 

-H 


0- 

C3 

C3 

a 

O 

a 

O 

o 

O 

a 

a 

o 

O 

a 

o 

VO 

a 

O 

a 

a 

o 

o 

o 

C3 

o 

a 

o 

o 

a 

o 

VO 

X 

a 

CM 


VD 

X 

a 

CM 


VO 

CO 

OJ 

VO 

o 

VD 

VO 

0- 




0- 

CD 

CD 


X 

CD 

A 

CTv 

C3V 

• 

a 

* 

0, 

0, 

0. 

0. 

a 

0. 

(3 

0. 

O 

a 

a 

0 

O 

a 



m H 

U 

m xj 

C CD 

0) 

w 

H C 

H 

<C * 


91 


c 

o 

•H 

4J 

m 

•H 

(D 

> 


<3 


tf. 

O 


voMDc%joQcocNJoor\i'^r\jc\icDrvj 
cn'^c'i^ococvjor;)-N:3'crih“VO<d‘CDC\J 
^tnfnoi\sD-*-aschOsJcor*^cD acMoo 
r^CTxfnr^^MDQLD^MDCNJCDCNJVfDr- 


QQr-'ir-(\jt\jron<:f«stLnLnh-ooa 

1 > t I I I I I t I I I I I I 


CO^f\i(^JC^JOOCD(^ICO^:^'OC)*=d''^0 
Ln(T\oroco-=;tincT\CDCDONoa'^Ln 
'^aD^r-aro'vO(MT-rnr^*ij'Lnv4Dh~ 
crNv-Lncr\rQ vjd 


c^aoO’r-‘-e-(Moj(^(D^'^Lnr— 00 

I -r— 

{ 1 » I ) 1 1 I I I I i I I 


-P 

U 

tu aaocoooDoocovQ^^^iMoo 

^ cr\PJ'vOcoc\iar^ocoaDCD’sj‘c^rO‘^ 

^ OOcDvocNjrvJvocNJv^rocoroinvoh- 

LU Q v-ir)VX3C3\4DO\CncJ\t^Cnr-“a\CD«r-r-* 


0 

m 

c 

•H 

0 

C 

o 

•H 

■P 

U 

0 

P 

P 

O 

u 


c 

o 

*p 

4 ;) 

0 

e 

‘H 

Xi 
□ I 

I Pl 



<c u 


LJ 

UJ 


p 

m 

CL 


m 


0} 

00 



rsj 

> < 
0 

rP 

JO 

0 

H- 


o>roLnr*“Lnr^crvh-T-a\OPO'^ tnco 
cncnLnroh-OvJcotncncDa\r^*^CNJo 

I^MDU0<:3"CO(^r\Jf\J(\irvJ«r«*r--T-T--T-« 

I 1 I I I i I I I I I I I t I 


C}OJCDCD-^r\Ja^'^'^CDcr3'J3VD 
rnr-rjcTvONtDrocnr^cTxincDC^i^ 
cn r^cor^<ov 430 NC\ico'r-c\a^'r-h- 
v-oc^h-rnar^cnvornoc\i\socrj 

iooMncnr-^<^j-inaarMcoT-^ 

vOT-h-otno’^fn^crvi^’^CNJo 

tnLn^^nrof\jrvjr\j^T-T-T-T- 

I 1 t t I t I I I I t I I I 


OJOCDC\iOa'vOC\JM3'sOVOMD 
CD r^(MLna^cDc\jN-r-ovD 
LnrJfO’5d'l^'^Ln'JorncD*^in 
pjrnh-cT\oa-^o4cr\aoT^o 


T-T-CrNr^OsJO^OCDCJNCMOD 

^vDOJooNCDl^^Ln^'sj-cn 

OJr-v-'!- 

I I I I I I 1 I I I I i 


OCMVDVD'^CviCDCDOJCO O 

(NJCTNO^^^r^COLnCO CD 

^^rvjcncovx)0(Msj*'«t co 

3K j<c v£)t-ofO-r-incr\0^un'?J*co 


VOOCT\t^CDin'JDr“h-«r- LD 

uicTN(>JaNr^vouoin*=d'in v- 

rOr-v- 

I I I I I I 1 I i I I 


h-aoooaoaaooaaoa 

vDaaaaaaoooooooa 

\£)MD h-?^f^h"l^aDCDCOCOaDa\CT\0 


oooooaoaaoaoaaT- 


0 

c 

0 

Q 

p 

•H 

•p 

P 

p 

n 

0 

H 

X 

0 

■> 

0 

1 

ai 

0 

JZ 


p 

«P 

c 

0 

•d 

0 

-0 

’P 

0 


p 

C 

0 

D 

p 


c 

C 

3 

•H 

0 


U 

0 

c 

0 

0 

3 


H 

(f) 

0 

0 

> 

•H 


P 


'H 

CJ> 

P 

P 

P 

0 

X 

c 

0 

0 

P 


0 

rp 

c 

rH 

l-H 

< 

ij: 


Table IV. 3 2=^ Pair EC Corrections in Be - Effect of a Variation 


tlJ 


a 


u 




'«^“aDCT^c\^o^^^-oa^ OLnvo cnoco 

^c\]c\j«::j'in-<-inc\JLncMcncD^"Ofn 

acoa3h-io^cD^roaDaQDCr\f^cN4 

O^^U^Ln^-•^-'^0■;d•^VDx-vxDL^lT)C7^ 

CDQNr\JLnoDf\iLncr\fr)h*-f\ivoMDh-QD 

rsirvjnrnfo^^-'^Ln!r)VD\oh“CDcr\ 

I I I t I I i I I I t I i I 1 


cr\l^a\T-rr)co h“LnrQcDcr\cr)^'r-MD 
-r-co^rn-r-r^cocrvvorn'^vot^x—'v}' 
cot^a:jinr^LnaDh-o\ir)Or-“OvO' 
CQVOCDh-C\lC^^V£)t^CD<:J^COf^CrN 

r^Lnl^o*^cDrncO'^T-cnN-Lnv.ocr> 
r^oDONT-cMmLnvDOD OT-roh-'r-Ln 
cocDCDo^□^o^a^c^\□^aooa^‘ 

<r*“ ▼*“ “v* x“ T“ c\j cNi c\i (M rsj rvj 

I t I t I t t I < I 1 I 1 t ( 


oaaooaoof^J-f~c^^-o^vD'X} 

OMDLnOSCDaDO'^ir)*d‘t^'<“f^’r*-rM 

o\vDh-T“inir)CM^Oooo’^crsco 

h-'^^c^^co(^JT-'^MDh"0^■«— oo^OO 

v£)^-'^"aCT^^-^*f^0D'vD^~CDCDvO^ 

rNJ'^M3’^T-^OcrNCJ\’c-‘^coa:i-»~*^ 

co'd-'«-r--LncncNJoa\cr\CDr^vDVDLn 

I ( t t I I I I I t I I t ^ I 


aaooaooor-fn^vouDcn 

c\Jo\*=^rnvi3h--t^ON«t--(:nvD'vO^^ 

fOCTNlDrvJCMLnVOOCDT-VD’r-fMO 

aD-c-CT^acD^^-cocDLna^^-^^^ 
?{i 5--uDO^cnaQo^^f^coacr\M3 
cof^coOJ'sDu^T-r-(^J’^^-^“CO^^^ 
^OCDVX}^OvJ-r-OC^COh-MDLnLO 

I t « I I I I I I 1 » I » * 


aooao^coaf^^^^'-o^" 

oacQcnc^co^J^t^ 

Ln^^iD'sDiD'T-h-r^cncTsO 
T-^h^oorO't-r^cD'-ocoh- 

VX3COCna-v-‘Cr\f^CDCTN'r--'r 

co-^-ooc^]oo^o^-c^JLnc3^<^^‘ 

vncoocrvoor^'^LDin^rom 


:♦« i|« * 


I I t t 1 


t I I 1 t t 


oooLnoJ^ar^ 

Ln**^covacrNi/iLnh- 

LOh-aoN^omr-'^ 

cocD'^OJrncnLnir) 

»••••••« 

9|C 5<C VX3COrOCMr^^rM<^ 

C3\LnoNOt^’^cn(\j 

^mr^rn h-h-OLn 
coLoaco^tnin^ 
<M ^ T- 1 

I I I I I I I I 


ro n 

CO 

ro Osi 

C\J 

in vrj 

0^ 

« » 

• 

a <M 

^ a 

a\ a 

r\4 

•r* 

CM 


r- 


I I 


t^aoaoaoaoQoaooo 
uDoaaaooooooooag 
v£)Sarvj'=a-vncDa(M^votDr\ju30 
\r\ f*” □□ CD CD CO CO 0^ O^ O 

aooooaoooooooo»- 


CD 

c 

CO 

0 


'H 

4^ 

-P 

U 

0 

m 

H 

X 

0 


(0 

n 


f 

0} 

CD 

X 

T— 

-p 


c 

0 

•H 

to 

T» 

-P 

QJ 

•H 

M 

C 

(D 

D 

-P 


C 

C 

0 

•H 

0 


u 

to 

c 

0) 

Q) 

D 


H 

to 

CD 

(U 

> 

•H 


-P 

>5 

*H 

cn 

H 

M 


CD 

X 

C 

to 

CD 

-P 


w 

H 

c 



<s: 



98 



99 


various a’ s are presented in Table IV. 4 and Table IV. 5. 

Similar results for B"*", and are presented in 

Tables IV. 6, IV. 7 and IV. 8 respectively. The EC corrections 
obtained with the gradient-corrected X® potentials for a 
pair of a values, viz. 0.768 and 0.6667 in the case of Be 
atom are given in Tables IV. 9 and IV. 10. 

IV. 4 DISCUSSION 

I V . 4 . 1 X and the Electron Gas Model s An Analysis 

The basic tenet of the Xa and like methods is the 
approximation of the electron density distribution in atoms 
by the uniform charge density of an electron gas. In a boms 
the electron density distribution is non-uniform and the 
validity of the above approximations is, therefore, subject 
to question. However, the electron gas model can be extended 
to non-uniform charge distributions, as well, provided the 
density is either high [71] or slowly varying[72]. This is true 
for the bulk of the electron distribution in atoms. Regions 
which are neither too near the nucleus nor far removed can 
be described as possessing a slowly varying charge distribution 
and the electron gas approximation would have reasonable vali- 
dity in these regions. In regions very close to the nucleus 
the electron density is high, but the density gradient is 
large and the changes are rapid. In this region of large nega- 
tive potential energy, the electron gas approximation leads to 
spurious densities, thus causing an unphysical behaviour of the 
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Table I\/,6 Pair-wise and Total EC Corrections 

and Total Sround State Energies in B'*’ 


Con tri- 

ct 


ADProximation* 


button 


c 

D 

E 

F 

1=2 
p air+ 

0 , 6667 
0,7200 
0.7600 
0.8200 

1 .0000 

-24.2430 
-22.7720 
-22.1115 
-21 .4717 
-21 .1873 

-23.7216 
-22.6367 
-22.0534 
-21 .4796 
-21 .2350 

-1 2.231 6 
-1 2.9635 
-13.5459 
-14.4633 
-1 7.3843 

-12.7395 
-13.5037 
-14.1 1 24 
-15.0719 
-18.1332 

1 s 2s 
p ai r"^ 

0.6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-31 5.2612 
-270.5920 
-255.2404 
-225.5220 
-162,4568 

-307.2124 

- 274.7944 

-253.5426 
-225.7670 
-1 64.11 72 

-47.2238 
-49.2000 
-50.8184 
-53.4224 
-62.11 76 

-51 .1420 
-53.2940 
-55 .051 2 
-57.8912 
-67.4096 

•5 2 

2s 

pair"'" 

0.6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-7356.0910 
-6750.1579 
-6357.8803 
-5852.1823 
-4 764.1 236 

-7151 .8010 
-6647.6688 
-6309.7828 
-5859.6859 
-4852.1684 

-2981 .7319 
-3010.3846 
-3033.9545 
-3072.5932 
-3210.5439 

-6771 . 9813 
-6869.8808 
-6950.0663 
-7080.4938 
-7536.3983 

Total 

Corre- 

lation 

Bnergy+ 

0 ,6667 

0.7200 

0.7600 

0.8200 

1 ,0000 

-7695.3764 
-7051 .5219 
-6635 . 2321 
-6099.1 759 
-4967.7677 

-7482.7350 

-6945.0998 

-6565.3787 

-6106.9325 

-5037.5205 

-3041 .1873 
-3072.5524 
-3098.3189 
-3140.4789 
-3290.0548 

-6835.8627 
-6936.6784 
-7019.2298 
- 71 53.4568 
-7621 .9410 

Total 

energy 

in 

Hartree 

0.6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-24.299134 

-24.297950 

-24.296782 

-24.294232 

-24.278033 

-24.297007 

-24.296886 

-24.296284 

-24.294309 

-24.278730 

-24.252592 
-24.258161 
-24.26141 3 
-24.264645 
-24.261 256 

-24.290539 

-24.296802 

-24.300622 

-24.304775 

-24.304574 


*In Approximations 
solution 

"^All energy values 


A and B instabilities encountered 

-5 

in units of ID Hartree. 


in 


the 
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Table 1\J ,7 Pair-wise and Total EC Corrections „ 
and Total Ground-State Energies in 


Contri- 

a 


AoDroximati on 


button 


C 

D 

E 

F 

1 2 

1 s 

pair'*' 

0,6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-38.9484 

-36.3067 

-37.9967 

-37.7481 

-37.9537 

-38.7947 

-38.2308 

-37.9603 

-37.7498 

-37.9930 

-26 .4460 
-27.4584 
-28.2444 
-29.4475 
-33.0034 

-27.4872 
-28.5407 
-29.3591 
-30.61 26 
-34,3225 

1 s 2s 
pair+ 

0.6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-334.6036 
-308.5242 
-290.9968 
-267.4646 
-21 1 .4250 

-330.7684 
-306.5632 
-290.0350 
-267.5098 
-21 2.4394 

-85.2072 
-87.3006 
-88 .991 6 
-91 .6882 
-100.3378 

-91 .8206 
-94.0788 
-95,9062 
-98.8252 
-108.21 28 

o 2 

2s 

pair+ 

0.6667 

0.7200 

0.7600 

0.8200 

1 .0000 

-6569.9495 
-6273.5244 
-6071 ,9774 
-5800.2808 
-5179.1428 

-6480.7426 
-6224.8707 
-6046.8764 
-5801 .5575 
-5217.4015 

-3778.8860 
-3808,4640 
-3832.6371 
-3872.2675 
-401 1 .1085 

-8478.3847 
-8598.5605 
-8695.4047 
-8851 .382'5 
-9372.6277 

Total 

corre- 

lation 

energy"*" 

0.6667 

0.7200 

0.7600 

0 .8200 

1 .0000 

-6943.5014 

-6620.3550 

-6400.9709 

-6105.4934 

-5428.5215 

-6850,3057 
— 6569 .6645 
-6374.8716 
-6106.81 70 
-5467.8338 

-3890.5391 

-3923.2230 

-3949.8733 

-3993.4032 

-4144.4496 

-8597.6926 
-8721 .1798 
-8820.6700 
-8980.8207 
-951 5.1630 

Total 

energy 

in 

Har- 

tree 

0.6667 

0.7200 

0.7600 

0.6200 

1 .0000 

-36.462605 

-36.464643 

-36.465460 

-36.465360 

-36.453685 

-36.461673 
-36.4641 37 
-36.465199 
-36.465373 
-36.454079 

-36.432076 

-36.437672 

-36.440949 

-36.444239 

-36.440845 

-36.479147 
-36,485651 
-36.489657 
-36.4941 1 3 
-36,494552 


*In Approximations A and B instabilities encountered in 


the solution 


■‘‘All 


energy values in 


units of 


10 


-5 


Hartree , 
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Table IV.S Pair-wise and Total EC Corrections 

and Total Ground-State Energies in N+3 


Contri- 



AcDroximation 


bution 

Ui 

C 

D 

E 

F 


0,6667 

-53.4573 

-53,3448 

-40.5062 

-41 .9643 

1 2 

0.7200 

-53.0404 

-52.9797 

-41 .621 6 

-43.1 201 

•1 s 

0.7600 

-52.8562 

-52.8236 

-42.4780 

-44.0079 

pair"^ 

0.8200 

-52.7432 

-52.7396 

-43.7732 

-45.351 5 


1 .0000 

-53.0161 

-53.0502 

-47.4550 

-49 .1 742 


0.6667 

-362.9558 

-360.2948 

-116.7482 

-1 24.9784 

1 s 2s 

0.7200 

-340.1164 

-338,6610 

-118.6924 

-1 27.0574 

pair+ 

0.7600 

-324.3942 

-323.6072 

-120.2622 

-128.7394 


0.8200 

-302.8232 

-302.7378 

-122.7762 

-1 31 .4266 


1 .0000 

-248.6756 

-249.3746 

-130.6930 

-139.9506 


0.6667 

-6717.2496 

-6655.0769 

-4504.8491 

-1001 2.7630 


0.7200 

-6512.0927 

-6475.7992 

-4537.9332 

-10157.5890 

o 2 

2s 

pair+ 

0.7600 

-6370.5758 

-6349.9286 

-4564.6477 

-10272.9550 

0.8200 

-6177.3783 

-6174.9527 

-4607.941 2 

-10456.5670 

1 .0000 

-5728.9502 

-5755.2325 

-4757,6843 

-11060,9780 


0.6667 

-7133.6626 

-7068.7165 

-4662 .1036 

-10179.7050 

Total 

0.7200 

-6905.2496 

-6867.4398 

-4698.2472 

-10327.7660 

Corre- 

0.7600 

-6747.8262 

-6726.3593 

-4727.3879 

-10445.7030 

lation 

0.8200 

-6532.9445 

-6530.4300 

-4774.4806 

-10633.3450 

energy^ 

1 .0000 

-6030.6421 

-6057.6573 

-4935,8323 

-1 1 250.1030 


0.6667 

-51 .1 38337 

-51 .1 37687 

-51 .1 1 3621 

-51 .1 687970 

Total 

0.7200 

-51 .141303 

-51 .140924 

-51 .119232 

-51 .1 755280 

energy 

0.7600 

-51.142778 

-51 .142564 

-51 .1 22574 

-51 .1797570 

in 

0,8200 

-51 .143579 

-51 .143554 

-51 .125995 

-51 .1845830 

Har- 

tree 

1 .0000 

-51 .133556 

-51 .133826 

-51 .122608 

-51 .1857510 


*In Approximations A and B instabilities encountered in 
the solution 


+ 


All 


energy values in 10 Hartree, 
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distribution [100], Nonetheless, this is not a serious 
shortcoming, tor the kinetic energy term dominates the 
i nt er- ele ctron interaction term close to the nucleus end 
if the kinetic energy as treated exactly as is done in the 
X ct method, the use of electron gas-like local potential 
approximation is reasonable. However, this justification 
is of no avail at large distances from the nucleus. In 
these regions the electron densaty is low and the inter- 
clectron term becomes significant in comparison to the 
kinetic energy term and the electron gas approximation breaks 
down. Besides, the use of the electron gas approximation 
entails two serious shortcomings. Firstly, the uniform 
charge description has the effect of smoothening the den- 
sity variations throughout the atom and therefore is incapa- 
ble of depicting correctly the atomic-shell structure. Also, 
the approximation fails to portray correctly the effect of 
the long range of the coulomb interaction. This is due to 
its inability to account for the shielding of the intera- 
ctions between the particles [109]. 

Several corrections to the electron gas approximation 
have been made in order to take into account the non- 
uniformity of the charge distributions in atoms [14]. The 
Xct method itself is an outcome of this. Use of any a value 
other than 2/3 in the Xa potential takes it beyond the 
electron gas approximation [11D]; values of a higher than 
the Gaspar-Kohn-Sham (GK5) limit of 2/3, characterise more 
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non-uniformity of the charge distribution. For systems 
with higher Z, better description is obtained with a values 
nearer to 2/3, As Z increases, the high density region in 
the atom expands and the uniform charge approximation finds 
greater justification [111] , In the limit of high Z the 
uniform charge distribution GKS value of n = 2/3 is obtained. 
Although, higher values of a may account for non-uniformity 
of charge distribution there are no rigorous theoretical 
grounds for choosing a , Due to the absence of a varia- 
tional bound on a , the energy decreases with increasing a. 
Several prescriptions have been given for the choice of a , 
the main amongst these being the choice for which the Xa 
energy reaches the HF limit or the total energy limit [20] . 

The resulting self-consistent potential is an EC potential 
although a demarcation of the “exchange-only*’ and “corre- 
lation-only” parbs cannot be made. 

The electron gas system most relevant to atomic pro- 
blems is the high density electron gas subject to a uniform 
positive field. This has been studied by Gell-Manp and 
Brueckner [71] through MBPT techniques. The uniform electron 
gas study gives very inaccurate estimates for atomic corre- 
lation energies [112]. To account for the discrepancy, the 
effect of non-uniform charge density was treated as a 
perturbation to high density electron gas by Ma and Brueckner 
[ 72] . Their analysis was based on the treatment of the 
electron-electron interaction to all orders in the ring 
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diagrams and the external field producing the non-uniform 
density to second-order. The ring diagrams yield terms 
which take into account the coulombic screeninn and are 
therefore extremely important for non-uniform charge dis- 
tributions. The slow density variation was accounted for 
by expansion of the EC potential in powers of density 
gradient. For atoms where regions near the nucleus with 
high density gradients exist, this expansion leads to 
infinite energies. The gradient-corrected EC potential 
also fails to account for the sum rule [69], which is a 
very stringent criterion on EC potentials, Brueckner has 
also shown, that at intermediate electron densities the 
perturbation theories diverge, and the electron gas approxi- 
mation becomes unacceptable. An extremely interesting 
result of Brueckner' s electron gas study is the observation 
that the ring diagram terms are the most divergent terms 
in tha perturbation series. This is in contrast to the many— 
electron theory result of Sinanoglu [41], which discusses the 
relative unimportance of the third and higher order rings. 

Our model which is a superior version of Sinanoglu' s many- 
electron theory has overtones of Brueckner's theory by virtue 
of the use of electron gas-like local potential approxima- 
tions to provide the reference functions. In the light of 
this it would be worthwhile to examine the many-body EC 
effects in the Xa model vis-a-vis Brueckner's electron gas 
theory and Sinanoglu' s many— elec tron theory. 
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The studies of Singal and Das [74] and Freeman [95], 
using the CC methodology provide an excellent insight into 
the many-body correlation effects in electron gas models. 

The importance of the ring diagrams in accounting for pair- 
correlation energy of an electron gas has been re-emphasized 
in their study. In view of the extensive studies on electron 
gas and their relevance in atomic problems we have under- 
taken this study of the X ct model with specific reference 
to EC effects. 

I V , 4 . 2 General Trends in EC Corrections 

The EC values reported in the previous chapter portray 
the general trends of diagonal diagrams-wise and pair-wise 
contributions for beryllium atom (vide Tables III,1 and 
I II. 2). Firstly, as depicted by the EC corrections in all 
the Approximations A to F, the inner-shells in beryllium 
atom are best described by the Xa local potential approxi- 
mation. Secondly, the electrons in the inner regions appear 
to be well localized as depicted by the small value of 
inter-pair EC corrections. However, the outer regions in 
the atom, which in the case of beryllium atom represent the 
valence 2s orbitals are poorly represented by the Xa local 
potential approximation and need to be corrected. This is in 
contrast to the results reported in the literature with 
reference to the HF method (vide Table IV.II). All the cal- 
culations that have gone beyond the HF method, including the 



Table IV. 11 Comparison of EC Energy Corrections 
in Be Atom 


Pair-wise 

Correction 

Contribution 

* 

to the 

Total 

Correction 

Total 

* Encrgy+T^^ 

Ref. 

1 2 

1 s 

1s 2s 

2s^ 


-4395.00 

r64B .00 




107 

-421 2.00 

-497.00 

-4488.00 

-9197.0000 


108 

-4208.30 

-497.00 

-4438.10 

-9143.4000 


109 

-4182.70 

-586.40 

-4535.10 

-9304.2000 


1 1 0 

-3758.00 

-457.00 

-4357.00 

-9572.0000 


35 

-4031 .00 

-547.00 

-3050.00 

-7628.0000 
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-9400.0000 

-14.66700 

1 1 2 

-4.41980 

-57.0836 

-9628.1753 

-9746.7623 

-14.663598 

Present 

Work 


Experimental non-relativistic Total energy = -14.668452 Hartree 

Reference [9^ 

* All values in units of lO”*^ Hartree 
+ 


Values in Hartree 
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CC calculations of Cizek [35] have indicated the importance 
of both Is and 2s pair correlation corrections. Comparison 
of our results with the exact experimental non-relativis tic 
energy reveals that with the inclusion of all the pair exci- 
tations an accuracy of 99.967% (for the a value of 0.769) is 
achieved. With the exclusion of the 1s and 1s2s pair 
corrections an accuracy of 99.962% is still obtained; this 
suggests that it will be sufficient to consider only 
valence-shells in the many-electron treatments based on the 
Xa model. 

The hierarchy of Approximations C,D and E all involve 
selective summations of various ladder diagrams. In Approxi- 
mations C and D the EC correction is over-estimated. As is 
well-known in the electron gas approximation to atoms, non- 
inclusion of higher order ring diagrams leads to over- 
estimation of correlation energy. The experimental corre- 
lation energies for oxygen and potassium are about 60% of 
that obtained from the perturbation result which does not 
include higher order ring diagrams [112], A similar result 
has been reported by Singal and Das [74], Inclusion of ring 
diagrams introduces the right trend in the EC corrections as 
is seen from the values obtained in Approximation A which 
includes rings and ladders to all orders. Also as expected, 
the second-order R5PT results give about 50% of the corre- 
ction obtained by the full calculation. It is worth noting 
that for the a =s 0.768, the higher order ring diagrams seem 


to be well-behaved. 
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IV, 4, 3 The Role of 'the Exchange Parameter Q 

The result of our calculations provide some insight 
into the role of the exchange parameter a . With increasing 
the EC corrections are generally seen to decrease in 
magnitude in all the Approximations except in those of E and 
F, The total energy with the EC corrections tends to decrease 
in magnitude with increase in a , although for higher orders 
of ^ the corrected total energy remain almost a constant. 
These results are again entirely in consonance with the 
electron gas results. The work of Herman et al [113] indicates 
that inhomogeneity causes the local one-electron potential to 
increase in magnitude and hence drive the a value above the 
GKS limit. Conversely, higher a values provide descriptions 
reflecting larger inhomogeneities which in turn decrease the 
total energy. Therefore, use of higher a values imply incor- 
poration of more correlation corrections and this causes a 
decrease in EC corrections. However, in Approximation E, only 
the particle-particle ladder summation is involved and this 
leads to incorrect energies and trends. The same is true of 

second-order RSPT corrections. Approximations C and D over- 
estimate the EC corrections, as expected. The inclusion of 
ring diagrams rectified this problem and like in the electron 
gas model correctly accounts for the screening of the charge 
at large distances. 
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The effect of a on the description of the EC 
potential merits a brief discussion. With the use of 
values of higher than 2/3, the EC potential in the 
model can be written as, 

''XC = "^^3 (IV. 3) 

As X increases, the EC potential becomes more negative and 
the repulsive interactions between electrons are decreased. 

In other words, the electrons tend to get correlated. This 
is reflected by the decrease in EC corrections with increase 
in Cl , However, there is no method for estimating the limit 
to which the negative potential can be carried to and 
therefore, any quantification of a with reference to non- 
uniformity is devoid of meaning. An interesting feature 
that emerges from our calculations is the ill-conditioning 
of the solutions for the cases where « values approach the 
GKS limit. This parallels the divergence behaviour of 
correlation energy of electron gas in perturbation theory. 

As in perturbation theory, the origin of this erratic 
behaviour rests with the most divergent terms arising from 
the ring diagrams. Exclusion of the ring diagrams stabili- 
ses the equations, but leads to over-estimation of EC 
corrections^ 

In all these calculations, the Latter correction 
potential has been used Tor obtaining the SCF Xct orbitals. 

The effect of the inclusion of the Latter correction potential 
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on the wavef unctions is to push the outer maxima slightly 
outward, as well as to increase the amplitude of the tails. 
This tends to inhomogeneise the charge distribution in the 
outer regions further, with the concomitant need to use 
higher a values [114], It can be surmised from the above 
discussion that a plays a very important role in deter- 
mining the structure of such local potential models. From 
the comparison of the EC corrected total energies with the 
exact non-relativistic experimental energies, it appears 
that the best EC corrected energy value is obtained for 
an a value near 0.760, This value of a has also been 
reported to yield the reference state energies which are 
closest to the HF energies [24], 

IV. 4. 4 



Tables IV. 6 to IV. 8 give the results of our calcula- 
tions of EC corrections in the Be isoelectronic series, 
and under all the Approximations outlined 

earlier. Five different o values have been used in the Xa 
potential. The trends displayed appear to be similar. 
Approximations A and B display an ill-conditioned behaviour 
for most cases. Approximations C,D,E and F, though better 
behaved, under— estimate the EC corrections. In all the 
later mentioned Approximations a consistent increase in the 
magnitude of EC corrections with increase in the atomic 
number, Z, is observed. Once again with increasing a a 
decrease in the magnitude of the EC corrections is observed. 
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As iTiBntioned earlier, an increase in Z along this series 
causes the contraction of the orbitals and homogeneises 
the charge density distribution to some extent. This tends 
to favour lower a values; however, the large density 
fluctuations encountered in the inner-most regions coupled 
with the incorrect description of regions farther away from 
the atom tend to undermine the local potential approxima- 
tion, The divergence of the ring diagrams in perturbation 
theory manifests itself through the ill-conditioned behaviour 
of the cluster-coefficients and energies as portrayed by the 
results of Approximations A and B in Tables IV.h to IV.S. 

The exclusion of these divergent terms arising from the 
ring diagram as in Approximations C,D,E and F assumes well- 
behavedness of the coefficients and energies. The EC corre- 
ctions obtained in the second-order RSPT approximation seem 
to display the correct trend. This fortuitious result is 
a consequence of the highly restricted nature of the summa- 
tion and bears no theoretical justification. The important 
point that emerges from these calculations is the total 
dominance of the valence-shell EC corrections. That the 
electrons in the inner-shells are well-correlated in the 
local potential model is borne out in all the cases 
discussed. 

I\/.4,5 Gradient Corrections to Exchange 

Herman et al [113,115] have suggested the use of a 
gradient correction for exchange within the X ct method. In 
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this method the exchange potential has the form 

t a + 3 G(p ) } (p ) ( IV.4) 

where 

G(P) = [i (^)2 . 2 -^] (1U.5) 

The potential is the X o potential for a = 1 and f'’ is 

an adjustable parameter. The introduction of gradient 

correction into the Xa potential was studied in the case 

of beryllium atom for three different values of B at 

a = 0.667 and a = 0.768. and the results are presented in 

Tables l\J ,9 and IV.IO. As can be seen from these tables, 

the gradient correction does seem to stabilise the equations 

for a = 0.667, although poor estimates of the EC correction 

are obtained. As pointed out by Herman et al [1 1 3 ] » 

lower B values are seen to be favoured for beryllium atom. 

Although the gradient corrections provide little improvement 

where EC corrections are concerned, they are found to be 

absolutely necessary to assure well-behavedness of the 

solutions for the case of a = 2/3. Use of higher a values 

appears to partially account for the atomic density gradients 

2 -4/3 

in these models. The presence of the (V p) p term in 

the gradient expansion (vane s as 1/r as r o ) w ould 
cause divergence of energies. We have taken into account 
this problem, through the introduction of a convergence 
factor [lIS] , so that the product of the convergence factor 

and G(P) remains equal to G(P) for a much larger range of r 
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than before. This is especially important at very low 
values of r, where there is a delicate balance between 
very large kinetic energies and potential energies and 
even a small change in the EC potential causes divergences. 
After incorporation of the convergence factor, eqn. (I\/.4) 
assumes the form, 

v^j- = a [1 + tanh { -I" G(P)} ] (I\/.6) 

where the hyperbolic tangent plays the role of a convergence 
factor, restricting the inhomogeneity term in the square 
brackets to the range + 1 . At very large values of r, 
tanh G(P)} approaches -1, such that v^^ approaches 0 

as r -»■ » . At very small values of r, tanh G(p)} 

approaches +1 such that v^|- approaches 2a v^^ as r 0 • 

At intermediate values of r, including the entire effective 
range of the Is orbitals, the deviations of { G( p ) } from 
0 are sufficiently small that tanh { ~ G(p)} = ^ f G(p)} 
to a high degree of approximation. For the cases involving 
higher a values in the Xci potential, gradient corrections 
do not give any significant EC corrections. The results of 
our calculations (Approximations A,B,C and D of Tables I\/.9 
and IV.ID) emphasize the importance of ring diagram terms 
for gradient-corrected potentials as well. The inherent 
shortcomings of the electron gas model are shared to some 
extent by the Xa method as revealed by the instabilities 
encountered in the cluster equations for a ->■ 2/3. However, 
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the exact treatment of the kinetic energy coupled with an 
appropriate scaling of the EC potential removes the diver- 
gences characteristic of bhe electron gas appr oxiniation to 
non-uniform charge distributions. 

IV. 5 CONCLUSIONS 

The import of our study rests on the analysis of the 
local potential Xa approximation. Owing to the use oF the 
uniform charge density approximation, the Xct model shares 
the problems inherent in the electron gas theories. The 
ill-behaved terms arising from the ring diagrams indicate 
the in appropriaten ess of the unmodified electron gas 
approximation for atomic structure studies. However, with a 
proper scaling of the local potential all the divergences 
inherent in the modified electron gas approximation can 
be eliminated. Our study prescribes the limits of validity 
of the local potential with reference to a variations 

and gradient corrections. In this context, it may be pointed 
out that, unlike in Brueckner’s study where only ring dia- 
grams were found to be important, or in Sinanoglu's study 
where only ladder diagrams were found to be important, in 
our study with the X“ model, it has been found necessary 
to take into account both the ladder and the ring diagram 
terms. Also, in contrast to Sinanoglu's studies with HF 
reference functions, where all the pair correlation corre- 
ctions were found to be important, our study reveals that the 
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only important pair corrections are the valence- shell EC 
corrections. This coupled with the simplicity achieved at 
the computational level, strongly recommends the use of the 
Xa local potential model for many-body atomic structure 
calculations . 



CHAPTER V 

EXPECTATION VALUES OF SOME ONE-ELECTRON OPERATORS 
WITH CLUSTER-EXPANDED WAVEFUNCTIONS 
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V,1 GENERAL 

This chapter reports the results of our preliminary 
investigation on the evaluation of one— electron operators 
for Be isoelectronic series obtained over the cluster- 
expanded wavef unctions ♦ Two of the hierarchy of approxi- 
mations to the full CC equations are examined in this 

2 

context. The expectation values of r,r and 1 /r are 
obtained over these wavef unctions . 

The calculation of the expectation values can be done 
in either of two ways. In the first method, the cluster- 
expanded wavef unctions are obtained first and the expecta- 
tion value of the operator is evaluated over these fun- 
ctions. Alternately, as in MBPT, it is possible to incor- 
porate into the Goldstone diagrams the operator directly, 
by replacing one interaction vertex by the corresponding 
operator in all orders [116,117]. This later method, when 
adapted to CC methodology has the drawback of invalidating 
the basis for the truncation of the expansion series after 
a few terms. In the limit of ultimate sophistication both 
methods are expected to yield identical results. In the 
following we use the simpler of the two methods, viz., 
evaluation of the property over the cluster expanded fun- 
ctions, with a view to gauge the effect of inclusion of 
higher order terms on the properties. 



V.2 THE EXPECTATION VALUES 


The CC expansion in the Approximations E and F outlined 
in the previous chapter is used to obtain the refined wave- 
functions. It may be recalled that in Approximation E, the 
particle-particle ladder diagram term is included resulting 
in a procedure equivalent to the selective ladder diagram 
summation to all orders in RBPT. Approximation F corres- 
ponds to the simple second-order R5PT expansion. Taking 

IT S 

into consideration that only the cluster coefficients, 

are significant in the expansion, the exact wavefunction may 

be written in terms of the cluster expansion formalism as 


]<!->= I $> + I t^® a* S io> (V.1) 

' ' Q a r a s r 

ci pr s 

where $ is the reference wavefunction and f the wavefunction 
obtained through the cluster expansion. The Cl expansion. 
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( V.2) 


may be compared with eqn. (V.l) in order to obtain the 
relationship between and t^®. As can be seen 
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For any operator 0, the expectation value expression assumes 


the form 
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( V.4) 


where, the c's are defined through eqn . (V.Z) and is the 
normalisation factor given by, 


= <'?['?> (V.S) 

For one-electron operators 0, the second and the third terms 
on the RHS of eqn. (\/.4) vanish owing to orthogonality of 
the orbitals; the first term in the equation corresponds to 
the ground state expectation value of the operator, i.e. 

< 6 = < $1q 1 5 > ( V,6) 

The last term in eqn. (V.4) can be written for the three 
cases, where the index sets (aPrs) and { a' p'r's' } may 
differ by none, one or more index labels. 

Case i) the index sets, {aprs} and {a'P'r's') are 
identical; the last term in eqn. (V.4) assumes the form, 

< -<ajdla>-<p|o)0> + 

a P ' a ' ° 

<r|lo|r> + <s|0js> . 


(V.7) 



Case ii) the index sets, 
differ by one label; 


{ aCrs} and { «' e*r’s'} 
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Let the differing labels be identified as i in the first 
set and j in the second set. The matrix element then has 
the form. 




Oi 5 


r' s' 
1 


> = <ilQ I j> 


(V.B) 


Case iii) the index sets differ in more than one label; 
the matrix elements vanish by virtue of the orbitals being 
orthogonal. 




ts 


a3 


0 


.r's' 


> ai D 


(\/.9) 


The expressions in eqns, (V.T) and (V.8) are used to 
evaluate the expectation values after taking into considera- 
tion appropriate numerical factors corresponding to spin 
summations. 


V.3 RESULTS AND DISCUSSION 

Tables V.l to V.B present the one-electron expectation 
values obtained in the two Approximations E and F for Be, B"*", 
and N"^^ systems. As can be seen from the tables, all 
the expectation values obtained over the cluster— expanded 
wa V ef unc ti on s are uniformly higher than the reference fun- 
ction expectation values. It is well known that the theorem 
of Miller and Plesset [33], which states that one-electron 
properties calculated from the HF wavef unctions are correct 
to "first order" can be extended to Xa wavef unctions as well. 
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Table V,3 <r > , <r> and <1/r> Values for B (Approximation E) 
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Whereas, two-electron properties are expected to be 
influenced considerably by EC corrections, most one- 
electron properties are not expected to show a marked 
deviation from the reference state expectation values. In 
view of this, it appears that the inclusion of ladder dia- 
grams to all orders (Approximation E) corrects the higher 
estimates obtained in second-order R5PT approximation 

(Approximation F) , As can be seen from the tables, with 

2 

increasing a values the < r > and <r > values decrease, 
while the <1/r > values increase steadily. At higher 
values of a a more contracted charge density description is 
provided by the Xa model and this results in the above 
mentioned trends. 

Due to the importance of these operator expectation 
values in the estimation of magnetic susceptibility and 
diamagnetic shielding, it will be of interest to examine 
the relative importance of contributions of various diagram 
terms in the evaluation of these expectation values. For 
this purpose, the direct method of incorporating the one- 
electron operator into the diagrammatic expansion is likely 
to be superior and further work in this direction is 

in progress. 



CONCLUSIONS 


The work presented in this thesis represents an 
attempt to describe many-electron interactions more accu- 
rately. It is different from other attempts, in the 
sense that it starts from a local potential-based model 
for studying these interactions and that it uses the CC 
method which originates from the powerful linked-clus tar 
theorem. This rigorously size-extensive, CC method, has 
provided an extremely useful framework for assessing the 
merits of local potential IPA's. The diagrammatic lan- 
guage employed in the method has Facilitated the under- 
standing and analysis of the terms encountered. Selective 
inclusion of diagrams has also enabled a comparison with 
various many-electron models. 

The CC method based analysis of the Xa local potentia 
model has revealed interesting facets of the underlying 
electron gas-like approximation. The limits of validity of 
the Xa approximation with reference to a variation and 
gradient correction have been outlined. It has been shown 
that by appropriate scaling of the uniform electron gas 
potential, through the parameter ^ , it is possible to 

obtain independent particle descriptions that serve as good 
starting points for many-electron treatments. It has been 
found that the best a-valued potential is the one that 
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mimics closely the HF potential; this is perhaps no mere 
coincidence, since the HF potential is the true potential 
of an independent particle approximated many-electron 
function. The superiority of the Xa model over the HF 
model in the computational realm is retained in the many- 
electron treatment as well. 

An intereating feature revealed by our study is the 
relative unimportance of the inner-shell EC corrections in 
beryllium atom. It appears from this that in many-electron 
treatments based on local potential models, consideration 
of valence-shell EC corrections alone would provide a good 
representation of many-electron interactions in atoms. It 
IS hoped that recognition of this fact will be particularly 
useful for large system calculations. The study also reveals 
the futility of employing gradient corrections to the poten- 
tial in order to account for non-uniform charge distribution 
in atoms. 

It is inferred from the CC analysis of EC effects in 
beryllium and its isoelectronic series that terms arising 
from both ring diagrams and ladder diagrams contribute to 
the EC correction. This is understandable since the X® model 
shares some features of both the electron gas and HF models. 

In the electron gas model, which serves as the basis for X« 
potential approximation, the ring diagrams alone contribute; 
and in the HF model, whose shell-structure like characteristic 
is well reproduced by Xa model, only ladder diagrams are 
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found to be important. In this connection it may be pointed 
out that for suitable values of ct , the Xo model overcomes 
energy divergences encountered through the application of 
electron gas description to non-unifoxtn charge distribu- 
tions; at the same time it has a computational simplicity 
not shared by the HF method, In this sense, many-ele ctron 
studies with reference Xa wavef unctions, can draw upon the 
merits of both electron gas and HF based many-elec txon 
theories. 

The application of CC methodology to go beyond the 
local potential IPA’s is not devoid of shortcomings. 

Firstly, the non-linear equations axe not well-behaved; 
this may well be due to the choice of the Xa model to 
provide the reference state, A detailed analysis of the 
terms that induce instability in the equations merits 
study. Also, it is not known, whethex a limited virtual 
orbital basis will suffice to describe the EC effects in 
larger systems correctly. Limited basis calculations on 
larger atoms, with "valence only" excitations would throw 
more light on the use of this approach for many-body cal- 
culations . 

In conclusion, it may be stated that our study has 
made an effort to analyse the nature of Xa approximation 
and gauge its use as a starting point in many-body calcu- 
lations. The results of our analysis clearly indicate the 
advantage of using the Xa local potential model in many- 
body atomic structure calculations. 



BIBLIOGRAPHY 


1 38 


1. H.A. Bethe and E, Salpeter, •’Quantum Mechanics of 
one-and two-electron Atoms” (Academic Press, New York, 
1957) . 

2. D.R. Hartree, Proc. Camb. Phil. Soc., 2A, 89 (1928). 

3. J.C. Slater, ’’Quantum Theory of Atomic Structure” 

(McGraw Hill, New York, I960), Vol. I. 

4. J.C, Slater, Phys. Rev., 509 (1930). 

5. P.A.M, Dirac, Proc. Roy, Soc. (Lend,), A 112 . 661 (1926) 

6. J.C. Slater, Phys. Rev., ^4, 1293 (1929). 

7. V. Fock, Z, Phys., 61, 126 (1930); 62, 795 (1930). 

8. J.C. Slater, Phys. Rev., 8J[, 385 (1951 ). 

9. D. Pines, ’’Elementary Excitations in Solids” 

(W.A, Benjamin Inc., New York, 1963), 

10. L.H. Thomas, Proc. Camb. Phil. Soc., 23,, 542 (1927). 

11. E. Fermi, Z. Phys., 48, 73 (1928); 4£, 550 (1928). 

12. P.C, Hohenberg and W, Kohn, Phys. Rev, 136 . B 864 (1964) 

13. N.H. March, ’’Self-Consistent Field in Atoms” (Pergamon 
Press, London, 1975). 

14. P. Rennert, Acta Phys., 3Ji, 219 (1974). 

15. W, Kohn and S.J, Sham, Phys. Rev., 1 40 . A 1133 (1965). 

16. P.A.M, Dirac, Proc. Camb. Phil. Soc., 2^, 376 (1930). 

17. P. Gombas, ’’Die Statistische Theorie des Atoms 

und ihre Anwendungen” ( Springer-Verlag , Vienna, 1949). 

18. P. Gaspar, Acta Phys,, 263 (1954), 

19. J.C, Slater, Adv. Quantum Chem. , S, 1 (1972). 

20. J.C, Slater, ’’The Self-Consistent Field for Molecules 
and Solids; Quantum Theory of Molecules and Solids, 
Volume 4” (McGraw Hill, New York, 1974). 



139 


21. K.H. Jonnson, Adv. Quantum Chem., X* (1973). 

22. R.P. Feynman, Phys. Rev., 5^, 340 (1934). 

23. J.C. Slater, J. Chem. Phys., 2389 (1972). 

24. K. Schwarz, Phys. Rev., B5., 2466 ( 1972). 

25. E. Clementi and C. Roetti, At. Data Nucl. Data 
Tables, 14, 177 (1974) . 

26. P-0. Lowdin, Adv. Chem. Phys., 2, 207 (1959). 

27. H.F. Schaefer, III, "The Electronic Structure of 
Atoms and Molecules" ( Addrson- Wesley , Massachusetts, 
1972) . 

28. A. Dalgarno, Adv. Phys,, H, 281 ( 1 962). 

29. H.J. Silverstonc and 0. Sinanoglu, J. Chem. Phys., 44, 
1899, 3608 (1966). 

30. E.A, Hylleraas, Z. Phys., 48., 469 ( 1928). 

31. W. Kutzelnigg, in "Methods of Electronic Structure 
Theory", edited by H.F. Schaefer, III (Plenum Press, 

New York, 1977), pp 129-188. 

32. E.A. Hylleraas, Z. Phys,, 347 (1929)« 

33. C. MjSlller and M.5. Plesset, Phys. Rev., ^ fit8 ( 1934). 

34. H.P. Kelly, Phys. Rev., 1 31 , 690 (1963). 

35. J. Cizek, J. Chem. Phys., 45» "^256 (1966); Adv. Chem. 
Phys., 11. 35 (1969). 

36. R. McWeeny, Rep. Prog. Phys., 4^, 68 (1980) and 
references therein # 

37. I. Shavitt, in "Methods of Electronic Structure 
Theory", edited by H.F. Schaefer, III (Plenum Press, 
New York, 1977), pp. 189-275. 

38. R. Courant and D. Hilbert, "Methods of Mathematical 
Physics", (Interscience, New York, 1953), Volume 1, 
pp . 175-6, 

39. P-0. LSwdin, Phys. Rev., £1, 1509 (1955). 

R.K. Nesbet, Phys. Rev., 102., 1632 (1958). 


40. 



140 


41* 0, Sinanoglu, Proc, Roy. Soc. (Lond.), A 260 . 

379 (196l)j J, Chem. Phys., 36, 706, 3198 (1962). 

42. H. Primas, in ’’Modarn Quantum Chemistry”, edited by 
0, Sinanoglu (Academic Press, New York, 1965), 

Vol. 2, pp 45-74. 

43. K.A. Brueckner, Phys, Rev., 2J.* 1353 (1955)} 1 00 . 

36 (1955), 

44. R.J. Yaris and J.I. Musher, J. Chem. Phys., 41 . 

1701 (1964). 

45. W. Kutzelnigg, Topics Curr, Chem., 31 (1973), 

46. R.K. Nesbet, Adv, Chem. Phys., 9., 321 (1965). 

47. J.C. Slater, Phys. Rev., 32, 349 (1928). 

48. R. Jastrow, Phys, Rev., £8,, 1479 (1955). 

49. L. Szasz, J. Chem. Phys., 1072 (1961 ); 

J. Math. Phys., 3, 1147 (1962). 

50. J.S. Sims and 5, A. Hagstrom, Phys. Rev., A ^,908(1971); 

J. Chem. Phys., 5^, 4699 (1971). 

51. D.C, Clary and N.C. Handy, Phys. Rev., A J_4,1607 ( 1976). 

52. J.D, Talman, Phys. Rev,, A 1805 (19BD), 

53. P~0, Lowdin, Int. J. Quantum Chem., 867 (1968). 

54. L. Brillouin, Actualites Sci . et Ind,, 159 (1934). 

55. J.S. Binkley and J.A. Pople, Int, J. Quantum Chem,, 
i, 229 (1975). 

56. J. Galdstone, Proc, Roy. Soc. (Lond.), A 239,267 (1957). 

57. B.H. Bxandow, Rev, Mod. Phys,, 771 (1967). 

58. J, Paldus and J. Cizek, Adv. Quantum Chem., £, 105(1975). 

59. I. Lindgren, Int, J. Quantum Chem., 512, 33 (1978), 

60. B.H. Brandow, Adv. Quantum Chem., 10,, 187 (1976). 

61. S. Raimes, "Many Electron Theory", (North Holland 
Publishing Company, Amsterdam, 1972). 



141 


62. H,P, Kelly, Adv. Theoret. Phys., 2, 75 (1968); 

Adv. Chem. Phys., 14., 129 ( 1969).“ 

63. E.S, Lhang, R,T. Pu and T,P. Das, Phys. Rev,, 

1 (1968). 

64. A. Hibbert, Rep. Prog. Phys., 1 21 7 ( 1975). 

65. F. Coester, Nucl, Phys., J, 421 (1958). 

66. H, Kummel, in "Lectures in the Many-Body Problem", 
edited by E.R, Caianiello (Academic Press, New York, 

1962), pp. 265-277. 

67. H, Kummel, K.H. Luhrman and J.G. Zabolitsky, Phys. 

Rep. , C 26, 1 ( 1978 ) , 

68. J.C, Slater, "The Calculation of Molecular Orbitals" 

( Wiley-Interscience, New York, 1979). 

69. 0, Gunnarson and R.O, Jones, Phys. Scripts, 21, 

394 (1900). 

70. E.P, Wigner, Phys. Rev,, 40, 1002 (1934). 

71. M. Gell-Mann and K.A. Brueckner, Phys. Rev,, 106, 

364 (1957), 

72. S.K. Ma and K.A. Brueckner, Phys. Rev., 165 . 16 (1968). 

73. K.S. Singwi, A. Sjolander, M.P. Tosi and R.H. Land, 

Phys. Rev., B i, 1044 (1970). 

74. C.M. Singal and T.P. Das, Phys. Rev,, B 6, 3675(1973); 

B 12, 795 (1975) . 

75. 0. Gunnarson and B.I. Lundqvist, Phys. Rev., B 13 . 

4274 (1976). 

76. H. Stoll, E. Golka and H. Preuss, Theoret. Chim, 

Acta, 52, 29 (1980) . 

77. C.F. von Weisacker, Z. Phys,, 26_» ^^31 (1965). 

78. M.S. Gopinathan, M.A. Whitehead and R. Bogdanovic, 

Phys. Rev,, A 14, 1 (1976), 

79. J.L. Gasquez and J. Keller, Phys, Rev., A 16, 1358(1977). 

80. J.A. Alonso and L.A, Girifalco, Phys. Rev., B U, 

3735 (1978). 



142 


81. R.N. Zare, J. Chem. Phys.. 45, 1966 (1966); 47, 

356 (1967) . ““ 

82. J. Paldus and J. Cizek, in '‘Energy, Structure and 
Reactivity", edited by W.D, Smith and B.M. Walter 
(John Wiley & Sons, New York, 1973), pp . 198-212. 

83. J. Hubbard, Proc. Roy. Soc. (Land.), A 240 . 539(1957). 

84. J. da Providencia, J. Nucl. Phys,, 401 (1963). 

85. T.P, Zivkovic”' and H.J. Monkhorst, J. Math. Phys., 1 9 . 
1007 (1978) . 

86. J, Cizek and J. Paldus, Int. J. Quantum Chera,,5, 

359 (1971). 

87. I. Hubac and P. Carsky, Topics Curr, Chem., 75 
99 (1978). 

88. G.H. Weiss and A. A. Maradudin, J. Math. Phys., 3., 771 
(1962) . 

V 

89. J. Paldus, J. Cizek and I. Shavitt, Phys. Rev., A 5, 

50 (1972). 

90. B.E. Adams, K, Jankowsky and J. Paldus, Chem. Phys, 
Lett., iX, 144 (1979) . 

91. H.A. Bethe and J, Goldstone, Proc. Roy. Soc, (Lond.), 

A 238, 551 (1957) . 

V 

92. J. Paldus, J. Cizek, M. Saute and A. Laforgue, 

Phys. Rev., A U, 805 ( 1978) . 

93. D. Mukherjee, R.K. Moitra and A. Mukhopadhyay , 

Molec. Phys., 3GI, 1861 ( 1975). 

94. F.E. Harris, in "Electrons in Finite and Infinite 
Structures", edited by P. Phariseau and L. Scheire 
(Plenum Press, New York, 1977), pp. 274-320. 

95. D.L. Freeman, Phys. Rev., B 5512 (1977), 

96. C.F. Bunge, Phys. Scripts, 21., 328 (1980). 

97. R.K. Nesbet, Phys. Rev., 1 55 . 51 ( 1967). 

98. C.F. Bunge, At. Data Nucl. Data Tables, 293 ( 1976). 

99. R. Latter, Phys. Rev., ii, 510 (1955). 



143 


10Q. F, Herman and S, Skillman, “Atomic Structure 

Calculations” (Prentice-Hall Inc., New Jersey, 1963). 

101. E.U. Condon and G.H. Shortley, ''The Theory of Atomic 
Spectra”, (Cambridge University Press, London, 1959). 

102. R.N, Zare, JILA, Rept. No. 80 (1966), 

103. B. Carnahan, "Applied Numerical Methods" (Wiley 
Interscience, New York, 1969). 

104. F.B. Hildebrand, "Introduction to Numerical Analysis" 
(McGraw Hill, New York, 1956). 

105. J.M. Ortega and W.C, Rheinboldt, "Iterative Solution 
of Non-linear Equations in Several Variables" ( Academic 
Press, New York, 1976). 

106. K.M. Brown, Comm. ACM, 10, 728 (1967); SIAM J. Numr. 

Anal., _6, 560 (1969). 

1Q7. B.Y. Tong and L.J. Sham, Phys. Rev,, 1 44 . 1 (1966). 

108. N.H. March and J.S. Plaskett, Proc. Roy. Soc. ( Lond . ) , 

A 125, 419 (1956) . 

109. D. Pines, Phys. Rev., H, 626 (1953). 

110. U. von Barth and L, Hedin, J, Phys. C, 2, 1629 (1972). 

111. N. Elyasher and D.D. Koelling, Phys. Rev., B 1 5 . 

3620 (1977). 

112. K.A. Brueckner, in "Atomic Physics", edited by 
B, Bederson, V.M, Cohen and F.M.J. Pichanick 
(Plenum Press, New York, 1969), pp. 111-130. 

113. F. Herman, J.P. van Dyke and I.B. Ortenberger, Phys. 

Rev. Lett., 22, 807 (1969). 

114. E. Kmetko, Phys. Rev,, A 1, 37 (1970). 

115. F. Herman and K. Schwarz, in "Computational Solid 
State Physics", edited by F. Herman, N.W. Dalton and 
T.R. Koehler (Plenum Press, New York, 1972), pp. 245-252. 

116. H.P. Kelly, Phys. Rev,, 182, 84 (1969). 

117. N.C. Dutta, T, Ishihara, C. Matsubara and T.P. Das, 

Phys. Rev. Lett., 22 , 8 (1969). 

110, H.J. Monkhorst, Int. J. Quantum Chem,, S H, 421 (1977). 



k good notation has a subtlety and suggestiveness 
which at times make it seem almost like a 
live teacher^ 

Bertrand Russell » 
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A1 THE LANGUAGE DF SECOND QUANTISATION 

The wavefunction for an N-Fermion system can be 
written in the farm of a Slater determinant (?•) of one- 
electron functions (<!)}. The determinant, on expansion in 
terms of cofactors yields 

«(1,2, N) = — I ^A^) 0,(2, N)(-1)'^'‘-^^ (A1.1) 

^ IeO ^ 


The subscripts refer to the orbital index, the number in the 
brackets to the electron coordinate index and n^i^ to the 
ordinal number of the orbital i in the determinant under 
consideration; the summation runs over all the orbitals in 
0 • Eqn, { A1 . 1 ) can be rewritten in terms of the annihila- 
tion operator, which is defined as, 

n .+1 

0(1,2, ..., N) = (-1) ^ 0^(2, ...,N) for <j)^eO(A1,2) 

and 5^ 0(1 ,2, . N) = 0 

Eqn. ( A1 , 1 ) then becomes. 


0(1,2, N) = T $,(1) a. 0(1', 2, N) (A1.4) 

/N X "" ^ 

where (1',2, N) is a function of electron coordinates 

2 to N and 1 ' is a dummy variable, A second expansion of 
yields, 


0(1,2, N) 


1 

/ifOCTl 


I , (1) «.(2)5,S^ »(1 

ij J J 


2', 3, ,,.,N) 
( A1 .5) 
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The adjoint of the operator is defined as the creation 
operator, which satisfies 


*A = 


"*aI 


( A1 .6) 


for all and It follows from the above that , 


B 


4. n . +1 

0(2,3, ..., N) = (-1) "■ 


0^(1, 2, ..,,N) for ^0(2, ...N) 

( A1 .7) 


0(2,3, N) =0 


for c 0 (2, . . .N) 


( A1 .8) 


The relationships between the a’s and a 's follow from their 
definition 


3 . a . Hh 3 3 . 

X J J X 


■‘t + 

a. a + a, a. 

i J 3 3. 


§ . a i + a' . a . 
13 3 1 


0 

0 


^ij 


( A1 .9) 
CA1 .10) 
{ A1 .11 ) 


where 6. is the Kronecker symbol, which implies 
13 

5= 1 for i = j 
=5 0 for i 3 

The Hamiltonian for an N-electron system may be written as 


H = Z + V 


= I z(vi) + I v(v,v) 
U y 


( A1 .12) 
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where z is the one-electron operator, v the two-electron 

operator and y , v are electron indices. For N-electron 

determinantal functions 0 and 

A B 

y 

= N < I z( 1 ) I ( A1 ,1 3) 

V I ^ ^ v(y,v) 1« > 

v<v “ 

= y2 N(N-1) <0^1 0(1,2) I 0g> (A1.14) 

Rewriting (A1.13) and (A1.14) in terms of eqn. (A1.4) and 
(A1.5) respectively, 


N z(1) |«B > = I <$.(1)1 J(1)| 4 (1)> «g> 


IJ 


= I M > ®i I (A1.15) 

^ 3 


N( N-1 ) 


2 <0 |v{1,2)hg>= y2 I <^.(1)<|>.(2)10(1,2)| <!>.(1)<).(2)> 

^ ijkl 1 J 1 J 




y2 I <ij!0!kl> <^1^1 (A1.16) 

ijkl A 1 j 1 kr B 


where any (j)|^ is represented symbolically by k. 
Eqn. (A1.12) can therefore be written as : 


H = £ <ilz|j> at a. + y2 I <ij|vjkl> at §t a. a. (A1.17) 

1 J ijkl i J X K 


xj 
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The second quantisation formalism represents a convenient 
way of handling the antisymmetry requirement and with its 
use the need to use a determinantal representation is 
obviated [ 1 ] . 

A2 FEYNMAN DIAGRAMS FOR PERTURBATION THEORY 

The second quantxsation formalism can be used to 
represent a general R5 expansion matrix element, 

<VRVR...RV> (A2,1) 

Jk. A 

where, V is the perturbation and R defined in terms of the 
projection operator is • 

A 

In general, if the perturbation V contains a one-and a 
two-body operator, it can be written in terms of creation- 
annihilation operators as, 

V s= V ^ ^ 2 

ik 

y2 5! <ijlv^|kl> at at a, a. (A2.2) 

ijkl 

A 

where 0^ and 0^ axe the one- and two-body components of V. 
The operator V can be represented diagrammatically in terms 
of the one- and two-body operators as shown in Fig. A2.1 . 


[l] R. Manne, 


Int. J. Quant. Chem., 


Symp. JJ., 1 75 (1977) 
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one-body diagram 

Fig. A2.1 Feynman-diagram 
operators 


two-body diagram 
for one- and two-body 




When the general matrix element eqn. (A2*1) is written in 

terms of eqn* (A2.2), a sum of 2^ matrix elements are 

□btaxned. Each of these is of the type < V R R -♦* 

^1 ^2 

R \l^> , V'^here any k. is 1 or 2. The expanded form of the 

n 

matrix element is very complicated; it is apparent that a 

A 

general term in this expansion contains products of V.j and 

A 


^2 matrix elements and creation and annihilation operators^ 
The operator part of the term has the form 


<aT a : 0 - a, R a a a a R R a a R 

ijlk mnqp ru 


* • # > 


Since the operator R rs diagonal xn the given basis, the 
operator indices occur equally in pairs, i.e,, corresponding 
to every creation operator aj^ there must exist an annihila- 
tion operator Only such of those summations of 2 terms 

where this criterion is met need be considered. The 
condition on occurrence of the operator indices in pairs 
can be represented by linking the two operators. For 
example, a single pair of operators obeys the contraction 


rule : 




1 50 

( A2.3) 


and two pairs of operators have the contraction relationship 



^ik ^ij^ 


( A2.4) 





- 6 . .) 

ij 


{ A2.5) 


where, the signs follow from the anticommutotion relationships 
between the operators. Di agramraatically this contraction 
implies the joining of the directed lines such that a 
creation operator line is joined to an annihilation operator 
line and further that the indices are matched^ The operator 
R = P/(E^-H^) generates the energy denominators of the RSPT, 


Typical second order RSPT diagrams and the corres- 
ponding matrix elements are given in Fig. A2.2 . Goldptone 
showed that when the 2^^ terms of a general RSPT matrix 

^ A A> A 

element <VR\/R,,.R\/> are expressed diagrammatically , 
all the Unlinked diagrams cancel mutually leaving non- 
vanishing contributions only from the totally linked dia- 
grams [2] . A diagram is considered linked if all the 
electron lines are connected to each ot,her either directly 
or indirectly via interaction lines or other electron lines. 
Further, all the diagrams representing the matrix elements 
will be closed due to absence of any uncontracted operators. 


[2] J. Goldstone, Proc. Roy. 5oc. (Lond.), A 239 . 267 (1957) 
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E<i|Tjl jXklT J l><aja^Haj5aj_> l><aj[ajE^^> 



* II 

I £ <illYgllcl><«il42lpq><ii^a3,^Sa2,4a^aj,> 



i 






Fig.A2.2 Typical second-order RSPT diagrams 
and matrix elements . 
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Diagrams I, III and V of Fig. A2.2 are examples of such 
closed-linked diagrams. Diagram II of Fig. A2.2 is an 
example of unlinked diagram where there are two disconnected 
pieces. Such diagrams do not contribute to the energy 
expansion. The general RSPT term may therefore be written 
diagrammatically . The rules for constructing the diagrams 
are as follows [ 3 ] ; 

1. For a given order of perturbation, all the topologi- 
cally different linked diagrams are drawn. 

2. These diagrams are labelled with orbital indices, one 

corresponding to an orbital in the ground state (hole 

lines) and the other corresponding to a virtual orbital 
(particle lines). 

3. The numerator of the term corresponding to the diagram 

is a product of a combination of one- and two-electron 

integrals. The denominator is given by a product of the 

terms ( \ t “ I ^r ) between each of the two 

hole lines particle 
lines 

vertices . 

4. For every pair of ’equivalent' lines in the diagram the 
term is multiplied by a factor of 72. An equivalent pair of 
lines is defined to be two lines going in the same direction, 
beginning at one vertex and ending at another, s.g. III Qf 
Fig. A2.2 has two such pairs. 


[ 3 ] B.H. Brandow, Rev. Mod. Phys., 771 (1967). 
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5. A sign factor of (-I)'*' is associated with the term 
which has in the corresponding diagram 'h* hole lines and 
* 1 ' closed loops e.g. Ill of Fig. A2.2 has two hole lines 
and two closed loops. 

6 . The term in each diagram is independently summed for 
all particle states and hole states. Exclusion principle 
violating terms [3] arising from these independent summa- 
tions are to be included. 

Application of the above rules to the fourth-order 
RSPT term is illustrated in the following section. 

A3 FOURTH-ORDER RSPT AND CANCELLATION OF UNLINKED TERMS 

Use of Feynman’s graphical techniques highlights the 
existence of the size-extensivity problem in limited basis 
fourth-order RSPT. Limiting, the basis to ground and double 
excited functions only, the fourth-order RSPT energy expre- 
ssions given by eqn. (1.46) can be expressed graphically. 
Expansion of the terms in eqn. (1.46) yields the terms 
represented by diagrams A-N of Fig. A3.1 [4]. The diagrams 
M and N arise from the term <V R <0 R V> R 0 > in eqn. 
(1,45). These diagrams called as bracketed diagrams are 
unlinked and by virtue of their being proportional to N , 
render an unphysical character to the perturbation expan- 
sion. As can be seen from the expressions in Table A3.1, 


[ 4 ] R.J. Bartlett, andG.D. Purvis, 
U, 561 (19T8). 


Int. J. Quant. Chem., 
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Fig.A3.1 Double excitation and bracketed 

diagrams of fourth-order perturbation 
theory . 



[Cable A3.1 


[terms corresponding to double-excitation 
and. bracketed diagrams of fourth-order 
perturbation theory 


Diagram Term 

A <<*Plbs><Mllu«><u«llxy><xy||«p>/D^p^3B^p^B^p^y 

B Ig-L <r8llap><apll6e><6Cllti5><nE;llrs>/i)^pj,gI)^gj,3D^5^3 

0 <rsllop><ap|l6e><rslluw><te||aw>/D^j,^I)^gj,^Dgg^^, 

D ^ 2 <rs| laP><»||rs><ap|li>e><6e| |uw>/D^p^^D^p„D.g^„ 

B Z <rsllaP><up|ls6><xy||ra><a6||xy>/I)„pj,gD_j^j,^^B^^^y 

B -I Z <rsllap><uwl|rs><ypl|w£,><a5||uy>/B„p,.,B„p^„B^^^ 

0 -I Z <rsl|ap><up||a6><aailn5><n5l|ru>/D^p^^D^j,3,^j])„5^^ 

H -I Z <rs||ap><aPlln6><u6lls5><nsliru>/D^p^^Dj,„j,gI)^j^^ 

1 Z <rsllap><up|ls6><wa|lue><ae||rw>/B^^^D^^j,^I)^g^„ 

K Z <rsl|ap><upl|B6><aw|leu><c6|lrw>/D^p^^I)^g^^J)gj^^ 

1 £ <rsllapXiiplls6><w^llre><ael 

aGuw 

M+N - ig Z <aPllrs><apl|rs><6el|M><6elluw>/B^gj,^Dgg^^D^g^g 

ajP»6j6»TijC represent the occupied orbitals» and 
rjSjUjWjXjy represent the virtual orbitals and 

<ij| |kl> = / dT^ f d^g <!jj^(l) <t'|(2) ■-“^= - *^1^^^ 

|ri-rgl 

^ioki = - e;L • 


there are no corresponding cancellation terms for M-I\l 
from diagrams and this introduces serious errors. 

However, inclusion of quadruple excitat.ian eliminates 
the unlinked bracketed terms, H and N of Fig. A3.1, in 
this order. The fourth-order terms originating from inclu- 
sion of quadruple excitations are represented by the dia- 
grams 0-W, in Fig. A3. 2 . The corresponding expressions are 
given in Table A3.2. The diagrams 0 to U in Fig. A3.2 are 
proportional to N and therefore well Dehaved, while V and 
W, the unlinked terms, are proportional to N [ 5 ]. However, 
these terms are identical to the unphysical terms arising 
from the diagrams M and N of Fig, A3.1 and therefore mutually 
cancell with them. This proves that even in fourth-order 
RSPT it is necessary to go to quadruple excitations to 
ensure the cancellation of unphysical terms. The higher 
order unlinked terms demand the inclusion of correspondingly 
higher n-tuple excited basis to ensure the well-behavedness 
of the theory. If only limited basis are used, tne unlinked 
terms escalate as the order increases and book-keeping be- 
comes impossible. This serious shortcoming of RSPT is 
easily overcome with the application of the linked-clus ter 
theorem to the expansion. The main premise of the linked- 
cluster theorem is that no unlinked terms figure in the 
summation; and the linked terms are well behaved with respect 
to l\l— dependence and cause no breakdown of the theory. With 


[ 5 ] K.A. Brueckner, Rhys. Rev., 100 . 36 (1955) 
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s 


u 



w 



Fig. A3. 2 Quadruple excitation diagrams 
of fourth-order perturbation 
theory. 



Table A3 ,2 Terms corresponding to quadruple- 

excitation diagrams of fourth-order 
perturbation theory 

Diagram Term 

^ 2 ^ l^’^^/^ccprs^epws^abru 

P+Q -ip <rsi |ap><uw| |&><PG| luwXaPl lrs>/D^p^gDpj.j;„B^^j,g 

E+s -|s <rs|lcip><uwl|6e><Se||rw><o;plluE>/D^pj,gEggj,^jjp„g 

T+n ^^2 <i;silap><™|l6B><rs||6e><apl|uw>/D_^pj,3Dggj,3D^p„ 

V+W <i'sl|ap>^<6el|uw>®/D2pj,3l)5g„ 


a,^ ,6»G jn, ? 

r , s,u,w,x,y 
<i3l|kl> = 


represent the 
represent the 

■f /<i'g 

Bi + Sj - 6], 


occupied orbitals and 
virtual orbitals and 


^(2) 
i. t) 





♦3.(2 
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the linked-clus ter expansion, therefore, neither unlinked 
bracketed diagrams like M and N of Fig. A3.1 nor the un- 
linked quadruple excited diagrams V and W of Fig. A3. 2 are 
ever considered, and the theory is rigorously size-extensive 
upto any order for any n-tuple excitation. 

A4 HOLE-PARTICLE FORMALISM 

The idea of a vacuum state with no electrons is tacitly 
implied in the preceding discussion. When dealing with 
ground states it is convenient to redefine the 'vacuum' to 
represent the occupied ground state instead of the 'true 
vacuum' consisting of no particles. This vacuum is referred 
to as the Fermi vacuum. Use of this representation results 
in a considerable decrease in the number of creation- 
annihilation operators that need be considered. The Fermi 
vacuum state is then 


1^0 > 



0 > 


( A4 . 1 ) 


where !o> is the true vacuum. 

In the Fermi vacuum representation, states occupied in 

I 0 > are referred to as hole states, while these un- 

I *o 

occupied are referred to as particle states. In this hole- 
particle formalism, the creation and annihilation operators 
and a. can be reexpressed in terms of the hole-particle 

i J 

* + 

operators b., b as. 



-t 


"f" 

(creates particles) 

A 

( destroy particle) 

b^ (destroys hole) 

b (creates hole) 

a 


1 60 


( A4 .2) 


a, g , etc. represent hole states while r,s, etc. represent 
particle states. i,j etc. are general indices and can be 
cither hole or particle states. The contraction relation- 
ships between a's and b's are as below ; 




a . 

1 



f 


b . 
1 


3 

1 



= 6 


ij 

= D 

= 6 . 

= 0 

= a 


ij 


= D 


if i = r, s, etc. 
if i = a »6 » etc . 

if 1 — ctf3 
if i = r , s, . . . 

if 1 = c(f g... 
if i = r, s, , . . 


i — 

A 

a 

1 

r i 

= 0 



F 

1 

A 

a 

3 

= 'iJ 

if 

i — f ^ f 



= 0 

if 

i = r, s, 

A 

a . 

1 

"It 

a . 
j 


if 

i = r , s , 



= 0 

if 

i = a , B , 



= 0 


if i = r,s, .... where 6^^ = <i|j> 

ifi=a»g» ••• 



0 


( A4.3) 



A5 NORMAL PRODUCT OPERATORS AND WICK'S THEOREM 


A normal product of operators (c, ,c , c^ ) 

^ ^ -^1 ^2 -^k 

where are either creation or annihilation operators is 

defined such that, all the creation operators are arranged 

to the left of all the anniliilation operators in the product. 


r” A 

N [c c . 

K2 




p it 


{ A5.1 ) 


where p is the parity of permutation 


kj k,-5 ••• k 

1 2 n 


(Xj X-^ ••• X 

y 1 2 


The contraction of two operators c^ and C 2 is then defined as 


'^l ^ *^1 *^2 “ ^ '^1 ^ 2 ] 

All contractions vanish except for the case ; 


( A5.2) 


®i ^’ij 


i Lt 

or b . b . = 6 • 

1 J ij 


{ A5.3) 


The normal product with contractions io defined as 


N [S, Cj ... ... ... ••• ... 


■■■ ^"' 2 '” iiL,”’ 


(-1)’^ c, c 


%S-”SV 


N [ C.J ^2 ' * * ^1^ -'1 


• # • 


( AS. 4) 
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where p is the total parity of permutation. 


Wick's theorem [6] is stated here without proof : 


« • « 






( A5.5) 


where the summation extends over normal products with all 
possible contrac kions of two operators, of four operators, 
etc. either second-quantised creabion, anni- 

hilation operators or hole-particle operators. 

A 4. A 

The Hamiltonian H = Z + V can be expressed in the 
normal product form, using Wick’s theorem. 


z= I <i|z|l<> 5+4,^ (A5.6) 

ik 

The operator product in this expression can be written 
in normal product form as : 


a 


t 

i 





k 


= N [a+ ] + 


"^^1 i,k ea , B . . , 


( A5.7) 


Therefore , 


Z = I <ilzlk> N [at a^] + I <ilzji> 

ik iea,3 


( AS .8) 


The second term in the above equation is nothing but 

< If I ^ ^ . Eqn. (A5.8) may therefore be written as 

*0 o 


z = 


< I ^ 




( A5.9) 


[6] G.C. Wick, Phys. Rev., BD, 268 ( 1950) 
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where , 


- 


I <ilz|l<> N [a a ] 

3.K 


0 = ^ <ijjOikl> at a'*' a a, 

i. 'J' I 1 j 1 I 

1 jkl 


{ A5.1 0) 

( A5.1 1 ) 


The operator product in this expression can be written in 
the normal product form as s 


al a T a, a, 
1 .1 1 t 


. A. <4* jt, 

a . a, 

J !■ 


N[a.^ 

a , 

J 

A 

"l ^ 

J + 

NElf 

, 

a 

J 

A 

h 

•f 




-.t 

a 

J 

“I 

^ A 

c] " 

N[a1; 

J 

1 

A 

"l 

^ -I 

"kl 

— I 

4- 



N[a+ 

F 

J 

) 

A A 

<] 

N[a-I 

r+ 

a 

J 

1 

A 

"l 

^k] 

+ 



nr 

M [a*2 

if 

a 

J 

-1 1 
A A 

ai aj 

<] 








N [a^ 

a 

J 

A A 

<] + 

N[a; 


< i 

.ik> 

i , k e a 


• * * 

<i|l 

> 


+ 

N[at 


<J 

!i> 





i ; 

, le a 

e 

Ju 



J ,l£a 

, 3 

< • * 

<j!k 

> 

J: 

, kea 

6 

[< i|k 

> <j 

ll 

> - 

<i|l> 

< J 

IkV] 









if 

,1c 

n f B 


( A5.1 2) 


Substituting (A5.12) in (A5.11), we get 


V = 


1/2 I <ij ! Oj kl> { N[a2 a? a^] + N [aj + 

ij kl 


- N[aj - N[a^ a^] ^ 6^^ “'il'jk} 

( A,B.1 3) 


J 

= + '^M + < *nl^i 



where , 


% = y2 I <ijlv|kl> [\![a^ a a ] (A5.14) 

ijkl 1 j X K 

A ^ A *1* A 

= I <i|g|t<> NCa^ (A5.15) 

ik 

with 

<i| g|k> = ^ 2 <iz\C\ k? > - < ie jO i E:k> ( A5.16) 

? 

the subscript A represents the antisymmetrisation of the 
matrix element. 

The Hamiltonian in normal product form is therefore, 

H=<$^!h1$^>+ I <iif! k> N[a!^ + 

X k 

yz I <ij!v|kl> N[at a^ a a. ] (A5.1T) 

ijkl 1 j 1 k-J 

where 

<ilf|k> = <ilzjk>+ <i|glk> (AS. 18) 

A6 FURTHER EXAMPLES OF DI AGR AMMATI CS 

The applicability of Feynman-Golostone diagrammatics 
extends well beyond perturbation theories. The diagrams and 
the rules for setting up matrix elements for non-perturbative 
theories are exactly analogous to those described earlier. 
The diagrams for the wavef unction , unlike the energy diagrams 
are open, and as discussed in Chapter II, a linked— cluster 
theorem can be proved for the wavefunction expansion as well. 
The elementary diagrams encountered in the non-perturbative 



method sre the and ”T” diagramsf these are presented 

in Fig. A6 . 1 . Linking of these diagrams invokes the same 
set of rules prescribed earlier. Examples of such linking 
are illustrated in Fig. A6.2. Since in both "H” and ”T'’ 
diagrams particle-hole operators occur in pairs, the linking 
results in none, one or two- open-paths. 
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APPENDIX B 

LISTINGS OF SOME PROGRAMS USED IN THIS 
THESIS 


I . 

II . 

III. 


Routine for 

Program for 

Program for 
Expectation 


General Matrix Element Evaluation 

the Solution of Linear CC Equations 

2 

the Evaluation of r,r and 1/r 
Values 
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n 


V 

c 


c 

r 

r 


c 


i 


CALCULATE THE BENERAL XNrEGRA!. H ( JAr JB? JCy JIO 

^ \ ^ V rn “V /n ^ # V /*\ /^ ri* i>n t ^ /*v ^ J*' ^ ^ ' 


Thist sub rout tri<'» evaulates the ^^enersi matriK 
eJemeiit over the W3v©fijnnt:ions froiri the XaJwha f-roiaraiTi* 
The «>uhroiiti)»e has^ jii coramori the rahiaJ wsverurirtlons of 
the marri Xalrhri f-roiiratru Part of the >*ro«sraiTi if> adapted 
frowi the roiitiries of 7are<vid© “ JILft Report SOif 19A6) 
The 3o iswmbola rioceaisare are s?ierieroted bw a ceparate 
prodram arid input rnto th'i<i prodra«i4 The routine ('alio 
the iiited routine iri which the antodratjon is done 

/'V ^ ^ /♦v /K PI / V /*% /*v /•v in. / V x'\ ^ ^ m ^ le* /*s ^ eii /n •-* e 


# i PI 


The 3 j coef fictents are c'orivorted irito the reauired CG 
matrijc for multiplication of tfie radial intedraJs* 

RRAIi *f (((CCi<TrJ»K)yK=l» 3 )»v}~ 1 .» 4 )»T-^lf 4 > 

HO 22 I=t»4 
DO 22 J==.l»3 

22 KA<IfJ)-=040 

KA<1»1)=1 
KA<1»2>“2 
KA( 1 i- 3>"=6 
KA ( 2 !» 1 ) ®3 
KA<2f2)=«7 
KA<3i-J >“4 
KA<3y2)-=9 
KA<4yl)«ri 
KA<4y2)®9 
00 198 JK®ly3 
DO 199 I®1»4 
DO 199 J®ly3 
DO 199 M®1»4 
DO 199 N®t»3 
JA®KA<I» J> 

JP®KA<MfN> 

IF (JA.EQ^O) GO TO 199 
TF <JB.EQ40) 00 TO 199 

CONCJAr JB» JK>®CG<:iyMy JK) 

19«-> CONTINUE 

190 CONTINUFi: 

DO 110 J®2y441 
SNN<tyJ)®SNL(ly J) 

SNN<2f J)«SNl (2yJ> 

DO ill K®3yS 
I I I SNN(Ky^^)“'SNL<3y J) 

SNN ( A y J > ®SNL ( 4 » J ) 

DO 112 K®7y9 
TIP SNN(Kf J)'-=SNL <5y J> 

1J0 CONTINUE 

ML<1>““0 
ML(2)®0 
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Ml <3>-l 
ML (4)=-0 
Ml. <5 1 

Ml-(<&)- 0 
Ml.<7)“"3 
ML<B)=0 
Ml <9)-=-“l 
11(3 114 

DO 114 JH-1»9 
DO 114 JC=-ti-y 
DO 1J4 

MEt: ( .lA y JD i. JC y JD ) ~0 . 0 

1 F<<MI (JA) fMI.(JB)~Mt (vJO) -Ml. (jrO)*NE»0> GO 30 tl4 
DO 115 JK~lyi 
JL -JK-J 

CAIL INTEGCJAyJPyJCyJIiyJI yREStiLrylf) 

■I :>9 COO •CGN < JA y JC y ,)K ) sKCON ( JD y JD y JK > 

115 vri: < JA y J''WB y JC y JD ) =M£- E < JA y JB y JC y JD ) 3 KE01H.T*C00 

!14 CONTINUE 

WR i TE < 2? y # ) ( ( < < UEE < J A y JB y JC y JD ) y JP' J y 9 ) y JC“ 1 y 9 ) y JB « .1 y 9 3 y J A- 1 y 9 ) 

Rf 3 URN 

END ^ 

<' 

rh 4 .i:> feubroutirm uistriK e|lf?m<!jrits usinii thf' 

e HK^thod dJven in ''Thp Ih*»or'J of Alomac Sr-ectra' bn 

x< E*UtCoridori and 0*H»Short'»l'.j* 


SUBROUIMNE INTEG i JA y JB y JC y JD y JK y F\*ESUl f y H ) 
hilcjd oubroutine uomr'Uloa Rf.*K 21 rite«iralat 


C 


c 


70 


^«FSUI..3 

COMMON 

COMMON 

COMMON 

COMMON 


R < JAy JBy JCy JI. 1 ) 
XNUM y RSOAL E 

V y SNLO 

XJ y SNl 

RU2 y RUX 


RUEXCH y XI 
NKKK y EE 
, RIIINL? y RUFNL2 


y RSATOM 

y R y RU y 
y NNL7 y WWNL y 
y RUINl.3 y RUF'NLl 

Cf)MM0N/K0NS1 /MEOH y NCSROS y Z y NBtOCK ^ ^ ^ , o ,• 44 1 ^ 

DIMENGION X(441)yRSC0RE<441 ) yRU2<441 ) yRU<441 ) yNNLZ<5) yKl441) y 
IRSVAl E< 441 ) yRU3<44t ) yXX <441 ) yWWNL(5) y V<441 ) yRBAT0M<441 ) * 

XJ<441)yNKKK<5)yEE(5)yRUEXCH<441)yDENM(441)y 
*< BNI (* 3 * 441 ) ySNl 0<441J yA(5y5) y0l.NS<441) yXNUM<441> 

DIMENSION RUINIJ (441 ) yRUFNl 1 < 441 ) yRUINL2<441 ) yRUFNl 2(441 ) 

EOU I VALENCE ( RSCORE y XNUM > y < RS VALE y DENM ) y ( X y R ) 

DIMENSION SNN(9y441) 

COMMON/RAD/SNN 


XI<l)'-0*0 
DO 70 I~2yMESH 
XJ(I)“1*0 
CONTINUE 



/3 


r, 

74 


V 


10 

ir. 

?o 


30 






45 

50 


55 

60 


IF<JI\) 20A0i»/4y/l 
mo 73 J - 1 y JK 


DO 7? I«2¥Mi nil 
XX<I)- Xl( I )*ua) 

CONriNOt" 

CONT INUE 

xi(r> NOW f^ONrAiNO nif- Ar<RrtY i<< i )!{t!«iJK 
CONrXNUF 
XJ<l)-0»0 
BN} 0 U ) -0 ♦ 0 
OLNB( I )-0»0 
DO 75 I -r»MI r>H 
HNI.OC 1 )“BNN^JA!» I ):{iHNN< JOm 'I. » 

1)1 NB CD- SNN ' vID !» r ) 'KBNN ( JD / I' ) 

V>t T or ARNAYO Ul" RF PI- Ari-JH Y UBI-Ji WAVE" 
XJt l>' J ./<XI ( 1 >Mk I )> 

XJU) NOW OONIAINO ARl^AY 1 /!?( .E )*!!{JK-} 1 

CONTE nue: 

RBror«:r c J )' 0*o 

DO 15 ME. BH 
RSCURI.d ) - 01 NB(1) 

CONTINUF 

}E(,)K1 r<n30*?o 

DO ?5 J-2s*MlSI-l 

KBCOKh ( £) 'KGCOKI- < J 1 1 J ) 


coNnwir 

CONI I NOE 

CALI BME‘SN < HlHCOICE' y ECHOAl l;‘ y 1 1 y NDl OOK y 40 ) 
E<BVALE; ( I .)--=0»0 
DO 35 1--2yMLf-)E} 

RBVALr < I ) -OX-SOAU: < 1 ) *BNI, 0 ( T ) >KX J < J. > 

CON rt NUE 

CAL 1. S T NT < F<S VAl E:. y II y RE: SUI I y NDLOCK ) 
E<SCORL (1) ■•0»0 


DO 40 T-'JiyMLSI-l 
K«CnRE<I)«-OLNB( l))ftXJ( I ) 

CONT ) NUE 

CALL E>JNr (E^BCOE^E.■yllyTOT■EN1 yNKl OCK> 
CALL BMF‘BN(R‘SnOPI yEiiSOAl L i.|-lyNDL0CKy40) 
EilOMAI. Er <1>”-=0,0 
DO 45 i-n^yME-SEI 

KSMALE ( T ) =-SNLO ( T ) >!c < TO E t NT -RBVAl E ( 1 U 

CONE I. NUE 

3:r(v)K) 50 y AC) » 50 

DO 5Ei J-2yMFSE-l 

FisvAiE (j:)-RSMAir.c.T.>*x,r<.n 


CON I T NUE 

CONTINUE o MKinrio 

CALI STNT C RBOALE y H y K*F-BUL2 y NPLOCK ) 


F UNCI ION 


1 71 


ru’ODUCis 
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Rl” SUL F"- »1 1 0^ ( I'll Slil I fM'..SUI >*) 

1000 Rt'TURN 

2000 «RnF<?i»20on jk 

200 ! I-ORMA I < ' fROGRAM lINTFG-HO Nl'G* ,) 'yJ3) 

S roR 
TNU 

SUP'ROUT J NF SMPGN t Y s- Y ] y liX » NP » N ) 
rUMLNSION Y(441>v Y!<44J) 

ASUM-“0 . 

PSUM--0 , 

Yi ( n.= 0» 

H-PX/;?, 

I=-J 

G= U/4* 

kmax-=ni;n/2. 

( NfN MUSI BE AM I, MEN INILOER 

BO GOO J=1,NP 
HO 400 K-L'-KMAX 
J»I+2 

BSUM-BSUM+H)K(Y<1- 2) 1 4*jf:Y( 7 -I >4Y( I ) ) 

YT< t> -ASUM-iBSUM 

400 Y i < 1 ■ J ) =YI< J‘- 2 ) i G)j! ^ 15 ♦ )KY ( T~? ) { 0 ♦ JKY < t J ) " Y ( T ) > 

ASUM==Y L < 1 ) 

BSUM="0* 

G“ G4G 

000 CONriNUE 

re: turn 
k:Nri 

SUPF<OUIINE SrNTadbVAl UE!i.N1,N) 

G SJNJ~EXTE:NriEB BlME"*BONS RULf rNIEGE^ATlON FOF<MULA 

n OBTAINS INTEGRAt E'ROM 0 10 1NI-INITY OF X ANB CALI.S IT MALUC 

ra MENS ION X<44l) 

CURSUM™0.0 

F?UNSUM--=0*0 

HMESH-H 

I---1 

C 

DO 10 J-^Ii-NlN 

C>fc5!i NIN 18 THE NUMBEIR Ol-’ FORIY E’OINT BLOCKS OOER WHICH THE 3NTLGRATE0N 
r L8 BEING CARRIEB 

cuE?8UM---x ( :r > 
r-iii 

no E> K^l*19 

CURSUM=CUR8UM+ 4 ♦ 4cX < I ) 4 7 ♦ >!cX < 3 + 1 ) 

I.-I.+2 

U rONIINUE 

0 

CURSUM-CUkSUM i 4 * 0#X < T ) I X a + 1 > 
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RUN'JOM-. kUNHUM+HM^ SHJicCUJkSjim 

ni)k3UM-“0*o 

r-i-i+f 

10 CDNUNUk 

VALUF-’-kllNOUM/.i^O 

Kl- TURN 

riNIi 



M ^ ‘y ^ J \ *4 t / % 


1 74 


c “• - Thia f^ro^^ram aoJvc's the ■lirK**^?' i-att of the-' 

0 couf'led oltjvt<'r o^tuot lonu* Kho I titosr eonefti otiu 

<•' oriaind from orifc- oi^en- >• ath rj33^r»i»a oro ftial 

o in tht; aorial ordfe^r I bo M(N"'M) arid then 

c 

c aro ii‘-»Led two oeen-path diai^ram oeuatjona 

a jn the oonijriHsnd aewu(^rit ufl oiMCN-M)'!-! trj 

(' <N~M;)KN“M>5j:M5>:M« Tht? Jotlr'r art- thon redoced iiMiid the 

c aymmeiry t(AI.!'BT iihiyS)" b(BT yAl. s'BfIO ♦ 

I.' Tfhe irivorae of the eoef f i c i ent unntt i;; ia raJ.t. uiated arid 

r multiplied h*.i th<-* coeffjcierit vector* 


* \ / V /•i /n ifn c\ /*> /*v /*\ e j % /•» 




^ % / li ^ 


COMMON/L ABEL/IND <?92?9v9) 

INTEGER AL.Bl i^R^Gi^Ur WunLi-EI-'yZI 
COMMON N»M»NMi.M1 

BlMtNGTON XX<2I0) i-TT(B!»2r9?9)!-V<9y?i.9ii9)»EN(9) 
niMf.NOION X<l:l9)iiML<9) 

BIMLNSION yiL(9i-9) 
niMENClON AX(li9) 
miMENSTON A<119? I 19) i.B(tl9) 

ORENCUNi r“20iim:oicr-^»sK' »rii E-'LJN.UAr' > 

0|:*E;N < UN rT==21 ? be V1CK« ' BSK ' y FI LE-- ' MET * BAT ' > 
0fl“N<UNIT“22*BEVi:Cl>^liSK'' yFILE-'MEl *OUT' ) 

N - 9 
M 2 
M I « M-t t 
NM~N-M 
N0LB=-21 0 
NEG- n9 
BO 10 l«lyNE« 

EU'I )“0* 

BO 10 J^JyNEO 

A< lyvJ)-<^0. , „ 

c - • • Inr-'Ut data for the iHiear eouatioria " the vyvL and Fri 

Q ...™ Matrix elementa calculated from the Xalcha prodram* 

JO CONTINUE 

WRin: <22y999) 

999 FORMATdOXy '/// COUPLEB CLUSTER EQUATION RESULIo /// ) 
RrAirK20?9000) 

9000 FORMAT <lHly 

1 60H 

2 > 

RLA0<20y3K) ALPHA y BE! A 

fiEAIi(20y3t!> <<(< 0 (TyJyKyL)?l. -l»N)yK™lyN)yJ=l»N)yI=“I.yN) 
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Rh:AD(20y3tt>( <VL<IyJ)yJ--:| »N)y l=-J yN> 

RLAri(20y*)r.0 

REAlK^Oy^f) (I•:M< f ) y 1. -rj ,N) 

WRIKE 1 22 y 9000) 

WRirk"<22»99) Al. RHAyKfcTA 

99 KDRMAI <5Xy ' AI.PHA yplS* 9 » 1 OX y " REI A yl J.b.9//) 

WRT.TI:; (22^9001) 

9001 rORMAK" THE: OlvBT fAI.. I:NE-R0.1F..S I'N THE SCQUE NCr 1 55»2S«2F-<,y;r:Jy,5H:»' 
WRirEi;(22y>K)ON<X>yJ.:^,|y9) 

» ' « . 

p I nd&'jj T.ri‘ji Lhp 1^2 (’Of?fr i c i, c*nt<» m r-'rorHo <n'd<“'r 

K=- 1 fj 

jiin Ifi AL^- lyH 
HO BT-lyM 
HO 15 K-H1yN 
HO J5 B=''M1»N 

jr < (Al.4E0«BT»ANH,Ka.l' ♦«) ♦OR«AL J.I ♦HI) 00 70 IT 
GO TO 14 

iS TNlKALyBI yRyG) -.K 

K=-K+l 
GO TO IS 

:I4 CND<Af..yBTyRyrf)--0 

15 CONTINUE 

<-'••• . . .. i. ... . 

C" “ oricj oP6>n--f'atft eciurfii otiiSf with tl jnd i2 co6'('('irx©i‘tl% 

„ 

DO 20 AL»-lyM 
HO 20 R=MlyN 
IsHl 

HO 21 7I“ly2 

B< t)-«B<I)+VI..<RyALi- 2 » )KV(R» 7 J' y Al.yZ.1 >iM<Ry 71 yZIyAL) 

21 CONI TNUE 

A<.E *l)"A<lyl)fEN<R) -f N(AI.) 

HO 22 S™MlyN 
HO 22 71=1 yM 
If- <S4E04R> go 10 22 
.5 TNDEX2(ALyG) 

A(|yJ)=A(XyJ)T24)l{0(Ry7ry0y71) -V(Ry7J yZlyi'5)"ML(RyG) 

22 CONTINUE 

HO 23 BT=lyM 
HO 23 71=1 yN 
irtHT^EOtAL) GO 10 23 
J-nNHEX2<BTyR) 

A<.l yJ)=A(IyJ)“ ( 2^ #0< HI y 71 y Al y 71 )-0< HTy 71 y 71 y Al. )“0L <B1 yAL. > ^ 

23 C0N1INUE 

HO 30 BT=lyM 
HO 30 S"“"MJyN 
J-INHEX2<B1 yS) 



Ad)- J)™A(l¥>J>+2»!{!V<h'?lt<‘l( i-AI.i-S^“M(I«!-BTi-,‘NAL ) 

CONI INUE 

r«0 40 bT-UM 

DO 40 S«M1 t-N 

DO 40 U«MJ »N 

J“INDrX4<DTvAL!>S?in 

Ad» J)“Adi-J)+:^’t:0(R?Dr-0»K)“V<l?i-DT»0»U) 

CONTINUE 

DO SO DT--1?M 

no SO DI.-IkM 

DO SO S-MIdN 

j=3 NriEX4<m ^ni s-s^fo 

Adi- J) =Adi- rDTfAl .H3- V(D!.i-tn »G»Al )) 

CONTINUE; 

DO 52 

DO S2 U«M1»9 
DO S3 

J = |NDEX4<ALrDLi.r<yU) 

Ady J)«AdyJ>+2»J|{<2«Wa;i|, yZI y UyiT I )-M<riL y ?/ y ZI »U) ■ 01 (DLyU) ) 

CONIINUE 

CONI INUE 

DO S4 DI.”lyM 

DO 54 U=»MlyN 

DO S4 ZI~lyM 

J~ tN»C:X4<AL yDl.yUyR) 

AdyJ)»-AdyJ) -( 2 * 3|{V( DL y Zl y U y Z1 )"0( DU y Z J y ZJ yU ) -OL <DL y U ) ) 

CONTINUE 

CONTINUE 


twu ope‘rt"-»»< 3 th d 3 f}«Jre}m 6‘«ii3tioris> with t1 and 12 toef f i cj ©rii*i» 


I=-1.4 

DO 60 AL=lyM 
DO 60 BTs-lyM 
DO 60 R»MiyN 
DO 60 S-“MlyN 
k ” j.ride5<4 <al ybiv ryi'>) 

|,f < <AU*EQ«BT*AND.R.I.E*CK0R»AI. ♦LTyDI) 00 TO 71 ? 

00 TO 60 
CONTINUE 
!• 1+1 

Dd )--0<RySyALyBT) 

Ad y I)“AdyI>+EN<K'ULN<0) "EN< AL>“EN(DI ) 


DO 61 U~MlyN 
DO 61 ZI=JyM 


J”-lNDEX4<Al..yBT«UyS3 

r< I ^ J ) ^- A ( i ) + < y ♦ ^ U ^ ^ ^ ^ ^ ^ ^ 

CONTINUE 
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6'J 




64 


70 


00 


90 


HO 62 U“H1. 

no 62 Zl-lcM 

J^JNLiEX-KBTyALyUi-K) 

IF (U^FIQ.Si GO 10 62 

A<1'» <2»!{(0<fb7.t »U)-7I )--M<Sp7r»/.1i.U>"yi.<SyU>) 

CONTTNUE' 

no 63 Dl~.UM 


no 63 

J=INnEX4(liL »B1 »Ri»S) 


IF (ni.*EO*Al.) 00 TO 63 

A<Iy J)"'A<Ti*J)“<2**V<Al y/J.?ri| )-7’! ) “•M< Al.^ ?t s- 7 1. s-IU, ) -VL(Al i-ni. ) ) 
CONTINUE 
no 64 nL~-J»M 
no 64 

J~lNnEX4tnLi-AL»S?fO 
IE (ni..*EQtBT) no 10 64 

A<I» J)“A<Ii»J>‘ (2*3|{U(B1 y/li-ni, »/r) "VOFn 1-7.! i-^IyUL) • 01 - fB1 i-ni.) ) 

CONTINUE 

no 70 U-MlfN 

no 70 BL-Ii-M 

J=.lNnEX4<nL»BT»UyS) 

A<I»J)“A<T»J) ♦2*«V<H?nL»AL^U)- V<Fi*-ni ^Ui-AL) 

J==INnEX4<nL»AI •Ui.R) 

Adr J)-«A<TvJ)+2«)KV<Si-ni »BT^U> -V(S?nEylhrn > 

J™INnEX4(nL»BTrS»U> 

A<I»J)==-A<.TrJ)“0(R»nLyAI.?U) 

J-« I NnEX4 < BE f BT f R y U ) 

A ( I y J ) ™ A < I y J ) "• U ( S y ni y U y AL > 

J-TNnEX4<nLyALyRyU) 

A(IyJ)«A<IyJ)-V<Syni yBIyU) 
j!.= INnEX4(nLyAEySyU^ 

A<IyJ)-'='A(JyJ)~0<f^yBL oUyBD 

CONTINUE 

no 80 U«MlyN 

no 80 W-MlyN 

J~INnEX4(ALyBTyUyW) 

A<IyJ)=A(IyJ)40.S*V<Ry0yUyW) 

J~INnEX4<BTyAl.yUyUI> 

A ( I y J ) “ A < .r y J ) 4 0 * fj 0 ( 0 y r< y U y W ) 


CONTINUE 
no 90 BL^lyH 
no 90 EE*'“lyM 

I NnEX4 ( PL y EP y R y S ) 

A<Iy J)-"A<IyJH0>S5!{V<ALyBTyDl..yEr^ 

J-=lNnEX4<nLyEPyHyR) 

A ( I y J > --A < I y J > +0 ♦ HJftM < B1 y AL y Bl y Ef ) 

CONTINUE 

no 91 U-=MlyN 

J'~INnLX2<Al yU> 
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Ads' J)-A<'l's>,J) ) 

J^-TNIiEX2(Br i-in 

Ady J) "•Ad»,jMV<S»ri:?UcAl.> 

<ri CONTiCNUF 

HO 92 DL- UM 

J«.:TN»EX?.(I.iLyF'1 

AdK J)--Adi- J)”V<liL > 

J-=-lNriHX2(liL ^8) 

Ad» J)"Ad !- J)- V<r>l i-K»BTf »AI } 
72 CONdNUt 

60 CONdNOr 


i9 
c- •• 
c— 


tt50 

C" " 
C“ “ 


r-— 


caJ I thf» ro<.itiiie f'or ih<? uolulion oT the* ir>'j'r>t('iTi oP e<!Mtsi a oi'i<» 
CALL LXr_0<AyNLai-r{»X> 

che('K iP Ihe tho a«itatiOii lu / ot'onracy 

DO 19 l-=.|di9 

AXd)-'0* 

HO 27 J"~1».U.9 
AXd1“AX< I )^A<Xf J)!KX<J> 

»irF-==AXd )•••»< A) 

TF(DlFF»Lr« l*r“7) GO TO 19 
WRITC <22y)f:)D.i;FFH 

CONTIhJUF. . , 

rpc'onv^* ry i oi'i of' Lho roctucod niiitiHtjr "to tho ori^ilri'M iiiiivibaT 
coc'ffj 010015 4 tor the case Be serjea the ronversion laFos 
tile 119 coefficients to 210 coePf 'iCjeriLs4 
Iwl4 

BO ISO ALr-.li^M 
BO 150 

DO 150 R~Mt»N 
DO ISO S»M1»N 


A = T41 

TX~INDEX4<ALi»DI s-R^G) 


XXd) =XdI) 

CONTINUE 

call the roMtine that cvaloates the oichanse-correlad ot. 
corrections usir.i1 the cncf f i c i ents and ms^trix elements* 
CALL COREN<XXi*Mi»ENyNyrOt-M) 

output from this pt'osram to be used as .input in the 

nonlinear program* , • a\ r-i oi 

WRITE <21 i-JK) < < (Odi- J f9) ^K=--1 ?9) J“-l j9 > » X ~1 1-9 ) 

WRITE (21 y#) <<0l (lyJ)?.!"! y9)? t-“'ly9) 

WRITE<21y*) <EN<T) yl-l r9> 

WRITE <2ty*) (XX< T ) ? 1 - 1 y 2.1 0) 


STOP 

ENn 

orderina the tl (‘oeff icients scauentiall'J 
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rUHCTTON INrCXStALi-H*^ 

COMMON N»M»NM»M1 

in»e:x:? « ) + <k M) 

NETURN 
f ND 

ordf'rins the t2 cueff j.<.' ients '5G«u(*iiij 3.I ly 
FUNCTION ■lNDtX4<ALyBTKR?S) 

COMMON Ns-M^NM^MI 
COMMON/LA0F L / TNIK ? » IN 9 y V ) 

TNIFOE-R ALfUTi-R^S 

IF <<AL*E(NlBr*ANIi,Rar.*{J).ORtAt ♦l.'KUK) 00 10 1,1 
I NliFX4= T NP ( HT y AL y IN R > 

RETURN 

T NPEX4- 1 NP ( AL « PT y R y S ) 

R6 TURN 
END 

routine for e>otvlii<i th€‘ liiiowr i:>y«iLein of ooucJtiarib* 
thiy routine is b^sed on the LU deromnosAiiori eJ «fon thiiti» 
SUBROUTINE I.TE0( AyNNyPyX) 

M 3 < 3 ter rro^roni to roordinate the uubroutineii} J 
E<ecompySolvey .tmeruv eiid SifiM* 

DIMENSION A<NNyNN) yB<NN) yXCNN) 

COMMON/ AREAl/UL < U9y119) 

CALL nECOMP<NNyA) 

CALL SOLUE(NNy]HyX) 

CALI. I MPRUV < NN y A y B y X y II IG ITS 1 
RETURN 


END 

0UBR0U1 INE riECOMF' < NN » A ) 

DIMENSION A(NNyNN) ySCALfSdlV) 
COMMON/AREAl/UL < 1 1 9 y 1 1 9) 
COMMON/AREA2/irS<3 19) 

N~“-NN 

inait lol i2:e IF'SyULyAND SCAI.ES* 

DO 5 1=1 yN 
1PS<X)=I 
R0WNRM«0*0 
DO 2 J-tyN 
Ul. <IyJ)--A(ly J) 

IF<ROWNRM-“ABS(Ul.i: lyJ))) Iy2y2 
ROWNRM=ABS<Ul, (.[yJ>) 

CONTINUE 
IFIROWNRM) 3y4y3 
SCALES<I)=1»0/R0WNRM 
GO TO 5 
CALL SING(l) 

SCALES<T)-0.0 

CONTINUE 

G3uss3.al elimxnation with eortiol 


voii riJ‘5 ♦ 
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HO 1/ 

BIG^OtO 
HO tl T-Ki-N 
lP=-tF'S(3 ) 

B T ZE:==AHS < UL < I P » K ) ) 5kf>CAI t: B t T P ) 
lF<ST2L-]fciC0) 10 

10 BTG-SIZr 

ir»xpiv==i 

.1.1 CONriNUl 

TF<BTG) n^i..l2?l,'? 

K> CALL S1NG(2) 

GO rO 17 

I/, XFdBXPl V“ K> 

\A J--.i;PS<K) 

1PS<K) ■ LPSdDXPlV) 
IPfidTiXPlV)- J 
1'-=» KP~1PS<K) 

PIVOT~Ul.<KPi-N) 

KP1-=KX1 

DO 16 I -KPlyN 

IP-XPS<I) 

»n. <TP*l\)/P]VOr 
UL(lPyK)""FM 
DO 16 J™KP1?N 

UL < T P y J ) = UL ( I F- y J ) iCMiKUl. ( KP y J ) 
1.6 CONTINUE 

17 CONTTNUF 

KP~1PS<N) 

l'F<UL(KPyN) ) 19»18yiy 
10 CALL SING<2> 

19 RETURN 

END 

SUBROUriNE SOLUF(NNyByX) 
DIMENSION B(NN)<*X<NN) 

COMMON/ AREA J /UL( n9y 1 IV) 
COMMON/ARE A2/ IPS < 1 J 9 ) 

N«NN 

NF’J«N+1 

IF’™IPS<1) 

X<1>~B(IP> 

DO 2 T~2yN 
IP™IPS<I) 

IM1»=I“1 
SUM“0*0 
DO 1 J"lyIMl 

1 SUM”SUMTUL ( II • y J ) 5lcX < J ) 

2 X<I)==B<1P)“SUM 

IP“-IPS<N) / 
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C 


3 

4 


J. 


2 


3 


4 


6 

7 

0 

9 

10 


I. i 
1 o 

J. Am 


X(N)-=X<N)/UL ( [Ps-N) 

HO 4 ]PA(:K-=r»N 
1-=:NP1- IBACK 

I GOES wvl 

1P“-TPS<T) 

TPl-"=Tfl 

SUM=:0*0 

no 3 J=-IP1-N 

SUM -SUM+UL ( IP s. J ) *X < J > 

X(.T) -<X< I )-SUM)/UI.C lp!< J ) 

ri:turn 

ENH 

SUBROUT tMF- TMPKUU ( NN A k B X « » H) IT S ) 

lU: ME- NS JON A<NNyNN)yB<NN; -/XCNN) vH<. IJV) !.riX( lJ.y> 

C0MM0N/ARELA1 /Ul < 1. 1 9 » J 17 ) 

EiOUBLH PRECISION SUM 
N"NN 

EPS-=1*0E~8 
irMAX~1<f- 
XNORM -0*0 
BO 1 I--1»N 

XN0RM-AMAX1 (XNORMy ABS(X ( t ) ) > 

XFtXNORM) 3»:?»3 
BTGTTS'*- Al OOlO(fJ-’S) 

GO HO 10 

BO 9 1TER-l»rTMAX 
BO U T=-| 1.N 
S(lM--0«0 
BO 4 ,i-lKN 
SUM=SUM-f A<I» J))f!X( J) 

SUM .-B(I)-SUM 
R<I) “SUM 

CALI SOLVE (Ni-RyBX) 


BXNORM :-0*0 
BO 6 I-“»l 
r^=X(l ) 

X<1)---'X(I)+BX<I) 

DXNOE<M-^AMAXia»XNORM»ABS(X< I ) -T)) 

CONriNUfc" 

i;F’<TTfcR"l) Sy7»8 

BIOITS- - AL0G10(AMAX.1 <BXN0RM/XN0RM?rPS1 ) 
XF (BXN0RM“F:PS*XN0RM) jo? lOyV 
CONTXNUF 
CALL SING < 3) 

RCTURN 


l£NB 


SUBROUTTNF 

F'ORMATCfOXy 

E-OFfMArdOXy 


SING ( X WHY > , , 

'MATRIX wriH ZERO IJOW XN BLCOMPOSL' «/) 
'SINGULAR MATRIX XN BE: COMPOSED* iE~RO HI, VIBE. 1 


N SOLVE*' 
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12 ! FORMA r < lOXf 'NO CONOL’RGflNCr IN LMI'RUV* MAJfUX TO NF ARI Y OINGUIAR*') 
GO rO (l»2r3)!.TWHY 

j PRINT 11 

no TO 10 

? PRCNT 1? 

no TO 10 

3 PfaNT 13 

10 Kt TURN 

END 

r This rou(<i(i<“ f'fjlcii tlK' corrc^J st] ori rorrrpoL loiiSt* 

SUPROUTINF- CnF»:EN<XXi-OfEN!.N'-rO*M) 

INTEGER ALyPTs-K^y 

PTMENSION XX<210)H N(V>).VC;yy»9i»9) 

DIMENSION T ('>-S2i-9f9) » TX<r?9) 

DIMENSION L^?-.2i-9y9> 

NM-N "M 
M1~M+1 

WRITE (22^28) 
l-=0 

DO 19 AL~1 s-M 
DO 19 R~M1 yN 
l-I + l 

TX<ALyR)-XX<l) 

19 CONTINUE 

DO 23 AL”JyM 
DO 23 PT“1yM 
DO 23 R“MlrN 
DO 23 S“MlyN 

I=NM)lc ( NM#M3|! < AL ~ 1 > +NM)|« < P I - 1. ) + < R-M- 1 ) ) -f ( S-M ) U4 
KALyPTyRyS) =XX( J ) 

WRITE <22y24) Ry Sy Al-yBTy XX< I ) 

E(ALyBryRyS)=-0* 

23 CONTINUE 

28 FORMAT( ' The* t2 coeTf ic ientb oDfc<^j nfc'M rrom linear eoualion'V 

24 FORMAT< ' < ' y 1 1. y U. * ' ! t? i ' y II y 11 y ' ' yF-/X yFlG»8) 

CELE“0y 

DO 100 AL-lyM 
DO 100 Billyh 
DELE==-0* 

DO 90 R™MlyN 
DO 90 S-»M1yN 

E < AL y PT y R y S ) "0 ♦ 53CM < R y S y Al y B T ) * ^ 2 . icT < AL y DT y Ry S ) - 1 < AL y BT y S y R ) ) 

DEL E=~DELE+E < AL y BT y Ry S ) 

EONE~EONE+(2«3!tU(RySyALyBT)~Va<ySyBI yAL) ))t:TX( ALyR)*TX(BT yS) 

90 CONTINUE 

CF:L E^CELEl DELE 
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7.6 


100 

46 


77 


5.^ 


79 


30 


c 

c 

c 


6 

7 

t4 


:io 

12 


WRITE <22i-26> AL»11U 

FORMAK '' ei5!eh3i'i«?e "f orrel at j,ori cor ori for iho* waj r ■' » lEr 's# ' y 
I12r 's'y3Xf ' ~ ' y F I'J * 9 y GX? "HARKREEB' ) 

CONI iNUE 

WRnE <22y4A> EONF 

FORMAT < 'Th<? one electron ('onlr i t»«jti on ~ ^ y5XyF15»8) 

WRLTl. <22y27) CELE 

FORMAT < ' The tota#l pja'harijie correlation correction ~'ft5X)‘ 
iF15*9y5Xy M^ARTF^Fl:H' ) 

ET0T=-CELFtE:0!Jt0*b 
WRITE<22y33> FTOT 

FORMAT THE TOTAI. ENtROY WITH EXCHANGE-CORRELATION CORRECTION' y 
15Xy yF10*9y3Xy 'HARIRt-E S' ) 

WRnE<22y29) 

FORMAT ( ' NOW FOLLOW THE FA.tRWTBE EXCITATION CONTRiminONS ') 

HO 30 Al = i yM 
HO 30 Pr"1yM 
no 30 R«M I » N 
DO 30 S-~MlyN 

IF<E<ALyDTyRy 8 ) .ErCUO*) GO TO 30 
WRITE(22y31) AL yDTyRySyE< Al yPTyRyS) 

FORMATdOXy' E( ' y 1 1 y ' y ' y 1 1 y ' y ' y 1 .1 y ' y ' y I J. i- ' > -'y5XyF10»9) 

CONTINUE 

CALL ANALd ) 

RETURN 

END 


The followinsf routine does the eartial-wave 
eKChande-correlaftion rorrectj orir>4 


ana I yi3i& of 


the 


SUBROUTINE ANAL(F) 
DIMENSION E<2y2y9y9) 
WRITE <22y<!>) 

FORMAT< 20Xy'// rARTIAI. 


WAVE ANALYSIS OF EX-COR ENERGY 


WRITE (22y/) 
FORMA KlSXy'. 

1 - 



WRITE(22y141 
FORMAT < llXy 
WRITE (27 y 7) 
DO 8 I--J y2 
DO 8 J"ly2 
IF (I 4 OT 4 J) 


rAIR'y20Xy'SG'y28Xy 


GO TO 8 


PWEN'-^Oy 
DO 12 K<Ty9 
DO 12 L=^3f9 

IFtK^EQ^A^ANDtl-^FO*/)) 00 TO 10 
PWEN«PW£Nf E ( I y J y K y L ) 


' FT ' y / ) 


GO TO 12 

PWENl-E<IyJyKyL) 

CONTINUE 


Z/'y/) 
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IF <I*EO* UANri*J*EQ* 1.) GO TO 2t 
IF < I ♦ Ef) ♦ I ♦ ANIi ♦ J ♦ EO 4 2 ) 00 TO 22 
ir (l^EtU'J^ANru 04^042) 00 HO 2;^! 

L' 1 WR T TE ( 22 » 3 I ) P WEN 1 H 'WEN 

GO TO B 

:>2 PWEN™2 4*PWEN 

PWEN1=24*PWEN1 

WRITE <22i.32) PWENlirrWFN 

GO TO 8 

2,J( WRITE (22y33) F'WENJ fPWEN 

>-51 rORMAT(lIX» i-15XyFli:549nriX*n,549i-/) 

,■52 FORHAKlOXy 'ly>2s' yi;"jXyFH“i49y30XyFrJ49,/) 

FORMAT < 1 1 X » ' ?.<i2 ' ? lOX y FI 5 » 9 y IfiX y FtO" ♦ 9y / ) 
WRITE (22y40) 

40 F0RMAT<15Xy 

J . _._-'y/) 

8 CONTINUE 

RETURN 
END 
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C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

L' 

c 

n 

(; 

c 

n 

G 

G 


G 


c 

c 

G 


109 

209 


088 

n 

c 


r 


PR0Gr<AM TO RVAl.UATir IIIF F-XF’lim A1 TON VALUGS OF THF OFLRATORO 
K 2i-R ANH 1/K IJOCNO HC CI.USfFR EXPANDEn WAMF FUNCTION 
fOR BERYLLIUM IHOELFCIRON fC OERTES OF ATOMIC OYSTEMS 

INPUT IS FROM TWO DION FIl KO ANH THE TTY 

THE T1Y INPUT (ACCEPT S I ATEMENT ) 18 ENTERFD IN 10 THE 

CONTROL FILE JN THE BATCH MOBE 

THE PJI.E 1I8KJ rUAl tDAT CONTATNINO I t I ry<3 • UALUES 

IS OUTPUT FROM THE CIUSIFR FOUATIONS PRUBRAM 

THE FILE tiSK1F0R04*DAT CONTAINING THE MATRIX ELEMENTS OF 

THF OPERAIORS OUER THE HIS NUMERICAL ORBITALS IS GENERAIEU VTA 

IHE MOHTFIEB “ ZARF " PROGRAM ♦THB:'SE. VALUES ARE f HOWEVER? 

ONIY OVER THE RABIAI PARTS* ANGULAR FAC TORS ARE INIROHUCFB 
WITHIN THIS ROUTINE. 

IHE OUTPUT 18 WRITIFN ONTO THE EILF 
USKJUNFINI .OUT 

WHICH IS APPENHEB 10 ON FVERY SUBSFOUENT RUN 
OF THIS PROGRAM 


DIMENSION T<2???3J9?S19> 

DIMENSION A0TNT(9?9) 

INTEGER AL?Br?RyS?R(9) 

LOGICAL HI?LO 

DOUBLE PREC I S t ON PROP 1 ( i ) » CASE 

DATA PROP!/' R ' 2 R 'f' 1 / R '/ 


OPEN < UNIT-4 y DEV 1 CE"- ' DSK ' ) 

0PEN<UNIT-=5yDEVICE'= 'DSK' yFIl E=“ 
0rEN<UNIT~6yriEVlCE"-''DSK' yETLE - 


" IVAL.DAT' ) 

' ONEINI . OU r ' y ACCESS-^ ' APPEND ^ ) 


A TITLE OP 10 CHARACIER LENGTH IS REIAD IN 


ACCEPT 109yCASE 
FORMAT <A10) 

FORMAT< ////////%/' " 

1 ' Jl! CASE? 'y AlOy ' 

READ<5yJl«) ( < < ( ( (ALyDT yRyS) ,S-3y9> yR-3y9) yBT-1 y2> yAL"l f 2) 

DO 9999 IPR“-ly3 . . 

READ< 4y3|f) ( (AOINT < I. y J> y T~1 1'9) y J=»l y9) 

WRJTE(6y888)PR0Pl< fPR) ^ __ 

FORMAT (///' - l~ LI PROPERTY t yAlOy =» > 

angular FACTORS FOR AO INTEGRALS ARE INTRODUCED 


DO 1 I«ly9 
R(I)-I 
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C0N7 INUr 


'•) 


400 

u 

c 

c 

n 


3 

V, 

06 

c; 

c 

67 


10 


200 


c 


on 2 

DO 2 J=-=t)-9 

lf'<K<T) ♦EOJ\U) )00 TO 2 
AOtNK J y J)-->0* 

CON 7 iNur: 

WIC I rr. ( 6 y 400 ) < ( AO I N7 < r . J ) y J ••• J y 9 ) y J! ' I y 9 > 

FOF?MAr (//' 3{c5}c3K)fC;|(>fi L -I.L INI OVf.'K* A0"0 >!(!|:)|C3f:3l{>|£ '//< 

IF (TPRyNE:* t)00 to 06 

THF: "i* VALOES AkC OONUr.RtTD 70 flir CORkEni'-ONDtNfi 
CT Coei-'F tCXENTS 


HO ,“7 Al."ty2 
DO 3 Dr-Jy2 
DO 3 R~-=3y9 
DO 3 S-=3y9 

7 (ALyBTyk*S)=^T(Al yDTyRyS)/4,0 
CONTINUE 

CONTINUE 

CALCULATION FDR RF FFRENCE S7A7F OTNGLE DFlt 


CONTINUE 

SUM-0 

DO 10 I=-<ly2 
SUM-SUM+2#A0INT(I» 


BASE=SUM 

WRI!E<6y200) SUM 

FORMAT</// " rXFlCIATtON OAI.UE=-- 'yri2.by 
I- FOR K-EKEfv-l.'NCE B DEI') 

KVALUA7I0N OF NtiRMAL ISATION FACTOR 


ANORM- 1 . 

I HE FACTOF<: RRE 
DO 20 AL-ly? 

DO 20 BT-=.<ly2 
DO 20 R^3y9 
DO 20 S"3y9 
F'RE ~6 * 


ACCOUNTS FOR GF'JNG 


'y 9 ri 0 » 6 O 
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J F ( Al » EQ . BT » OK JC ♦ T 0 . S > l-’KL <l , 

IF ( AL ♦ EQ , BT ♦ ANB . K J- 0 * H ) 1 . 

ANGRM^ANOKM-i f < Al . s- B F i- K ? S ) * I ( Al s- B f » K y S ) )iif -Rl. 

:>0 CONTINUE 

C EVALUATE THE CONTIiTMJTJ CMS FROM OfUFI? ULIEKhJNANIO 

BO iOO Al 
BO 100 BT-jy:> 

BO too K~ j;yy 

BO JOO G~3yy 

BO 99 lAL- J 
BO 99 iBr-=jy:> 

BO 99 IK~3y9 
BO 99 J S“3 y 9 

rA(>r(Al yBIyRyS)*? (lALylBIv [RylS) 


C COMF’ARE THE BFTS5 t ALyBTyF<yB > ANB i XALy ] B1 y I K'y XS > 

LO«<fi.HQ* [R4ANI.uSa;0y li>).OR*<K*EQ*]b\ANl[i»lUEay JR) 
HI-=<AL*rQ»XAf.4AND*Fir4FQ4 DBO^OR^ <AI. ♦EOyTHT^ANByBI A 0* JAL) 
jr< ♦NOT* (HI »0R4I.0) )G0 TO 91 
jr<LOtANB«Hi)oo ro 5.1 
IF <1.0)00 TO 41 

C NON- IBENTITY AMONG t RyS F ANB I ilRylS > 

C CHECK IF ONLY ONb JO FFLRENCE EX) BIG ANB BRANCH 

IF <R*EtnXR) GO TO 

ir (S^EO^IS) GO TO 33 

IF (RyEQylS) GO ro 34 

IF (S*EO*XR) 00 TO 35 

GO TO 91 
r? M=B 

N”I&‘ 

GO TO 80 

33 M=F? 

N=«tR 

GO TO 80 

34 M«-XS 
N-R 

GO TO 80 

35 M=-XR 
N~8 

GO TO 80 
41 CONTINUE 

C NON~IBENI riY AMONG < ALyB > ANB < 

C CHECK IF ONLY ONF BlI'FERF.NCr: EXISTS ANB BRANCH 
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IF 

<AI,. 

.E.O. 

J.AL ) 

00 

ro 

49 

IF 

<BT 

a;o* 

FBT) 

GO 

ro 

AS 

IF 

(Al. 

*E:.0* 

TBD 

GO 

10 

44 

IF 

<Bl 

4 EG* 

J AL ) 

GO 

TO 

40 

GO 

TO 

91 






N~rPT 

GO ui no 

4.^ M--AI 

GO ro GO 
44 

N-AL 

GO TO «<) 

41) M=l'AI 

NS-B7 

GO TO »0 

0 INIiEX BF7H nurNIlOAI. 

n.l CONrJNUlr 

CONIR-PAGK+AOlNli <|\>R)-fAOINr (G kOV' AOJN l <Al yAD-AOtNrd^r i-BF ) 
PRr -6 ♦ 

IF < Al. ♦ 1= 0 ♦ B T ♦ OR ♦ Rf EQ * 0 ) I 'FcC- -2 ♦ 

I F ( AL ♦ EO ♦ B r ♦ ANLi ♦ R « FQ ♦ 0 ) fl?!-'- 1 ♦ 

SUM-. SUM f r AC5KCON I R>KH-:£. 

GO TO 99 

C ONfc mrFFRf.NCL- IN INBITX SET 

00 SUM."BUM I F- AOjKAOIN r ( M 1- N ) 

GO ro 99 

0 

0 THE, INBF-X SIHS filFI CR BY M0F<F1 I HAN ONI.. 

0 

91 CONT ENUE; 

0 

99 CONJ 1 NUEE 

0 

J,00 CONffNUE 


C 


;roo 


NORMALISE TIIIE REOUl T ANB OUrPU'l 


SUM-' BUM/AMORM 

WFnTL(6y.r00) BUM 

FORMAT (/// ' e.Xf'EOTAIlON 


VALUE:..’' 
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1 " ALL CONFN?) INCLUDKm') 

9999 CDNIINUF 
S f OP 
F'ND 



